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I.  INTRODUCTION 


The  phenomena  of  fatigue  crack  closure  was  originally  investigated  by  Elber  (1)  in 
the  early  1970's.  His  discovery  that  the  crack  remains  closed  during  the  initial  portion  of  the 
load  cycle  led  to  the  definition  of  an  "effective"  load  range,  as  opposed  to  the  applied  load 
range,  for  predicting  fatigue  crack  growth  behavior.  In  effect,  the  applied  load  range  is 
reduced  by  the  closure  load  level,  thus  defining  an  effective  stress  intensity  factor  range. 
Elber  later  utilized  this  closure  information  to  develop  representative  models  for  predicting 
fatigue  crack  growth.  According  to  Paris  (2),  understanding  of  this  observation  has  the 
potential  for  contributing  significant  insight  into  explaining  observed  crack  growth  behavior. 
This  statement  by  Paris  places  significant  emphasis  on  conducting  a  thorough  and  complete 
investigation  into  all  aspects  of  closure. 

With  an  increased  understanding  of  the  closure  phenomena,  it  has  become  evident 
that  a  consistent  set  of  definitions  of  terms  such  as,  closure  load  level,  crack  opening 
displacement  profiles  and  the  inherent  effective  stress  intensity  factor  has  not  been 
established  because  of  the  complex  nature  of  closure.  This  is  a  major  problem  area  which 
has  caused  inconsistencies  in  generating  experimental  data  and  evaluation  of  results.  Several 
technical  reviews  on  closure,  such  as  those  by  Baneijee  (3),  and  Suresh  and  Ritchie  (4),  have 
identified  the  need  to  establish  a  cohesive  closure  data  base  utilizing  a  consistent  set  of 
definitions. 

The  selection  of  surface  flaws  (part-thru-cracks)  for  this  investigation  of  closure 
effects  on  fatigue  crack  growth  rate  was  based  on  its  occurrence  in  many  engineering 
applications.  Although  a  common  crack  type,  the  amount  of  crack  closure  information 
available  on  surface  flaws  in  comparison  to  thru-thickness  cracks  is  relatively  small.  This  is 
due  to  the  three-dimensional  complexities  of  surface  flaw  crack  opening  displacement 
(COD)  patterns,  the  difficulty  of  measuring  COD  experimentally,  and  the  lack  of  usable 
stress  intensity  factor  solutions  around  the  crack  boundary. 
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The  three-dimensional  aspects  of  closure  in  surface  flaws  impose  major  measurement 
restrictions  in  metals  since  only  the  crack  face  is  directly  visible.  Typical  closure  and  fatigue 
crack  growth  (FCG)  measurement  techniques  used  in  the  past  have  not  provided  real-time 
information,  but  required  data  to  be  analyzed  at  the  end  of  a  test  (such  as  with  heat  tinting 
and  crack  replication)  or  utilized  analytical  representations  of  crack  opening  displacements 
(such  with  compliance  gauges).  The  use  of  a  transparent  material  in  this  investigation  solves 
many  of  these  problems. 
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II.  BACKGROUND 


The  following  sections  present  a  background  review  of  prior  experimental  and 
analytical  investigations  of  closure,  surface  flaws,  and  the  use  of  closure  for  predicting 
fatigue  crack  growth. 


Fatigue  Crack  Growth  and  Closure 

The  fracture  mechanics  concept  of  plasticity-induced  fatigue  crack  closure  was 
originally  introduced  by  W.  Fiber  (1)  some  20  years  ago.  In  essence.  Fiber  found  that  cracks 
do  not  open  elastically  under  applied  load,  but  are  held  closed  by  some  mechanism  which  he 
termed  crack  closure.  He  determined  that  the  closure  load  of  an  aluminum  panel  with  a 
cracked  hole  could  be  as  much  as  50%  of  the  applied  loading  amplitude,  leaving  this  range 
of  the  load  cycle  ineffective  for  growing  the  crack.  Hie  closure  concept  initially  gained 
widespread  popularity  in  the  fracture  mechanics  community  because  of  its  potential  to 
predict  fatigue  crack  retardation.  After  years  of  research,  interest  in  the  closure  concept 
began  to  waver  in  the  late  70's  because  of  a  lack  of  progress  in  sorting  out  the  many 
variables. 

The  importance  of  closure  for  predicting  fatigue  crack  growth  has  gained  renewed 
enthusiasm  in  recent  years  as  a  result  of  new  frndings.  This  rekindling  of  interest  has  come 
about  because  of  a  better  understanding  of  the  implications  of  closure  and  because  of  an 
expansion  of  theoretical  closure  models.  As  was  stated  by  Paris  (2),  closure  is  "the  most 
central  question  still  to  be  resolved"  in  the  area  of  fatigue  crack  growth.  This  is  a  strong 
statement  coming  from  a  scientist  who  has  made  significant  contributions  to  validate  the  use 
of  the  stress  intensity  factor  (K)  in  crack  growth  analysis.  In  reference  to  closure,  Paris 
describes  it  as  "simply  good  luck"  that  the  usefulness  of  K  was  not  eliminated  because  the 
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closure  phenomena  had  been  overlooked.  In  a  recent  paper  by  I .  ris,  "Twenty  Years  of 
Reflection  on  Questions  involving  Fatigue  Crack  Growth,"  he  concluded: 


It  seems  that  crack  closure  reappears  as  a  key  factor  and  that  methods  must  be 
found  to  analyze  it  or  a  new  fundamental  way  of  avoiding  it  in  analysis  needs  to 
be  formulated.  It  is  viewed  by  this  author  as  raising  the  most  central  question 
still  to  be  resolved.  It  seems  difficult  to  avoid  concluding  that  crack  closure  is  a 
key  physical  phenomenon  in  the  fatigue  cracking  process. 


The  relationship  for  fatigue  crack  growth  (FCG)  rate  was  given  by  Paris  as: 

da/dN  =  C  (AK)"  (2.1) 

and  modified  by  Elber  to  include  the  closure  parameter  U  (discussed  later),  such  that: 


da/dN  =  C  (U  AK)“ 


(2.2) 


=  C  (AKeff)« 


(2.3) 


where  C  and  n  are  material  constants,  and  AK  the  cyclic  change  in  stress  intensity  factor. 


Surface  Cracks 

A  great  deal  of  fracture  mechanics  technology  is  based  on  experimental  evidence 
using  two-dimensional  through-the-  thickness  crack  specimens  and  stress  intensity  factor 
solutions.  Although  surface  cracks  (Figure  2.1)  represent  the  most  common  natural  defect 
shape  for  crack  initiation  and  propagation,  analytical  modeling  complexities  have  caused 
investigators  to  work  with  the  more  manageable  thru-crack  configuration.  The  complexities 
of  modeling  surface  crack  stress  intensities  arise  from  its  three-dimensional  characteristics. 
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Figure  2.1  Surface  flawed  fatigue  crack  specimen  loaded  in  bending 


such  as  constantly  changing  stress  field  and  stress  intensity  factor  around  the  crack  tip 
perimeter,  changing  crack  aspect  ratio  (a/c)  during  crack  growth,  and  a  lack  of  available 
closed-form  stress  intensity  factor  solutions. 

The  lack  of  a  useful  surface  flaw  stress  intensity  factor  solution  has  been  a  major 
reason  for  the  scarcity  of  experimental  fatigue  crack  growth  data.  In  recent  years,  this 
condition  has  improved  with  the  introduction  of  stress  intensity  factor  solutions  for  surface 
flaws  in  cracked  plates.  One  of  the  more  popular  solutions  was  presented  by  Newman  and 
Raju  (S),  who  utilized  finite  element  calculations  to  develop  empirical  relationships  for 
numerous  crack  configurations.  These  correlation’s  were  developed  for  both  tension  and 
bending  —  using  parameters  such  as  crack  depth  (a)  and  length  (c),  plate  thickness  (t)  and 
width  (b),  and  angular  location  around  the  crack  perimeter  (^). 

Orange's  (6)  report  to  the  ASTM  Task  Group  E24.01.05  on  "Part-Through  Crack 
Testing"  also  identified  a  general  lack  of  experimental  data  on  surface  flaws.  His 
conclusions  were  based  on  a  survey  questionnaire  of  the  major  aircraft  companies,  NASA, 
and  the  Atomic  Energy  Commission.  In  a  November  1984  telephone  conversation,  he 
reiterated  that  besides  the  introduction  of  stress  intensity  factor  solutions,  little  had  been 
accomplished  since  1975  to  change  his  original  recommendations. 


Closure  Models 

The  physical  representation  of  closure  is  based  on  five  fundamental  models:  the 
original  plasticity-induced  closure  model,  and  four  recent  models  which  have  developed 
from  microscopic  descriptions  of  crack  surfaces.  These  five  closure  models  constituted  a 
growing  awareness  of  closure  complexities,  and  provide  a  brief  insight  as  to  why  the  fracture 
mechanics  community  has  struggled  to  quantify  closure  relationships  to  crack  growth.  The 
four  recent  models  have  been  shown  to  be  of  particular  importance  in  the  near-threshold 
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region,  and  are  identified  as  oxide-induced,  roughness-induced,  viscous  fluid-induced  and 
phase  transformations-induced.  Each  model  type  is  represented  schematically  in  Figure  2.2 
and  briefly  described  in  the  following  paragraphs. 


Plasticity-Induced  Closure 

The  Elber  plasticity-induced  closure  model  (Figure  2.2a)  is  characterized  by  the 
development  of  a  plastic  zone  around  the  crack  tip,  and  subsequent  growth  of  the  fatigue 
crack  thru  this  plastic  zone.  This  zone  is  described  as  a  residual  "wake”  of  plastically 
deformed  material  around  the  crack  face  (Figure  2.3)  which  produces  the  crack  closure 
effect  The  approximate  size  of  the  plastic  zone  at  the  crack  tip  was  originally  modeled  by 
Irwin  (7)  as  circular,  and  later  by  Dugdale  (8)  as  a  narrow  strip.  The  plastic  zone  radius  (rp) 
was  represented  by  Irwin  for  plane  stress  as: 

tp-[K/aysl2/2n  (2.4) 


and  for  plane  strain  as: 


rp=[K/Oys]2/6n  (2.5) 

where  K  is  the  mode  I  stress  intensity  foctor  and  Oy^  the  material  yield  strength. 

When  a  three-dimensional  stress  fleld  is  modeled  by  using  the  Tresca  or  von  Mises 
relationships,  the  crack  tip  plastic  zone  is  more  complex  than  a  circular  or  narrow  strip, 
forming  a  kidney  shape.  Banks  and  Garlick  (9)  present  a  detailed  discussion  of  three- 
dimensional  plastic  zones  for  both  plane  stress  and  plane  strain  in  isotropic  and  anisotropic 
materials.  Since  K  is  a  direct  function  of  stress.  Figure  2.4  reflects  the  general  shape  of  a 
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2.2  Schematic  Illustration  of  the  various  mechanisms  of 
fatlque  crack  closure  (Ref.  4) 
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thni-the-thickness  yield  surface  where  the  most  yielding  occurs  at  the  free  surface  (plane 
stress,  O2  =  0),  while  a  smaller  yield  zone  is  present  at  the  middle  of  the  specimen  (plane 
strain,  62  =  0). 

As  an  expansion  to  his  original  closure  model,  Elber  (11)  presented  an  analytical 
relationship  to  account  for  mean  stress  (quantified  by  the  stress  ratio,  R)  influence  on  closure 
loads  in  aluminum.  This  approach  led  Elber  to  deflne  an  effective  stress  intensity  factor 
(AKet^  which  is  related  to  fatigue  crack  growth  rate  (da/dN)  by  an  effective  stress  intensity 
factor  closure  parameter  (U)  and  the  applied  cyclic  stress  intensity  factor  (AK)  such  that: 

da/dN  =  C(U  AK)“  =  QAKgff)"  (2.6) 

where: 

AKeff=UAK  =  K^ax.Kci  (2.7) 

and  C  and  n  are  empirical  material  constants  determined  from  the  Paris  da/dN  versus  AK 
fatigue  crack  growth  curve.  The  effective  stress  intensity  factor  closure  parameter  is  deHned 
as: 


U  _  _..^max.'  ^1 


^max  *  '^min 


Here  ^min  ^1  ^p  (figure  2.5)  are  the  maximum,  minimum  and  closure 

stress  intensity  factors,  respectively.  In  effect,  if  is  larger  than  the  usable  AK  for 
growing  the  crack  is  reduced,  thus  reducing  the  fatigue  crack  growth  rate. 

The  terms  opening  and  closure  stress  intensity  factor  (K^p  or  K^i,  respectively)  are 
often  used  interchangeably  in  the  literature,  although  there  are  important  differences  between 
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Figure  2.5  Schematic  of  applied  and  effective  stress 
Intensity  (fijC  ,,)  load  ranges  (Ref.  3) 


the  two.  Kop  is  used  to  define  closure  when  load  is  applied  to  the  specimen  (increased 
loading)  whereas  Kd  is  used  when  the  specimen  is  unloaded  (decreased  loading).  Although 
the  closure  stress  intensity  factor  (K^)  and  opening  stress  intensity  factor  (K^p)  have 
different  definitions,  they  have  been  reported  to  be  approximately  equal  (3).  As  a  general 
guideline,  most  researchers  agree  that  the  opening  load  should  be  used  to  determine  AK^ff 
for  fatigue  crack  growth  calculations  instead  of  the  closure  load.  This  is  based  on  the  crack 
tip  opening  (versus  closing)  being  a  primary  driving  mechanism  for  fatigue  crack  growth. 
All  values  reported  in  this  investigation  are  opening  loads,  unless  specifically  stated  to  be  a 
"closure"  load  for  comparison  purposes.  Otherwise,  the  term  closure  and  opening  load  may 
be  used  interchangeably. 

For  2024-T3  aluminum,  Elber  determined  that  U  was  related  to  the  stress  ratio  R 
(where  R  =  Pmin^max) 


U=:0.5  +  0.4R  (2.9) 

for  -0.1  <  R  <  0.7.  Since  its  original  introduction,  the  coefficients  of  this  formulation  for 
aluminum  have  changed,  but  the  general  approach  for  modeling  the  effective  stress  intensity 
factor  is  being  used  today. 

The  generally-agreed-upon  experimental  definition  of  closure  can  be  seen  in  Figure 
2.6,  where  applied  load  is  plotted  against  crack  opening  displacement  (V)  or  offoet 
displacement  (AV).  The  load  versus  displacement  measurement  is  typically  generated  using 
a  clip  gauge  attached  to  the  crack  face,  or  more  recently  by  using  an  optical  interferometry 
S3rstem  (12, 13, 14).  The  offset  displacement  definition  of  closure  is  considered  more 
sensitive  to  data  deviations  with  experimenters  using  the  interferometry  systems  (IS).  The 
effective  load  range  associated  with  equation  2.7  is  defined  by  (Pmax'  Pci)  (Pmax "  ^l)- 
Two  examples  of  the  usefulness  of  crack  closure  can  be  seen  by  reviewing  Figures 
2.7a  and  2.7b.  Figure  2.7a  shows  the  direct  relationship  of  stress  ratio  (R)  to  fatigue  crack 
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Figure  2.7  Effect  of  load  ratio  (R)  on  fatigue  crack  growth  rate  da/dN  (a) 
without  closure  and  (b)  with  closure  load  Influences  (Ref.  3) 


growth  rate  (FCGR)  -  as  R  increases,  FCGR  increases.  By  utilizing  closure  stress  intensity 
factor  (Kcl)  ^  calculate  an  effective  stress  intensity  factor  (AK^ff),  instead  of  an  applied  AK 
used  in  Figure  2.7a,  the  data  collapse  around  one  FCGR  curve  (Figure  2.7b).  Similar  effects 
are  seen  in  Figure  2.8  where  comparisons  are  made  between  AK  and  AK^ff  for  notched  steel 
specimens  in  the  near-threshold  region  for  short  and  long  cracks.  These  examples  are  an 
indication  that  closure  is  a  fundamental  mechanism  for  controlling  FCGR,  that  effective 
stress  intensity  factor  normalizes  the  FCGR  data,  and  that  applied  AK  does  not  totally  reflect 
the  crack  growth  driving  mechanism. 


Oxide-Induced  Closure 

The  oxide-induced  closure  mechanism  (Figure  2.2b)  arises  when  oxide  deposits 
collect  on  the  crack  internal  surfaces  and  wedge  the  crack  open  during  unloading.  This 
mechanism  has  also  been  used  to  describe  the  influence  of  a  moist  environment  on  fatigue 
crack  growth  rate  (FCGR).  The  model  is  particularly  useful  in  the  near-threshold  region 
where  closure  stress  intensity  factor  (K^^j)  levels  are  comparable  to  the  minimum  stress 
intensity  factor  (Kmin)  levels. 


Roughness-Induced  Closure 

The  roughness-induced  closure  mechanism  (Figure  2.2c)  arises  when  the  crack 
ffacture  surface  asperity  is  of  a  comparable  size  to  crack  opening  displacement  (COD)  levels. 
This  mechanism  can  also  be  present  when  significant  mode  II  displacement  exists.  Like 
oxide-induced  closure,  there  is  a  wedging  action  which  keeps  the  crack  surfaces  apart  during 
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unloading,  causing  K^j  to  be  greater  than  K'"'”.  Again,  this  characteristic  is  dominant  in  the 
near-threshold  stress  intensity  factor  region.  The  process  effectively  lowers  AK  and  thus 
reduces  FCGR. 


Viscous  Fluid-Induced  Closure 

The  presence  of  a  viscous  fluid  (Figure  2.2d)  in  the  crack  can  produce  a  wedging 
effect  due  to  hydrodynamic  pressure  when  the  crack  opens  and  closes.  This  effect 
accelerates  the  fatigue  crack  growth  rate  in  the  near-  threshold  region,  but  not  in  the  higher 
_K  regions  where  an  oxide-induced  mechanisms  (caused  by  the  moist  environment) 
dominate  the  closure  loads  effect. 


Phase  Transformation-Induced  Closure 

Materials  which  undergo  a  stress/strain  induced  phase  transformation  during  cyclic 
loading  produce  a  closure  mechanism  similar  to  plasticity-induced  closure  (Figure  2.2e). 
The  crack  tip  plastic  zone  produces  a  metallurgical  phase  transformation  which  forces  the 
yielded  material  into  a  compression  state  to  compensate  for  the  unchanged  elastic  material. 
As  the  crack  grows  through  this  transformation  zone,  it  experiences  a  closure  load  effect 
which  reduces  the  nominal  stress  intensity  factor  range.  Although  little  experimental  data 
are  available,  increased  temperature  is  expected  to  enhance  phase-transformation-induced 
closure. 
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Crack  Opening  Displacement  Measurements 


The  use  of  Newton  interferometric  measurement  techniques  for  defining  crack 
opening  displacement  (COD)  patterns  in  a  transparent  polymer  specimen  provides  a  unique 
experimental  capability.  The  technique  employs  the  principle  of  interfering  light  paths 
reflecting  off  two  surfaces  which  are  typically  separated  by  an  air  gap  (Figure  2.9). 

Packman  (16)  described  this  interferometry  technique  and  its  usefulness  for  measuring  stress 
intensity  factors,  fracture  of  materials  and  stress  corrosion.  He  used  the  term  crack  opening 
interferometry  (COI)  for  measuring  crack  opening  displacements  (COD)  in  order  to 
distinguish  it  from  other  experimental  techniques  such  as  clip  gauges,  strain  gauges,  etc.  The 
COI  technique  is  particularly  useful  for  crack  opening  measurements  since  it  can  distinguish 
very  small  displacements  (to  within  a  fraction  of  a  wavelength).  By  counting  the  number  of 
fringe  patterns  (t))  and  knowing  the  light  source  wavelength  (X),  the  air  gap  thickness 
produced  by  COI  can  be  calculated  as  (17); 

COI  =  2d  =  (2ti  - 1)  X  /  4  (2.10) 

for  the  destructive  fringes  (dark  lines)  with  a  normal  incidence  light  source  (where  r)  = 

1,2,...).  This  technique  has  been  successfully  applied  in  a  number  of  investigations,  and  will 
be  discussed  in  more  detail  in  the  body  of  this  report. 

Closure  Observations 

An  investigation  by  Fleck,  Smith  and  Smith  (18)  used  a  crack  mounted  displacement 
gauge  and  back  face  strain  gauges  to  measure  surfece  crack  opening  displacements  on 
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generating  fringe  patterns  from  crack  opening 
displacement  (Ref.  19) 


compact  tension  specimen  of  BS4360  B50  steel  (24%  elongation,  55  ksi  yield  strength).  A 
new  "push-rod  gauge"  was  developed  to  measure  cracic  internal  surface  displacements  in  the 
specimen  plane  strain  region.  The  investigators  found  the  push-rod  gauge  to  be  accurate,  and 
concluded  "that  a  greater  degree  of  crack  closure  occurs  at  the  surfece  than  the  bulk  plane 
strain  regions  of  the  specimen.”  Tbe  closure  parameter  U,  defined  as  (Kjq^- 

was  calculated  to  be  0.85  in  the  plane  strain  region,  and  0.75  at  the  plane  stress 
region  for  the  surface  crack. 

Pitoniak  (19)  used  PMMA  and  a  monochromatic  light  source  to  measure  crack 
opening  displacement  of  thru-the-thickness  cracks  in  compact  tension  specimens.  His  crack 
opening  displacement  (Figure  2.10)  measurements  revealed  that  at  zero  load  the  thru-crack 
was  closed  at  the  specimen  surface  (plane  stress  regions),  but  was  always  open  in  the  crack 
center  (plane  strain  regions).  That  ts,  at  zero  lc»d  the  COD  plots  show  the  crack  internal 
surface  to  be  open  approximately  1.3E-4  inches  (3.3E-3  mm)  at  the  max  displacement  point, 
while  closed  at  the  crack  face. 

Further  interpretation  made  here  of  the  closure  load  measurements  in  the  Pitoniak 
report  reveal  that  approximately  17%  of  the  maximum  applied  load  (Pop/Pmax  required  to 
open  the  internal  crack  tip  plane  strain  region,  30%  to  open  the  crack  face  middle  and  47%  to 
open  the  crack  tip  plane  stress  region  at  the  crack  face  (compared  to  50%  for  2024 
aluminum).  It  is  also  noted  that  the  point  of  maximum  displacement  is  located 
approximately  30%  of  the  crack  length  from  the  crack  tip.  These  frndings  are  significant 
since  they  not  only  show  an  internal  displacement  under  zero  load,  but  also  because  the 
measured  closure  load  ratios  (Pop^max)  similar  in  magnitude  to  those  seen  in  metals 
even  though  the  plastic  zone  size  is  considered  to  be  small. 

The  COD  pattern  observations  by  Pitoniak,  et  al.,  (20)  using  PMMA  material 
were  later  observed  by  Pelloux,  Fatal  and  McGee  (21)  using  fractc  graphic  analysis  in  metals. 
These  findings  are  shown  in  Figure  2.11,  where  the  mid-crack  tip  region  (plane  strain  region) 
is  shown  to  be  open  under  zero  load  while  closed  at  the  specimen  free-surface  (plane  stress 
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Figure  2.11  Scheioatlc  o£  the  closure  contact  areas  o€  thru> thickness 
crack  opening  displacement  patterns  using  fractographlc 
measurements  In  metals  (Ref.  21) 


region).  With  load  application,  the  crack  opens  outward  from  the  internal  crack  surface 
toward  the  crack  free  surface.  For  the  crack  surface  to  be  fully  open  required  an  applied  to 
maximum  load  ratio  greater  than  0.5  for  2124-T351  aluminum  alloy. 

Analytical  results  of  Newman  (22)  showed  the  closure  ratio  (K^i/Kq^^x)  to  be  in  the 
range  of  0.25  to  0.35  for  plane  strain,  and  0.5  to  0.6  for  plane  stress.  Analysis  by  Fleck  (23), 
using  an  elastic-perfectly  plastic  two-  dimensional  finite  element  analysis  developed  by 
Newman  (24),  determined  that  closure  load  behavior  is  insensitive  to  material  property 
changes  of  Oyg/E  and  v  (where  Oyg  is  the  yield  stress,  E  is  Youngs  modulus  and  v  Poisson's 
ratio).  When  closure  ratio  was  plotted  against  Oyg)^  at  the  plane  stress  and  plane 

strain  regions  of  the  crack,  definite  closure  load  trend  differences  were  noted.  The  plane 
strain  region  showed  a  decreasing  closure  ratio  with  increasing  load  ratio  (orjjjax^ayg). 
However,  the  plane  stress  region  showed  little  difference  for  various  load  ratios.  The  closure 
ratio  levels  calculated  by  Fleck  were  similar  to  those  calculated  by  Newman. 

Fleck  found  that  the  closure  parameter  (U)  values  for  thru-crack  measurements  were 
similar  to  those  for  surface  flaws  (0.84  for  plane  strain  and  0.74  for  plane  stress).  These 
values  compare  quite  favorably  with  the  PMMA  measurements  of  approximately  0.83  and 
0.70  for  the  plane  strain  and  plane  stress  regions,  respectively,  reported  by  Pitoniak,  et  al., 
(18).  Ray  and  Grandt  report  similar  results  for  thru-thickness  compact  tension  specimens 
(25)  and  surface  flaws  (26). 

Again,  it  is  of  interest  to  note  that  although  the  plastic  zone  size  for  PMMA  is 
reported  to  be  small,  the  closure  parameter  (U)  is  of  the  same  relative  magnitude  as  values 
reported  for  a  ductile  steel  (Fleck,  Ref.  18),  2024  aluminum  (Elber,  Ref.  11)  and  for  analysis 
(Ref.  22).  These  results  are  summarized  in  Table  2.1. 
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Table  2.1  Summary  of  Closure  Ratio  and  Closure  Parameter  (U)  Values 


Thru-Crack  Specimens 

Kr.l/Kmay 

U 

Elber  —  2024  Aluminum 

0.50 

0.50 

Fleck-  Mild  Steel 

0.27 

0.74 

Newman  -  Analytical 

0.50-0.60 

0.50-0.40 

Pitoniak  -  PMMA 

0.47 

0.53 

Ray/Grandt  -  PMMA 

0.38 

0.69 

Part-thru  Surface  Cracks 

Fleck  —  Mild  Steel 

0.25 

0.76 

Ray/Grandt  -  PMMA 

0.18-0.30 

0.91-0.78 

Closure  Problem  Areas 


Detailed  literature  reviews  of  closure  research  have  identified  important  questions 
concerning  the  state-of-the-  art  in  understanding  and  using  closure  relationships  to  predict 
fatigue  crack  growth  (FCG).  In  general,  these  reviews  have  identified  the  need  for  more 
fundamental  research  into  this  area. 

Banerjee  (3)  makes  the  following  observations  in  his  review  of  crack  closure: 

1.  Even  though  the  use  of  an  effective  stress  intensity  factor  (AKgff)  has 
normalized  fatigue  crack  growth  rate  data  over  a  wide  range  of  ^in 
and  R  values,  the  actual  closure  mechanisms  are  not  well  understood. 

2.  The  dependence  of  closure  stress  intensity  factor  (Kg|)  on  K^ax»  ^in 
-R  for  constant  amplitude  load  test  is  not  clearly  established.  As  an  example, 
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K^l  has  been  observed  to  increase,  remain  constant  or  even  decrease  with 
increased 

3.  Effects  of  cyclic  load  history  on  closure  have  not  been  established.  This 
includes  the  effects  of  preciacking,  variable  amplitude  loading,  frequency 
changes,  etc. 

4.  Even  though  crack  stress  intensity  varies  along  the  crack  perimeter  (and  is 
thus  di^rent  at  the  interior  and  at  the  exterior  of  the  crack  •  reference  Figure 
2.4),  the  influences  of  plane  stress  and  plane  strain  conditions  are  not  clearly 
defined,  well  understood  or  typically  accounted  for  in  the  literature. 

5.  There  is  little  distinction  in  the  literature  concerning  free>surface  and  mid¬ 
thickness  closure  behaviors.  The  lack  of  a  consistent  definition  has  led  to 
confusion  about  the  actual  crack  tip  stress  intensity. 

6.  The  extent  of  closure  and  residual  displacement  produced  by  closure  is 
rarely  repotted. 

7.  Plasticity  induced  closure  depends  on  the  formation  of  plastic  zones  and 
inelastic  deformations  in  the  vicinity  of  the  crack  tip. 

8.  The  fundamental  closure  models  (plasticity  versus  asperity,  oxide-induced, 
etc.)  reflect  observed  phenomena  in  separate  stress  intensity  factor  (K)  regions  - 
plasticity  in  the  high  K  region,  asperity  and  oxide  in  the  low  K  regions. 

9.  Asperity  and  oxide  induced  closure  operate  only  in  the  near  threshold 
region. 


Conclusions  by  Surech  and  Ritchie  (4)  follow  the  same  general  theme.  Their 
observations  are  summarized  here: 

1.  Understanding  the  role  of  closure  is  essential  to  understanding  the  behavior 
of  short  cracks,  variable  amplitude  loading,  single  overloads  and  block  loading 
above  and  below  the  threshold.. 
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2.  Crack  length  is  critical  to  evaluate  closure  because  of  wake  effect 
load  history  on  the  wake/closute  relationship. 


effect  of 
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III.  PROBLEM  STATEMENT 


A  great  deal  of  activity  has  gone  into  understanding  crack  closure  and  opening 
displacement  characteristics  since  Elber's  initial  description  of  the  closure  phenomenon. 
Most  of  these  investigations  have  used  thru-crack  specimens  because  of  the  availability  of 
stress  intensity  factor  solutions  and  the  relative  ease  of  making  closure  measurements.  In 
comparison,  relatively  few  closure  investigations  have  been  conducted  with  surface  cracks. 

The  inability  to  measure  three-dimensional  crack  opening  displacements  in  metals 
>een  the  major  deterrent  in  conducting  surface  flaw  investigations.  The  fundamental 
data  required  to  address  the  influence  of  plane  stress  and  plane  strain  effects  on  closure  are 
also  lacking  for  the  same  reasons. 

Many  of  the  three-dimensional  measurement  difficulties  associated  with  metals  can 
be  resolved  by  using  an  optically  transparent  polymer  to  measure  surface  flaw  crack  growth. 
By  employing  interferometry  techniques  to  accurately  measure  internal  crack  opening 
displacement  (COD)  under  various  opening  loads,  the  influence  of  crack  closure  on  fatigue 
crack  growth  can  be  established.  The  key  to  this  effort  is  being  able  to  directly  measure 
three-dimensional  COD,  closure  and  crack  growth  simultaneously  using  Newton 
interferometry  techniques  in  a  transparent  polymer. 

Objective 

The  objective  of  the  current  research  is  to  determine  the  influence  of  closure  on  the 
fatigue  crack  growth  of  surface  flaws  loaded  in  cyclic  bending  in  a  transparent 
polymethylmethacrylate  (PMMA)  material. 
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Approach 


An  experimental  investigation  was  conducted  to  determine  the  influence  of  fatigue 
crack  closure  on  the  growth  of  surface  flaws.  The  investigation  was  undertaken  to  broaden 
the  understanding  of  surface  flaw  crack  opening  displacement  (COD),  closure  loads  and 
growth  rate  characteristics.  The  approach  is  made  possible  by  combining  a  unique 
interferometric  measurement  capability  with  a  transparent  polymethylmethacrylate  (PMMA) 
polymer  which  exhibits  fatigue  crack  growth  properties  which  are  represented  with 
conventional  fracture  mechanics  analysis  procedures.  The  transparent  polymer  material 
provided  the  ability  to  evaluate  three-dimensional  crack  opening  displacement  (COD) 
characteristics  of  surface  flaws.  Optical  interferometry  measurements,  using  a  helium-neon 
laser  light  source  for  the  first  time,  were  used  to  define  the  crack  opening  displacement 
profiles  and  crack  opening  loads.  The  fatigue  tests  were  conducted  on  rectangular  plate 
specimens  (Ref.  Figure  2.1)  loaded  in  four-point  bending. 

To  evaluate  plasticity  effects  on  closure,  four  different  load  histories  were  tested:  (1) 
constant  stress  intensity  factor  along  the  crack  surface  (AK^j),  (2)  constant  stress  intensity 
factor  into  the  crack  depth  (AK^),  (3)  constant  load  and  (4)  constant  amplitude  blocks 
loading  where  subscripts  "a"  and  "d"  identify  crack  tip  boundary  locations  A  and  D, 
respectively  (Figure  3.1).  COD  and  closure  measurements  were  made  for  various  crack 
lengths  and  aspect  ratios. 

The  investigation  was  organized  into  six  tasks  (Figure  3.2):  data  generation  (tasks  1- 
4),  data  reduction  (task  5)  and  data  analysis  (task  6).  A  description  of  these  tasks  is 
presented  below.  The  approach  for  tasks  1-4  is  presented  in  Section  IV,  experimental  results 
are  presented  in  Section  V,  and  data  analysis  in  Sections  VI  through  VIII. 
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o  Accuracy  of  an  dastic  TOD  prediction 


Figures^  Block  diagrm  of  the  experimental  test  and  data  analysis 


Experimental  Data  Generation 


The  objective  of  tasks  1-4  is  to  generate  closure  load  (Kd),  crack  opening 
displacement  (COD)  and  fatigue  crack  growth  rate  (da/dN  and  dc/dN)  data  for  surface  flaws 
in  a  polymethylmethacrylate  (PMM A)  specimen  subjected  to  cyclic  four-point  bending 
loads.  The  specimen  were  loaded  with  a  four-point  bending  fixture  mounted  on  an  MTS 
closed  loop,  electro-hydraulic,  fatigue  machine.  Four  loading  histories  were  generated  to 
isolate  plasticity  induced  closure  load  effects  —  two  low  level  constant  stress  intensity  factor 
(AK)  tests,  a  variable  AK  test  (constant  load),  and  a  block  loading  test.  An  x-y  traveling 
microscope  was  used  to  track  the  crack  growth  in  the  "a"  and  "c"  directions.  Photographs  of 
the  crack  interference  fringe  patterns  were  taken  as  load  was  applied  from  zero  to  a  level 
greater  than  the  closure  loads.  Fringe  pattern  photos  were  taken  at  "c"  crack  length  intervals 
of  approximately  0.05  in  (1.23  mm)  for  crack  sizes  which  ranged  between  0.02  and  0.60 
inch.  Additional  measurements  of  "a"  and  "c"  crack  lengths  and  closure  loads  were  taken 
when  unexpected  occurrences  were  noted  during  the  experiment. 


Experimental  Data  Reduction 

The  objective  of  task  5  is  to  organize  the  experimental  data  from  tasks  1-4  into  a 
cohesive  set  of  crack  opening  displacement  (COD)  and  closure  load  definitions.  Because  of 
the  diversity  of  information  required  to  address  the  closure  question,  experimental  data 
reduction  was  evaluated  from  several  perspectives:  the  effect  of  crack  opening  displacement 
(COD)  formation  on  closure  loads,  the  ability  to  quantify  COD  and  closure  loads,  the  effect 
of  the  four  load  histories  on  closure,  and  the  influence  of  COD  and  closure  in  defining 
AKgf^  These  evaluations  are  covered  in  Chapters  V  and  VI. 
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Experimental  Data  Analysis 


The  object  of  task  6  is  to  use  the  experimental  data  to  evaluate  fatigue  crack  growth 
(FCG)  and  COD  prediction  procedures.  This  is  accomplished  by  the  following  data 
reduction  and  analysis  steps: 

A.  Categorize  and  evaluate  crack  opening  displacement  (COD)  and  closure 
loads  as  a  function  of  crack  size. 

B.  Evaluate  the  influence  of  COD  and  plasticity  induced  wakes  on  closure 
loads.  Measure  closure  loads  from  constant  stress  intensity  factor  tests 
(tasks  1  &  2)  and  compared  to  those  from  a  constant  load  test  (tasks  3). 
Determine  if  closure  load  differences  are  attributable  to  load  history 
and/or  plasticity  induced  wake  patterns. 

C.  Compare  the  influence  of  plane  stress  and  plane  strain  on  COD  and 
closure  using  results  from  tasks  1  &  2.  Since  the  constant  stress  intensity 
factor  tests  is  conducted  by  maintaining  the  same  K  value  at  the  crack  tip 
locations  A  and  D  (<j»  =  0  and  90  degrees  in  Figure  2.1,  respectively), 
relative  differences  in  measured  closure  and  COD  are  emphasized. 

D.  Compare  fatigue  crack  growth  rate  analysis  using  applied  AK  and  AK^ff. 
Compare  predicted  and  experimental  FCG  data  (da/dN  &  dc/dN).  Use  the 
Newman-Raju  (5)  stress  intensity  factor  solutions. 

E.  Evaluate  the  ability  of  AK^ff  to  normalize  FCG  rate  predictions  using 
measured  and  an  assumed  closure  loads.  Analysis  of  all  tasks  provided 
insight  into  how  closure  loads  influence  fatigue  crack  growth.  Utilize 
Paris  Law  material  constants  from  available  literature  sources  when 
applicable. 


F.  Compare  predictions  of  elastic  crack  opening  displacement  (COD) 
amplitudes  to  measured  COD  amplitudes. 
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IV.  EXPERIMENTAL  PROCEDURES 


This  section  covers  the  experimental  procedures  used  in  this  investigation.  Included 
are  discussions  of  the  test  equipment,  loading  histories,  material  selection,  and  preparation  of 
the  specimens  for  fatigue  testing. 


Test  Procedures  and  Equipment 

Development  of  the  Newton  interferometry  system  for  mapping  crack  surface 
displacement  profiles  was  a  major  part  of  this  investigation,  and  is  an  original  experimental 
setup.  This  technique  provided  a  precision  method  for  measuring  closure  loads  and  crack 
opening  displacement  (COD)  patterns.  This  detailed  information,  in  combination  with  a 
carefully  selected  test  matrix,  provided  a  unique  data  set  for  understanding  the  influences  of 
the  closure  mechanism  on  fatigue  crack  growth  of  surface  flaws.  A  photograph  of  the 
experimental  arrangement  (MTS  electro-hydrolic  controller,  load  frame  with  the  specimen 
mounted  in  a  four-point  bending  fixture  and  the  laser  interferometry  system)  is  presented  in 
Figure  4.1. 


Newton  Interferometry  System 

A  laser  interferometer  system  was  used  to  record  optical  interference  fringe  patterns 
under  various  applied  loads  during  fiitigue  cycling  of  the  specimen.  The  Newton 
interferometer  system  combines  a  lO-milliwatt  helium-neon  laser  monochromatic  light 
source  (X  =  6328  A),  a  reflecting  mirror,  beam  splitter  and  35mm  camera  mounted  to  a 
three-way  adjustable  base  plate.  The  laser  monochromatic  light  source  was  used  instead  of  a 


Figure  4. 1  Experimental  arrangement  fo  the  MTS  electro-hydraulic 
controller,  load  frame  with  the  specimen  in  a  four-point 
bending  fixture,  and  the  laser  powered  interferometry 
system. 
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sodium  arc  light  source  (X  =  5890  A),  as  proposed  by  Packman  (16),  to  increase  the  light 
intensity  reflecting  from  the  crack  surface.  A  schematic  and  photograph  of  the  laser 
interferometry  system  are  presented  in  Figures  4.2a  and  b,  respectively.  Light  from  the  laser 
is  reflected  90  degrees  at  the  mirror,  split  into  two  beams  at  the  partial-mirror  beam-splitter 
where  one  beam  passes  into  the  specimen  and  reflects  off  the  two  crack  surfaces.  The  two 
reflected  rays  from  the  crack  surfaces  destructively  interfere,  due  to  a  180*^  waveshift  off  the 
crack  back  surface,  forming  fringe  patterns.  The  rays  then  pass  back  through  the  beam 
splitter,  thru  a  bellows  and  into  the  35-mm  camera.  The  beam  splitter  is  used  to  minimize 
the  angle  between  the  light  source  and  the  light  reflecting  back  to  the  camera,  providing  a 
direct  line  of  sight  between  the  crack  surface  and  the  camera. 

Since  the  fringe  patterns  represent  contour  lines  of  constant  crack  surface 
displacement,  the  interferometry  system,  with  the  transparent  PMM  A  material,  provides  the 
unique  capability  to  measure  very  small  displacements  over  the  complete  crack  interior 
surfece  (16,20,25).  By  counting  the  interference  fringe  orders  as  the  crack  opens  under  an 
applied  load,  displacement  changes  of  6.22E-6  in.  (1.58E-4  mm)  can  be  resolved  between 
zero  displacement  and  the  first  fringe  order  formation.  As  an  example,  fringe  order  #2  in  the 
schematic  diagram  of  Figure  4.3  represents  a  constant  displacement  between  the  two  crack 
surfaces  of  1.87E-5  in  (4.75E-4  mm)  in  accordance  with  the  optical  relationship  for  total 
crack  opening  displacement  given  by  eq.  4.1  (17): 

d  =  (2Ti-l)X/8  (4.1) 

Here  X  is  the  wave  length  of  light,  and  d  is  half  of  the  total  crack  surface  separation  for 
fringe  order  of  t]  =  1,2,3,...  The  zero  order  fringe  defines  the  location  where  the  crack 
surfaces  are  closed.  The  maximum  displacement  point  occurs  within  the  second  fringe  order 
contour  line  (in  this  example),  and  is  known  to  be  less  than  the  displacement  corresponding 
to  a  third  fringe  order  or  3.11E-5  in  (7.91E-4  mm).  This  topographical  information  can  be 
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igure  4.2  (Continui 


used  to  define  COD  profiles,  as  well  as  identify  profile  trends  with  applied  load  and  changes 
with  crack  size. 

Example  fringe  patterns  produced  by  the  laser  interferometry  system  are  shown  in 
Figures  4.4a  through  4.4d  using  a  test  specimen  with  a  0.187  inch  long  surface  crack  (2c). 
The  photos  show  relative  crack  opening  displacement  as  load  is  increased  from  zero  (4.4a)  to 
half  of  the  fully-open  load  (4.4b),  3/4  of  the  fiilly-open  load  (4.4c),  and  at  a  load  level  when 
the  crack  is  fully-open  (4.4d).  As  can  be  seen  in  Figure  4.4a,  the  crack  surfaces  and  crack  tip 
(location  A)  ate  open  (fringe  order  #1)  with  no  applied  load  to  the  specimen.  Note  the  four 
crack  boundary  locations  A,  B,  C  and  D,  the  crack  firee-surface  length  (2c)  and  depth  (a).  As 
a  bending  moment  is  applied,  the  crack  surfaces  open  toward  the  crack  free-surface  (Figures 
4.4b  and  4.4c),  and  are  completely  open  in  Figure  4.4d.  The  crack  ffee-surface  and 
reflection  from  the  PMMA  specimen  are  shown  in  Figure  4.4d.  The  crack  maximum 
measured  height  at  the  fully-open  load  level  is  shown  in  this  figure  by  the  13th  fringe  order 
to  have  a  total  displacement  between  crack  surfaces  of  1.6E-4  in.  (3.19E-4  mm).  The  plane 
stress/strain  constraint  effect  is  shown  by  location  C  being  the  last  portion  of  the  crack 
surfaces  to  fully-open  (Ref.  Figure  4.3). 

The  sequence  of  fringe  order  photographs  shown  in  Figure  4.4  are  an  example  set 
obtained  at  four  different  loadings.  During  the  experimental  tests,  interferometry  and  crack 
size  photos  were  taken  at  approximately  0.050  in  (1.27  mm)  increments  of  crack  size  along 
the  specimen  free-surface  (2c).  At  each  of  these  crack  lengths,  10  to  15  photos  were  taken  of 
the  interference  fringe  patterns  as  load  was  applied  to  the  specimen.  Analysis  of  this  data  for 
all  load  histories  is  presented  in  Chapter  V. 

Visual  observations  of  optical  interference  fringe  orders  were  used  to  determine 
closure/opening  load  levels.  Opening  and  closure  loads  are  used  interchangeably  here  since 
less  than  5%  difference  in  load  level  was  measured  to  achieve  the  same  displacement  when 
the  load  was  increased  on  the  specimen  (opening)  or  reduced  (closure).  However, 
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measurement  values  presented  in  this  report  are  ail  crack  opening  loads  except  where 
comparisons  of  opening  and  closure  load  are  identified  in  the  text. 


Digitizing  Crack  Displacement  Profiles 

Accurate  digital  measurement  of  the  interference  fringe  orders  was  accomplished  by 
using  a  20  inch  x  20  inch  digitizing  plate,  a  photo  negative  enlarger  (Figure  4.5),  and  the 
3Smm  negatives  generated  from  the  laser  interferometer  system.  The  35mm  negatives  were 
inserted  into  the  enlarger  to  expand  the  fringe  order  photos  on  the  digitizing  plate.  The 
digitizing  plate  was  connected  to  a  Techtronics  terminal  with  a  self-contained  magnetic  tape 
for  storing  the  digitized  data.  A  software  program  written  for  the  Techtronics  computer  was 
used  to  store  and  analyze  the  fringe  order  numbers  and  x-y  digital  displacement  data.  This 
information  was  then  used  to  construct  a  three-dimensional  contour  map  of  the  COD 
patterns. 

A  sample  set  of  digitized  COD  profrles  for  zero  and  various  applied  loads  is  shown  in 
Figure  4.6.  The  discrete  data  points  are  connected  by  straight  lines  to  shown  the  profile  at  an 
applied  load.  Digitizing  the  fringe  orders  followed  a  procedure  of  initially  locating  the  crack 
maximum  depth  "a"  to  define  location  A  (Figure  4.3).  Fringe  orders  were  then  located 
relative  to  the  crack  tip  (location  A)  and  digitized  along  a  line  between  locations  A  and  B. 
Figure  4.6  presents  half  of  the  full  crack  opening  displacement  (COD)  profile  for  various 
loads  -  with  displacement  plotted  along  the  ordinate,  and  distance  from  the  crack  tip 
location  A  on  the  abscissa.  Loads  from  zero  to  levels  above  the  non  linear  closure  load 
region  are  presented.  A  detailed  discussion  of  these  patterns  is  covered  in  Chapter  V. 
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X-Y  Traveling  Microscope 


An  x-y  directional  traveling  microscope  was  used  to  measure  the  "a"  and  "c"  crack 
dimensions  (Figure  4.3)  during  fatigue  testing.  The  microscope  had  a  10  X  magnification 
factor  with  a  4.5  in  (114.3  mm)  focal  length  for  viewing  the  crack  through  the  specimen. 

The  microscope  was  mounted  to  the  load  frame  and  measurment  accuracy  calibrated  with 
known  reference  marks  on  the  specimen.  Frequent  recalibration  checks  were  conducted  to 
maintain  accurate  and  consistent  crack  length  measurements.  The  accuracy  for  measuring 
the  2c  surface  crack  length  (Figure  4.3)  was  estimated  to  be  0.001  inch,  and  into  the 
specimen  depth  to  be  0.002  inch.  These  differences  are  attributed  to  the  traveling 
microscope  scale  markings,  and  linear  travel  distance  per  revolution  of  the  traveling 
mechanism.  It  was  assumed  that  the  readings  could  be  missed  by  ±  2  marks. 

The  microscope  was  calibrated  by  comparing  actual  measurements  to  known 
reference  marks  on  the  specimen.  The  reference  marks  were  placed  in  the  crack  plane  and 
the  microscope  aligned  on  the  fatigue  machine  to  remove  parallax  error.  Measurements  were 
then  validated  by  expanding  photographic  negatives  of  the  known  reference  marks  and  crack. 
It  is  estimated  that  over  1700  individual  crack  measurements  ("a",  "c"  and  closure  load)  were 
made  using  the  x-y  traveling  microscope.  These  measured  values  of  crack  lengths  and 
closure  load  (along  with  applied  load  levels  and  calculated  AK, 
closure  stress  intensity  &ctor  ratios)  are  included  in  the  appendix. 


Fatigue  Machine 

The  experiments  were  performed  on  a  closed  loop  MTS  fatigue  machine  with  a  5620 
pounds.  (25  kN)  load  cell  and  four-point  bending  fixture  (Figure  4.1).  The  MTS  load  cell 
was  calibrated  for  four  load  setting  ranges.  Range  1  was  100  %  of  the  5,620  pounds  (25  kN) 
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capacity,  range  2  was  50%,  range  3  was  20%,  and  range  4  was  10%  of  maximum  capacity. 
Load  range  3  was  used  during  precracking,  while  load  range  2  was  used  for  most  of  the 
fatigue  testing.  Load  range  4  was  used  exclusively  for  measuring  closure  load  amplitudes. 
This  setting  represented  a  measurement  sensitivity  of  562  pounds  (2.5  kN)  at  10  volts,  or 
56.2  pounds  (0.25  kN)  per  volt.  The  fatigue  fixture  exposed  the  specimen  to  a  bending 
moment  (M,  in-lb)  of  0.8125  times  the  applied  load  P  (lb).  That  is,  M  =  (3.25/4.0)  P  = 
0.8125  P  in-lb.  (M  =  0.0206  P  Newton-m,  for  P  in  Newtons). 


Load  History  Selection 

Four  experimental  load  histories  were  selected  to  evaluate  the  effect  of  plastic  zone 
patterns  on  plasticity  induced  closure.  The  four  plastic  zone  patterns  were  generated  by 
fatiguing  the  specimen  to  four  different  load  histories:  (1)  a  constant  stress  intensity  factor 
along  the  crack  surface  AK^  =  constant  at  location  D,  Figure  4.3);  (2)  constant  stress 
intensity  factor  into  the  specimen  depth  (AK^  =  constant  at  location  A);  (3)  a  constant  cyclic 
load  applied  to  the  specimen  and,  (4)  a  variable  block  loading.  The  stress  ratio  (R)  was 
maintained  throughout  the  test  at  0.035  to  be  as  close  to  zero  as  possible  while  maintaining  a 
positive  load  on  the  bending  specimen.  These  four  load  histories  were  briefly  described  in 
Section  III  as  task  1-4. 


Plasticity  Considerations 

The  four  load  histories  were  selected  to  emphasize  plastic  zone  size  differences  at 
locations  A  and  D  (predominent  plane  strain  and  stress  regions,  respectively)  on  the  crack 
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boundary  (Figure  4.3)  by  controlling  AK^  and  AK^j.  These  considerations  can  be  seen  by 
reviewing  the  stress  intensity  factor  (K)  equation  for  a  surface  flaw  in  bending  (5).  That  is: 


K  =  a  V(jia/Q)  Hf(a,c,  (|>) 

=  6M/bt2  V(jra/Q)  Hf(a,c,<t>)  (4.2) 

where  M  is  the  bending  moment,  b  and  t  are  the  specimen  width  and  thickness,  respectively, 
"a"  is  the  crack  depth  into  the  specimen  thickness,  "c"  is  the  half  surface  crack  length,  Q  is 
the  elliptical  correction  factor,  H  is  the  correction  term  for  bending,  and  f(a,c,(j))  is  the 
functional  Boundary  Correction  Factor  —  where  "a"  is  a  function  of  a/c  at  a  given  angle  (j)  on 
the  crack  boundary  (Figure  4.3).  Removing  the  constants: 

K  a  M  V(a/Q)  f(a,c,<j»)  (4.3) 

showing  that  K  can  be  defined  by  controlling  M  and  measuring  the  "a"  and  "c"  crack  lengths 
"  note  that  Q  is  a  function  of  a/c.  Thus,  by  either  controlling  or  measuring  these  three 
variables  (a,c  and  M),  an  experimental  approach  can  be  established  to  investigate  how 
closure  load  and  crack  opening  displacements  vary  with  K. 

The  influence  of  plane  stress  and  strain  (<j>  =  0  and  90  degrees,  respectively,  Figure 
4.3)  constraints  on  plastic  zone  size  was  also  evaluated  as  part  of  the  test  matrix 
requirements.  This  is  typically  accounted  for  by  a  plastic  constraint  factor  (pcf)  which  is 
defined  as  (27): 


pcf  -  /  Oyg  (4.4) 

where  cfujax  maximum  applied  stress  and  Oy^  is  the  material  yield  stress.  By  utilizing 
the  Von  Mises  yield  criteria  and  an  assumed  Poissons  r?tio  of  1/3  for  metal,  a  pcf  value  of 
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3.0  was  calculated  for  plane  strain  to  account  for  internal  yielding  constraints,  and  1  for 
plane  stress.  Since  the  internal  crack  surfoces  are  not  in  a  purely  plane  strain  condition,  as 
calculated  by  this  approach,  Irwin  (28)  used  an  average  pcf  value  of  1.68  instead  of  3  as  an 
effective  yield  stress,  such  that  the  plastic  zone  radius  (rp't  is  defined  as: 

rp  =  (K/1.68ay5)^/2ji  (4.5) 

=  (K/ays)2/6n 

Experimentally,  effective  values  of  pcf  (29)  have  been  found  to  vary  between  1.5  and  2.0  for 
metals.  Using  the  Irwin  relationships  for  plastic  zone  size  the  following  definitions  can  be 
established  for  plane  stress  (equations  2.4  and  2.5)  where  rQ  is  defined  at  (j>  =  0  degrees  as 
(Figure  4.3): 

ro  =(Ko/ays)2/2n  (4.6) 

and  plane  strain  where  rgQ  is  defined  at  4*  =  90  degrees,  as: 

rqo  =  (K90/  CTys)^  /  6jt  (4.7) 

If  AK  is  then  maintained  at  a  constant  level  at  two  locations  (<J>  =  0  and  90  degrees)  during 
separate  tests,  such  that  Kq  for  one  test  is  replaced  with  K9Q  for  another  test,  it  follows  that: 

"■o  =  ^  '■90  (^-S) 
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Therefore,  by  maintaining  a  constant  plastic  zone  size  at  (^  =  0  and  90  degrees,  such  that  rQ  = 
tqq,  it  follows  that: 


K90  =  V3  Ko  (4.9) 

These  differences  between  plane  stress  and  plane  strain  plastic  zone  size  are  utilized 
in  the  constant  AKd  and  AKa  tests  (tasksl  and  2)  to  emphasize  the  relative  influence  of 
plastic  zone  size  and/or  K  level  on  closure  loads.  Plastic  zone  size  was  not  measured 
experimentally,  but  was  calculated  using  equations  2.4  and  2.5.  The  results  of  these  tests  are 
discussed  in  Chapter  V,  in  the  "Closure/Opening  Loads”  section. 


Isolating  Plasticity  Effects 

Four  load  histories  were  selected  to  isolate  specific  plastic  zone  influences  on  closure 
loads  at  crack  tip  boundary  locations  A  and  D.  Although  plasticity  effects  cannot  be  isolated 
for  all  areas  around  the  crack  tip  boundary,  relative  effects  of  dominent  influences  can  be 
evaluated  at  specific  locations  using  distinctive  test  cycles.  The  four  cyclic  load  histories 
selected  to  experimentally  emphasize  relative  effects  of  plastic  zone  differences  on  closure 
and  fatigue  crack  growth  are  defined  in  tasks  1-4. 

Task  1:  A  constant  AK^  load  history  test  was  conducted  to  produce  a  constant  plastic 
zone  size  along  the  crack  tip  ffee-surface  (location  D,  (J)  =  0  degrees.  Figure  4.3)  by  load 
shedding  as  the  crack  size  increased.  The  sinusoidal  bending  load-shed  history  produced  a 
constant  AK^j  along  the  crack  face  shov’n  in  Figure  4.7.  A  schematic  of  the  plastic  zone  size 
is  theoretically  represented  by  Figure  4.8.  A  constant  stress  intensity  factor  (AK^j)  of  600 
psi*in^^  (659  KPa*M^^^)  was  maintained  in  this  test. 
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Task  2:  A  constant  AK^  load  history  test  was  conducted  to  produce  a  constant  plastic 
zone  size  into  the  specimen  thickness  (location  A,  (^  =  90  degrees,  Figure  4.3)  by  load 
shedding  as  the  crack  size  increased  during  cyclic  loading.  The  sinusoidal  bending  load-shed 
history  produced  the  AK^  levels  sh  own  in  Figure  4.9.  Representation  of  the  theoretical 
plastic  zone  is  depicted  by  Figure  4.8.  A  constant  AK^  of  600  psi*in^^  (659  KPa*M^^^) 
was  maintained  into  the  crack  depth  in  order  to  make  direct  comparisons  with  data  generated 
in  task  1  (as  discussed  in  the  previous  section,  "Plasticity  Considerations"). 

Task  3:  A  constant  amplitude  sinusoidal  bending  load  test  was  conducted  to  produce 
a  continually  increasing  plastic  zone  size  along  the  crack  tip  boundary.  The  load  history 
produced  the  AK^  levels  shown  in  Figure  4.10.  The  continually  increasing  plastic  zone  size 
resulting  from  crack  growth  under  a  constant  load,  which  produces  a  continually  increasing 
AK,  is  schematically  depicted  in  Figure  4.11. 

Task  4:  A  block  loading  test  (Figure  4.12)  was  conducted  using  constant  moment 
segments  to  evaluate  the  effects  of  large  plastic  zones  on  closure  and  fatigue  crack  growth. 
Each  block  loading  segment  was  maintained  at  a  constant  bending  moment  until 
predetermined  AK^}  levels  were  reached.  The  initial  block  loading  segment  used  a  AM  =  644 
in-lb,  or  a  AK^j  level  which  increased  from  608  to  657  psi*in^/2  (40%  to  44%  of  a  published 
reference  fracture  toughness  value  of  1500  psi*inl/2^  Ref.  27).  The  load  was  then  increased 
20%  to  a  ABC^  of  790  (53%  of  and  the  specimen  fatigued  until  a  AKj  value  of  996 
(66%  of  K|jj)  was  reached.  The  load  was  then  decreased  20%  to  the  original  starting  load 
level  of  AM  =  644  in-lb,  or  AK^j  =  830  (55%  of  K^^^).  Testing  continued  at  this  load  level 
until  a  AKjj  of  1191  (80%  of  Kj^^)  was  achieved.  The  final  block  loading  segment  was 
accomplished  by  reducing  the  load  level  50%  to  the  initial  baseline  AK^j  leve>  ^^08 
psi*in^^),  and  testing  continued  until  a  large  crack  size  (a  »  0.5,  or  approx  half 

through  the  specimen  thickness)  was  reached  (AK^j  =  836).  The  intent  of  this  cyclic  block 
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CONSTANT  LOAD  TEST 
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BLOCK  LOADING  TEST 


Kguxe  4.12  Actual  block  loadins  pattern,  and  resulting  stress  intensity 
facttn  calculated  at  the  crack  tre  free-suttree  (AKd).  to 
generate  large  plastic>zone  wakes 
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loading  pattern  was  to  grow  the  crack  through  two  large  plastic  zones.  The  crack  was  grown 
at  least  20  plastic  zone  diameters  past  the  initial  crack  tip  location  before  changing  the  load 
level. 


All  four  tests  used  the  same  precracking  procedures  to  establish  a  consistent  baseline 
starting  AK  of  600  psi*in^^.  The  two  constant  stress  intensity  factor  tests  were  maintained 
at  locations  A  and  D  (tasks  1  and  2)  to  within  ±3.7%  (±22  psi'in^'^  or  ±24  KPa*M^^)  by 
regularly  measuring  the  "a"  and  "c"  crack  dimensions  and  calculating  K  levels  at  those  crack 
sizes.  Once  a  level  less  than  620  psi*in^^^  (680  KPa*M^^)  was  reached,  the  load  was 
reduced  to  achieve  a  level  which  was  greater  than  580  psi*in^^  (637  KPa*M^^^)  and  the 
fatigue  cycling  continued. 

Tasks  1  and  2  were  selected  to  emphasis  the  relative  differences  between  the  larger 
plane  stress  and  smaller  plane  strain  plastic  zone  at  locations  A  and  D  on  the  crack  boundary. 
That  is,  controlling  AK^j  along  the  crack  free-surface  (task  1)  to  emphasize  the  plane  stress 
(0^=0)  influences.  Likewise,  controlling  AK^  into  the  crack  depth  (task  2)  to  emphasize  the 
plane  strain  (t^-0)  influences  on  closure  relative  to  task  1. 

Comparison  of  closure  load  differences  due  to  a  constant  plastic  zone  wake  pattern 
and  a  continually  increasing  wake  pattern  (Figures  4.8  and  4.11,  respectively)  is 
accomplished  by  comparing  results  from  task  3  with  tasks  1  and  2. 

The  block  loading  test  (task  4)  was  conducted  to  compare  relative  effects  of  large 
plastic  zone  wakes  on  closure  loads.  That  is,  will  the  closure  load  levels  for  cracks  which 
grow  through  large  plastic  zones  be  significantly  different  from  those  generated  by  the 
constant  AK  tests? 

The  four  load  histories  provided  a  unique  set  of  data  for  evaluating  the  influence  of 
plastic  zone  size  on  the  mechanisms  of  closure.  These  test  results  (task  1-4)  are  then  used  to 
make  comparisons  between  the  individual  and  combined  influences  of  crack  opening 
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displacement  (COD)  and  closure  loads  on  fatigue  crack  growth  characteristics,  as  discussed 
in  Chapters  V  and  VII. 


Specimen  Material  Selection 

To  meet  the  research  objectives  of  measuring  surface  flaw  internal  crack  opening 
displacements  (COD)  and  closure  loads  simultaneously,  it  was  desirable  to  use  a  transparent 
material  whose  crack  growth  properties  are  described  using  conventional  fracture  mechanics. 
This  was  accomplished  by  selecting  polymethylmethacrylate  (PMM A),  a  polymer  which  has 
been  used  extensively  over  the  past  IS  years  for  fracture  mechanics  investigations 
(16,18,20,25,30,31).  Although  PMMA  has  been  found  to  be  a  very  good  polymer  for 
conducting  fatigue  crack  growth  investigations,  it  is  also  important  to  recognize  that  PMMA 
is  viscoelastic,  and  that  its  crack  growth  characteristics  can  change  under  certain  test 
conditions.  For  this  reason,  a  great  deal  of  attention  was  paid  to  controlling  and  evaluating 
the  influence  of  variables  which  can  change  the  polymer  crack  growth  characteristics  during 
the  fatigue  test  (such  as  test  frequency,  environmental  temperature,  viscoelastic  relaxation, 
etc.). 

One  of  the  more  comprehensive  reviews  of  polymers  was  conducted  by  Hertzberg 
and  Manson  (32).  Their  reviews  identified  the  major  test  variables  (test  frequency, 
temperature,  etc.)  which  need  to  be  considered  when  conducting  fatigue  crack  growth 
investigation  with  polymer  materials.  These  variables  have  undergone  a  thorough  review  for 
this  research,  and  essential  variables  controlled  during  the  test.  The  following  sections 
discuss  these  variables,  and  address  both  the  positive  and  negative  aspects  of  using  PMMA 
material  for  this  investigation. 
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Material  Selection 


Polymethylmethacrylate  (PMMA)  is  a  transparent  polymer  which  has  been  used  in 
numerous  fracture  mechanics  investigations  (16,18,20,25,31,32).  A  review  of  polymers  by 
Hertzberg,  Manson  and  Wu  (33)  found  that  PMMA,  polycarbonate  (PC),  acrylonitrite- 
buladience-styrene  (ABS)  resin,  polyethylene  (LDPE)  and  nylon  66  are  a  few  of  the  more 
widely  used  polymers  which  follow  the  familiar  fracture  mechanics  power  function 
relationship  between  fatigue  crack  growth  rate  and  stress  intensity  factor  (Figure  4.13). 

The  positive  aspects  of  using  the  transparent  PMMA  as  a  model  material  for  crack 
closure  studies  are  summarized  by  their  capability  to: 


Exhibit  conventional  fatigue  crack  growth  properties  described  by  fracture 
mechanics;  power  law  relationship  of  da/dN  versus  AK,  Figure  4.14. 

Provide  the  opportunity  to  measure  three-dimensional  crack  closure,  crack 
opening  displacement  (COD)  and  crack  growth  patterns  (a  and  c) 
simultaneously  since  it  is  a  transparent  material  (not  directly  available  with 
metals). 


Exhibit  closure  load  ratios  (closure  load  to  maximum  applied  load)  which  are  in 
the  same  range  as  reported  for  ductile  and  brittle  metals,  as  well  as  analytical 
predictions  (Table  2.1). 


The  limitations  of  PMMA  as  a  model  material  for  studying  fatigue  crack  growth  are 
also  important  to  the  success  of  this  investigation.  The  primary  variables  which  need  to  be 
controlled  and/or  evaluated  during  experimental  fatigue  crack  growth  testing  are: 

-  Test  frequency 

-  Environmental  temperature  and  humidity 

-  Viscoelastic  relaxation  with  time 


Flgoxe  4.13 


Fatigue  Click  giowtfa  FUe  curves  for  sevenl  ensineering 
polymers;  corves  for  polymetfayimetfaacryitte  (PNA4A), 
low  density  polyetfayiene  <LDPE), 
acrvionitriie-tetadienoe-styvene  resin  (ABS),  polycarbonate 
(PC),  and  nylon  66  0^.  32) 
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Figuxe4.14  Fatigue  crack  growth  me  curve  for  PMMA  (ref.  19) 


Each  of  these  variables,  and  their  relationship  to  this  investigation,  are  discussed  **'.  the 
following  sections. 


Polymer  Characteristics 

As  with  metals,  PMM A  fatigue  crack  growth  rates  are  dependent  upon  frequency, 
loading  wave  form,  environment  and  temperature.  Of  the  six  categories  reviewed  by 
Hertzberg  and  Manson  (32)  (molecular  characteristics,  chemical  changes,  homogeneous  and 
inhomogeneous  changes,  transition  phenomena  and  thermal  heating),  the  effect  of  local 
thermal  heating  due  to  test  frequency  is  the  area  which  requires  the  greatest  control  for  the 
present  investigation.  The  other  variables  were  controlled  by  using  the  same  material  batch 
and  annealing  procedures  for  all  the  tests,  as  well  as  conducting  the  tests  at  room 
temperature. 

Experimenters  have  shown  that  a  test  frequency  above  approximately  10  Hz  can  have 
a  significant  effect  on  PMMA  crack  growth  rates.  Atermo  and  Ostberg  (34)  showed  that 
higher  test  frequencies  cause  excessive  heat  build-up  at  the  crack  tip.  It  is  hypothesized  that 
above  10  Hz,  the  conductivity  of  a  polymer  is  not  sufficient  to  dissipate  heat  away  from  the 
crack  tip.  This  results  in  a  local  crack  tip  heat  build-up  which  causes  the  local  crack  tip 
region  to  reach  a  "glass  phase"  state.  In  this  "glass  phase,"  polymers  will  absorb  moisture 
and  exhibit  different  fatigue  crack  growth  characteristics.  It  was  observed  by  Hertzberg  and 
Manson  that: 


"...although  heat  transfer  from  the  plastic  zone  to  its  cooler  surrounding 
environment  might  limit  the  rate  of  crack-tip  heating,  fatigue  testing  at  high 
frequencies  should  nevertheless  produce  a  finite  temperature  rise  at  the  crack 
tip.  To  wit,  Atermo  and  Ostberg  [79]  recorded  a  maximum  increase  in  crack- 
tip  temperature  of  up  to  20  degrees  K  in  fatigue  testing  of  PVC,  PMMA,  and 
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PC  at  11  Hz.  With  a  significant  increase  in  temperature,  yielding  processes  in 
the  material  surrounding  the  crack  tip  should  be  enhanced  and  lead  to  an 
increase  in  the  crack-tip  radius.  This  greater  radius  of  curvature  at  the  crack  tip 
should  result  in  a  lower  effective  AK.  As  the  effective  AK  decreases,  the 
fotigue  crack  growth  rate  is  expected  to  decrease  accordingly." 


The  potential  for  absorption  of  moisture  at  the  crack  tip  was  also  controlled  during 
the  fatigue  test.  Although  effects  of  humidity  are  only  known  to  occur  when  the  polymer 
reaches  its  "glass  phase"  temperature  (a  condition  prevented  by  testing  at  a  low  frequency), 
precautionary  measures  were  taken  by  sealing  the  specimen  surface  around  the  crack  free- 
surface  with  a  desiccate  material. 


Relaxation  of  PMMA 

Thru-thickness  Cracks.  Pitoniak,  et  al.  (20)  reported  that  PMMA  crack  opening 
displacement  (CX)D)  profiles  tend  to  relax  (change  in  crack  internal  surface  displacement) 
under  zero  load.  He  observed  that  by  allowing  the  specimen  to  relax  for  1, 15  and  180 
minutes  after  fatiguing  the  specimen,  the  internal  COD  magnitude  would  decrease.  Utilizing 
information  provided  in  the  report,  this  displacement  was  calculated  to  be  approximately 
2.49E-5  in  (6.3E-4  mm),  or  an  18%  change  in  the  total  displacement  under  zero  load 
between  1  minute,  and  180  minutes  Half  of  this  relaxation  was  noted  to  occurr  within  the 
first  IS  minutes.  These  observations  initiated  a  similiar  investigation  in  this  experiment  of 
surface  flaws.  Results  of  this  efrbrt  are  discussed  in  the  next  section. 

Surface  Cracks.  An  extensive  evaluation  of  surface  crack  relaxation  in  PMMA  was 
conducted  during  this  investigation.  The  need  for  the  evaluation  was  based  on  prior  findings 
for  thru-thickness  compact  tension  tests  discussed  above.  The  surface  crack  evaluation  was 
accomplished  by  measuring  crack  relaxation  at  various  crack  sizes  as  the  specimen  was 
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fatigued.  The  review  included  the  recording  of  crack  displacement  changes,  as  measured 
using  the  laser  interferometry  system,  after  various  relaxation  hold  times  at  zero  load  (no 
external  load  on  the  specimen)  and  with  loads  above  the  closure  load  levels. 

Example  tracings  of  enlarged  35mm  photographic  negatives  showing  COD  profile 
patterns  before  and  after  the  relaxation  hold  times  are  presented  in  Figures  4.15a  and  b  for 
two  crack  sizes.  The  figures  are  enlarged  to  emphasize  the  differences,  but  are  proportional 
to  each  other.  The  solid  fringe  order  contour  lines  were  recorded  immediately  after  testing, 
and  the  dotted  lines  15  minutes  later.  Figure  4.15a  presents  the  smaller  crack  size  (a  =  0.12, 
c  =  0.14)  under  an  applied  load  of  17%  of  the  maximum  &tigue  load.  No  changes  were 
noted  under  zero  load  for  this  crack  size.  Figure  4.15b  presents  the  larger  crack  size  (a  = 
0.19,  c  =  0.30)  under  zero  load.  Note  the  larger  number  of  fringe  orders,  representing  a 
larger  internal  displacement,  but  only  minor  differences  due  to  the  relaxation  hold  time. 

An  evaluation  of  the  load  needed  to  return  a  crack  to  its  original  displacement  after  a 
20>minute-relaxatioQ  time  was  conducted  for  one  crack  size  (a  =  0.079  in  and  c  =  0.097  in). 
This  restoring  bending  moment  was  found  to  be  0.65  in-lb  (0.0734  N-m)  ~  that  is,  the 
moment  needed  to  compensate  for  the  relaxation  effect.  This  is  0.11%  of  the  maximum 
bending  monent  applied  to  the  specimen  during  fatigue  testing. 

Repeatability  of  the  loading  and  unloading  measurement  process  was  also  calibrated 
during  the  tests.  Upon  loading  a  specimen  with  a  crack  size  of  a  =  0.074  in  (1.88  mm)  and  c 
=  0.097  in  (2.46  mm),  the  first  fringe  was  formed  at  a  load  level  of  47  in-lb  (4.8  N-m),  and 
the  crack  free-surface  initially  opened  at  67  in-lb  (7.6  N-m).  Continuing  to  increase  the  load 
level  above  the  COD  elastic  region  (approximately  152  in-lb  or  17.2  N-m)  and  then 
unloading  the  specimen,  the  first  fringe  order  formation  occured  at  38  in-lb  (4.3  N-m). 
Continuing  to  unload  the  specimen  to  zero  and  then  reloading,  the  first  fringe  order  was 
formed  at  41  in-lb  (4.6  N-m)  ~  compared  to  42  in-lb  during  initial  loading.  This 
repeatability  was  demonstrated  approximately  20  minutes  after  the  initial  load  cycle.  These 
measurements  represent  a  repeatability  of  approximately  2%  for  an  "opening"  load 
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measurement  process,  and  approximately  10%  difference  between  measurements  when 
"loading"  versus  "unloading"  for  this  crack  size. 

In  general,  the  variability  of  opening  load  measurements  was  about  2%  for  smaller 
crack  sizes  —  where  the  opening  load  levels  were  greater  than  40  in-lb  (4.5  N-m).  For  larger 
crack  sizes  (where  the  moment  required  to  open  the  crack  was  25  in-lb  instead  of  40  in-lb) 
the  variability  increased  from  2%  to  approximate  10%. 


Specimen  Preparation 

A  great  deal  of  care  was  taken  to  develop  a  repeatable  crack  initiation  site  and 
precracking  load  cycle  procedure.  This  effort  was  undertaken  to  minimize  precracking 
effects  on  the  fatigue  crack  growth  experimental  results.  Primary  consideration  was  given  to 
repeatability  of  the  procedures  from  specimen  to  specimen.  Other  important  specimen 
preparation  procedures  discussed  in  this  section  include  annealing  and  end  polishing. 


Specimen  Size 

The  experimental  tests  used  a  0.75  x  3  x  7  inches  (19.1  x  76.2  x  177.8  mm)  specimen 
with  a  part-thru  crack  on  one  face  (Figure  2.1  and  4.16).  The  specimen  size  was  selected 
after  reviewing  recommendations  from  two  primary  sources.  Jolles,  McGowan  and  Smith 
(36)  indicated  that  a  plate  specimen  with  a  width  of  two  and  a  half  times  the  crack  length 
(2c)  experienced  the  same  stress  intensity  as  those  for  plates  of  infinite  width.  Anderson  (6) 
suggested  that  the  width  be  two  and  a  half  times  the  crack  length  and  the  specimen  length  be 
two  times  the  width.  For  this  experimental  investigation,  the  selection  of  a  maximum  crack 
length  of  1.20  in  (30.5  mm)  met  both  of  these  specimen  size  criteria. 
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specimen  center 


Specimen  Annealing 


The  specimen  annealing  process  followed  those  used  by  Perez  (37 j  lo  minimize 
potential  residual  stress  developed  in  the  specimen  during  machining.  The  PMM  A  specimen 
was  annealed  at  100  degrees  centigrade  for  24  hours  prior  to  precracking.  The  cool  down 
cycle  was  performed  at  a  rate  of  approximately  3.5  degrees  C  per  hour  until  the  oven  reached 
SO  degrees  C.  The  specimen  was  then  allowed  to  cool  down  to  room  temperature  with  the 
oven  door  closed.  Total  cooling  time  was  approximately  48  hours. 


Specimen  Polishing 

In  order  to  utilize  the  laser  interferometry  system  for  measuring  interference  patterns, 
it  was  necessary  to  have  a  flat  transparent  finish  on  one  end  of  the  PMMA  specimen.  The 
end  polishing  procedures  developed  at  the  Air  Force  Wright-Patterson  Materials  Directorate 
(38)  are  as  follows: 

1.  Select  the  end  with  the  smoothest  surface  before  beginning.  Grind  that  end  of  the 
sample  with  an  8  in  diameter  polisher  at  350  rpm  using  an  8-in  PSA  silicon  carbide  600-grit 
disc.  This  operation  should  continue  until  the  surface  is  relatively  smooth. 

2.  Wash  the  specimen  in  soap  and  warm  water  to  remove  all  residue;  flush  with 
methanol  and  blow  dry  with  hot  air. 

3.  Polish  the  specimen  on  an  8  inch  diameter  polisher  using  8  inch  PSA  nylon  cloth, 
3-micron  diamond  paste  and  oil  lubricant.  Polish  at  350  rpm  until  surface  deformations  from 
the  600-grit  polishing  operation  are  completely  removed. 
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4.  Clean  the  specimen  as  in  step  #2. 


5.  Polish  the  specimen  on  an  8  inch  diameter  polisher  with  an  8  inch  PSA  micro 
cloth  using  0.06-micron  colloidal  silicon  suspension  at  200  rpm  until  the  surface  is  smooth 
and  transparent. 

6.  Clean  the  specimen  as  in  step  #2. 


Crack  Initiation  Site 

The  development  of  an  accurate  and  repeatable  crack  initiation  site  was  an  important 
part  of  the  precracking  process.  It  was  initially  decided  that  a  material  removal  approach 
would  be  used  to  eliminate  buildup  of  excess  material  experienced  by  an  indent  procedure. 
Several  material  removal  procedures  were  considered:  (1)  feeding  the  specimen  into  a 
rotating  foil,  (2)  scribing  a  straight  line  on  the  specimen  with  a  knife  edge,  and  (3)  scribing  a 
mark  with  a  rotating  knife  edge.  The  rotating  knife  edge  was  selected  because  it  best 
represented  a  surface  crack  shape,  and  was  the  most  uniform  and  repeatable  crack  initiation 
method. 

The  crack  initiation  site  procedure  was  developed  to  control  the  cut  length  and  depth 
dimensions.  This  was  accomplished  by  using  a  lathe,  with  the  PMMA  specimen  mounted  on 
the  traversing  be^  and  a  knife-edge  cutting  tool  mounted  to  the  lathe  rotating  head.  After 
mounting  the  specimen  on  the  lathe,  the  bed  was  traversed  to  the  knife-edge  which  was 
brazed  to  a  circular  rod.  The  knife-edge  was  located  next  to  the  specimen  mid-surface  using 
a  thin  shim,  then  rotated  away  from  the  specimen  and  the  lathe  table  (with  the  specimen) 
traversed  toward  the  cutter  to  set  the  initiation  site  depth.  The  knife-edge  cutter  was  then 
rotated  into  the  specimen,  nicking  the  PMMA  material.  This  procedure  scribed  a  consistent 
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crack  initiation  site  of  0.005  inch  (0.13  mm)  deep  by  0.060  inch  (1.5  mm)  long  into  the 
specimen  surface  (Figure  4.17).  The  knife-edge  radius  of  curvature  was  0.0925  inch.  All 
specimens  were  scribed  at  the  same  time  to  ensure  consistency  of  the  crack  initiation  site 
(Figure  4.16).  Figure  4.17  also  identifies  an  "initial  crack"  site  (dotted  lines)  where  the  real 
crack  growth  began  versus  the  scribed  initiation  site  generated  by  the  knife-edge  cutting  tool 
(solid  line).  This  "initial  crack"  location  and  subsequent  growth  pattern  were  typical  for  all 
tests.  With  fatigue  cycling,  the  "initial  crack"  encompassed  the  initiation  site  during 
precracking,  and  grew  to  at  least  1.5  times  the  crack  initiation  site  length  along  the  free- 
surface,  and  8  times  the  initiation  site  depth  into  the  specimen  (as  discussed  in  the  next 
section). 


Precracking  Procedures 

The  precracking  load  shed  procedures  followed  ASTM  standard  practices  for  a  thru¬ 
thickness  crack,  plus  additional  criteria  established  for  this  investigation  of  surface  flaw. 

The  load  level  was  kept  below  80%  of  the  material  yield  strength,  and  typically  ran  below 
30%  of  yield  after  precracking.  The  load  shedding  cycle  was  established  to  meet  the  plastic 
zone  size  (or  stress  intensity  factor  levels)  requirements.  This  was  accomplished  by 
measuring  "a"  and  "c"  crack  lengths  with  the  traveling  microscope,  and  calculating  stress 
intensity  factors  (K)  using  the  Newman-Raju  (5)  solutions  to  determine  when  the  load  should 
be  reduced.  The  load  shedding  procedure  maintained  a  constant  stress  ratio  of  R  =  0.035. 
This  ratio  was  selected  to  be  close  to  zero  while  maintaining  a  positive  load  on  the  specimen 
throughout  the  test. 

A  concentrated  effort  was  initiated  to  define  a  precracking  load  cycle  which 


combined  both  the  ASTM  standards  (39)  and  additional  criteria  established  for  this  tests 
(items  2-4  below).  Since  there  were  no  ASTM  precracking  standards  for  surface  flaws,  the 
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Figure  4.17  Cross-section  schematic  of  the  symmetrical  crack  initiation 
site  showing  the  cutting  tool  radius  of  curvature 


standards  and  the  intent  of  the  ASTM  thru-crack  standards  were  followed.  This  was 
accomplish  by  writing  an  iterative  computer  code  to:  (1)  meet  ASTM  thru-crack  standards, 
(2)  initiate  all  experimental  tests  at  the  same  crack  tip  plastic  zone  size  along  the  crack 
surface,  (3)  ensure  the  crack  size  before  the  fatigue  crack  growth  analysis  began  after 
precracking)  was  sufficiently  large  to  encompass  the  crack  initiation  site,  and  to  be 
accurately  recorded  by  the  laser  interferometry  photographic  system  and  the  traveling 
microscope,  and  (4)  ensure  the  plastic  zone  size  was  decreased  after  each  load  shedding  step. 
In  effect,  the  iterative  computer  code  defined  a  precracking  load  shed  cycle  while 
maintaining  ASTM  and  other  criteria  for  this  investigation  of  surface  flaws. 

The  ASTM  load  shedding  criteria  states  that  the  load  shedding  level  shall  be  less  than 
20%  of  the  maximum  applied  load,  and  the  crack  should  grow  at  least  3  times  the  prior 
plastic  zone  diameter.  Additional  criteria  established  under  this  investigation  of  surface 
flaws  were  the  minimum  to  maximum  load  level  (R)  be  maintained  at  0.035,  the  maximum 
load  level  would  not  exceed  80%  of  the  material  yield  strength,  a  crack  length  "c"  would  be 
greater  than  0.010  inch  (0.25  mm)  before  the  first  load  shedding  occurred  to  produce  an 
accurate  crack  length  measurement,  the  load  levels  were  not  above  the  critical  stress  intensity 
factor  nor  below  the  threshold  stress  intensity  levels,  and  that  the  plastic  zone  size  was 
derveasing  from  the  prior  step  to  the  end  of  the  next  load  shed  step. 

To  meet  the  above  criteria,  an  assumed  value  of  crack  aspect  ratio  was  required  for 
the  iterative  computational  analysis.  From  experimental  tests  with  PMMA,  and  experiments 
by  Com  (40)  and  Mahmoud  and  Hosseini  (41)  using  various  metals,  the  a/c  ratio  was  found 
to  initially  increase  from  a  value  less  than  1.0  to  a  value  greater  than  1.0,  and  then  to 
decrease  below  1.0.  For  small  surface  cracks,  an  a/c  value  of  1.0  was  found  to  be  typical.  A 
sensitivity  study  of  a/c  effects  on  K,  using  the  Newman-Raju  stress  intensity  factor 
calculations,  established  that  aspect  ratios  in  the  range  of  0.82  <  a/c  <1.0  only  changed  the  K 
solutions  a  small  amount.  This  information,  plus  additional  calculations  discussed  in  the 
Appendix  computer  section,  led  to  the  selection  of  0.85  as  an  average  value  of  crack  aspect 
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ratio  for  the  precracking  load  analysis.  The  value  of  a/c  equal  to  the  initiation  site 
dimensions  (0.005/0.03  =  0.17)  was  not  realistic  since  the  "initial  crack"  always  started  on 
the  scribed  initiation  site  surface  (Figure  4.17)  and  propagated  in  a  semicircular  pattern  as 
described  in  the  last  section. 

The  crack  growth  analysis  in  this  report  only  uses  measurements  taken  after 
precracking;  although  all  precracking  and  crack  growth  experimental  data  are  presented  in 
Appendix  A  through  D.  Per  ASTM  standard  E647  (39),  the  crack  growth  data  begin  when 
the  crack  has  grown  more  than  three  times  the  plastic  zone  diameter  of  the  prior  load  leN-el. 
This  criteria  resulted  in  the  fatigue  crack  size  after  precracking  being  at  Irast  1.5  times  the 
crack  initiation  site  size  along  the  surface  length,  and  8  times  the  initiation  site  size  into  the 
specimen  depth. 


Minimizing  Precracking  Influences 

In  addition  to  utilizing  identical  specimen  annealing,  crack  initiation  site  and 
precracking  load  shedding  procedures  on  each  specimen,  the  final  data  set  was  reviewed  to 
ensure  the  crack  free-surface  length  had  grown  at  least  three  times  the  plastic  zone  diameter 
after  the  last  precracking  load  shedding  step.  This  procedure  was  followed  to  ensure  the 
experimental  crack  growth  data  used  in  the  analysis  followed  ASTM  precracking  standards. 

For  example,  the  final  precracking  load  shed  surface  stress  intensity  factor  (AK^)  was 
600  psi*in^^^  (659  KPa*M^^^).  For  the  constant  AK^j  data  set  (task  1)  after  precracking,  the 
crack  length  c  =  0.0365  inch  (0.93  mm),  and  the  plastic  zone  diameter  is  8.73E-3  inch  (0.07 
mm).  Then,  for  the  crack  growth  data  to  be  used  in  the  analysis,  the  crack  length  was 
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required  to  be  3  times  the  plastic  zone  diameter,  or: 


c  =  0.0365  +  3.  *  (0.0087)  (4. 1 0) 

=  0.0452  (4.11) 

Therefore,  all  data  which  had  a  "c"  crack  length  less  than  0.045  inch  (1.15  mm)  is  defined  as 
precracking  data,  and  not  used  in  the  fatigue  crack  growth  analysis.  This  same  approach  was 
followed  for  all  fatigue  crack  growth  data  sets. 
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V.  CRACK  OPENING  DISPLACEMENT  AND 


CLOSURE  LOAD  MEASUREMENTS 

The  geometric  crack  opening  displacement  (COD)  patterns  and  the  closure  load 
measurements  generated  from  the  experimental  tests  for  surface  flaws  in  bending  are 
presented  in  this  chapter.  Because  of  the  complex  nature  of  fatigue  crack  closure,  a  number 
of  new  definitions  were  established  to  describe  and  organize  the  experimental  data.  Three 
new  categories  of  COD  and  three  defrnitions  of  closure  load  were  identified.  To  differentiate 
between  these  effects,  the  chapter  is  divided  into  two  major  sections  --  the  first  describing  the 
geometric  COD  patterns  which  related  to  the  whole  crack,  and  the  second  describing  the 
closure  load  measurements  needed  to  define  effective  stress  intensity  factor  (^K^ff)*  The 
implications  of  these  results  for  defining  A^^ff  discussed  in  Chapter  VI,  and  their 
influence  on  predicting  fatigue  crack  growth  in  Chapter  VII.  Comparisons  between  predicted 
and  experimental  COD  data  are  presented  in  Chapter  VIh. 


Crack  Opening  Displacement  (COD)  Profrles 

The  experimental  investigation  of  surface  flaws  in  bending  has  identified  three 
primary  crack  opening  displacement  (COD)  patterns  defined  here  as  types  I,  II  and  III.  The 
uniqueness  of  these  COD  types  is  embodied  in  the  crack  surface  displacement  profile 
patterns  under  zero  load,  and  during  load  application.  Besides  the  pure  geometric  differences 
between  the  crack  types,  each  has  a  direct  impact  on  crack  closure/opening  load  level  and 
crack  tip  stress  intensity  factor,  and  thus,  the  definition  of  effective  stress  intensity  factor. 
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A  description  of  the  COD  patterns  as  they  apply  to  crack  size,  applied  load,  and 
fatigue  load  history  is  presented  in  this  section.  Crack  size  is  found  to  be  a  primary  factor  in 
organizing  COD  profile  pattern  data  as  will  be  described  in  the  following  sections. 


Typical  COD  Profiles 

The  experimental  investigation  identified  three  primary  crack  types  which  were 
observed  to  be  consistent  for  all  tests.  Each  crack  type  is  unique  because  of  its  influence  on 
crack  tip  stress  intensity  factor  (K)  at  various  crack  tip  boundary  locations.  The  crack  tip 
locations  identified  in  Figure  S.l  will  be  referred  to  extensively  throughout  this  report. 
Location  A  defines  the  maximum  crack  depth  location,  location  B  the  midpoint  of  the  crack 
free-surface,  location  C  is  approximately  12  degrees  from  the  crack  tip  free-surface,  and 
location  D  is  located  at  the  crack  tip  free-surface.  The  angle  f  identifies  points  along  the 
crack  perimeter,  and  is  measured  from  the  crack  free-surface  with  respect  to  line  AB  (not  the 
elliptical  angle).  The  crack  is  assumed  to  be  symmetrical  about  line  AB.  The  fringe  order 
lines  (constant  displacement  contours)  shown  in  the  Figure  are  described  in  Section  IV, 
"Newton  Interference  System." 

The  three  crack  types  are  schematically  shown  in  Figure  5.2  for  a  zero  applied  load. 
The  Figure  is  divided  into  two  parts  to  describe  the  differences  between  "crack  geometry" 
and  observed  "fringes"  patterns.  The  cross-hatched  area  illustrates  that  the  two  mating  crack 
surfaces  are  closed  (in  contact).  The  clear  area  (not  cross-hatched)  illustrates  that  the  crack 
surfaces  are  separated,  and  is  referred  to  in  this  report  as  a  geometric  "void."  The  dotted 
"fringe  order"  lines  represent  contours  of  constant  crack  opening  displacement  (COD) 
patterns. 

The  terminology  "void"  will  be  used  in  the  remainder  of  the  report  to  describe  a 
displacement  hump  which  is  internal  to  the  crack  mating  surfaces.  This  internal  void 
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displacement  is  always  open  once  formed  even  under  zero  load.  The  void  formation  is 
depicted  in  Figure  5.1  by  the  fringe  orders,  and  shown  by  interference  photos  in  Figure 
4.3a-d. 

A  type  I  crack  (Figure  5.2a)  is  defined  as  the  crack  geometry  where  the  two  mating 
surfaces  are  fully  closed  under  zero  load.  This  crack  type  is  typically  modeled  as  an  elastic 
COD  pattern.  The  type  I  crack  was  observed  for  smaller  cracks  having  an  a/t  value  less  than 
0.095  in  this  investigation. 

As  the  crack  grew  in  size  due  to  fatigue  cycling,  a  type  II  crack  was  formed  (Figure 
5.2b).  Hiis  change  in  crack  geometry  occurs  when  the  crack  internal  surfaces  separate  and 
form  a  displacement  "void”  close  to  the  crack  tip  at  location  A.  The  pattern  is  distinguished 
by  the  interference  fringe  order  contours  under  zero  load.  The  uniqueness  of  a  type  11  crack 
is  that  under  zero  load  the  crack  tip  boundary  is  fully  closed,  while  there  is  an  internal 
nonzero  displacement  void  area  separating  the  two  crack  surfaces.  If  the  void  displacement  is 
too  small  to  be  detected  by  a  fringe  pattern  (the  void  maximum  displacement  is  less  than 
6.2E-6  inches),  the  crack  initially  opens  at  this  location  when  load  is  applied  (as  shown  by 
the  fringe  order  formation).  In  this  investigation,  the  formation  of  a  type  II  crack  ranged  in 
size  from  an  a/t  of  0.068  to  0.095. 

As  the  crack  continues  to  grow  due  to  fatigue  cycling,  the  void  size  increases  in  area 
and  height.  This  void  growth  forces  the  crack  tip  at  location  A  to  separate  and  remain  open 
even  under  zero  applied  load.  This  opening  of  the  crack  tip  defines  a  type  III  crack  (Figure 
5.2c).  In  this  investigation,  a  type  III  crack  was  formed  when  the  a/t  range  was  between 
0.139  and  0.189  for  all  tests. 

Another  means  of  visualizing  the  three  void  types  is  by  slicing  the  crack  between 
locations  A  and  B  («^  =  90  degrees.  Figure  5.1)  as  schematically  depicted  in  Figure  5.3.  For  a 
type  I  crack  (Figure  5.3a),  the  crack  mating  sur&ces  are  completely  closed  under  zero  load. 
With  increased  crack  growth,  a  type  II  crack  develops  when  the  crack  surfaces  separate 
(forming  an  internal  void  area)  while  the  crack  tip  boundary  remains  closed  (Figure  5.3b). 
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TYPE  I 


As  the  crack  continues  to  grow,  a  type  III  crack  is  formed  when  the  void  size  is  sufficiently 
large  to  open  the  crack  tip  at  location  A  (Figure  5.3c). 

The  percent  of  crack  tip  boundary  AD  which  is  open  for  a  type  III  crack  at  zero  load 
depends  on  the  crack  size.  For  example,  the  crack  tip  boundary  between  locations  A  and  D  is 
open  17%  for  a  newly  formed  type  III  crack  (a/t  =  0.156  and  a/c  =  0.652  from  the  constant 
AK^j  test).  When  the  crack  reached  the  size  of  a/t  =  0.180  and  a/c  =  0.677,  30%  of  the  AD 
boundary  is  open.  With  continued  fatigue  cycling,  65%  of  the  AD  boundary  is  open  for  a 
crack  size  of  a/t  =  0.260  and  a/c  =  0.650.  This  opening  of  the  crack  tip  boundary  between 
locations  A  and  D  results  from  an  increase  in  void  surface  area  and  internal  displacement 
with  crack  growth. 

Table  5.1  lists  the  measured  crack  lengths  when  transition  occurred  between  crack 
types  I  and  II,  and  type  II  and  III  cracks.  The  crack  transition  size  was  found  to  correlate 
most  favorably  with  a/t.  Table  5.2  presents  the  measured  range  in  crack  size  for  the  three 
crack  types.  These  ranges  do  not  include  precracking  measurements  which  show  higher  a/c 
values  for  the  smaller  crack  sizes  (Appendix  E  for  a  complete  data  set). 

Table  5.1:  Crack  Size  when  Transition  Occurs  Between  a  Crack  Types  I  and  II,  and  Crack 
Types  II  and  III. 


Transition  to 

a 

c 

a/c 

a/t 

Type  II  Crack 

Constant 

0.051 

0.058 

0.872 

0.068 

Constant  AK^ 

0.060 

0.075 

0.840 

0.084 

Constant  Load 

0.071 

0.084 

0.845 

0.095 

Block  Load 

0.052 

0.075 

0.693 

0.069 

Type  III  Crack 

Constant  AK^ 

0.117 

0.180 

0.652 

0.156 

Constant  AK^ 

0.127 

0.153 

0.833 

0.169 

Constant  Load 

0.104 

0.149 

0.698 

0.139 

Block  Load 

0.142 

0.202 

0.705 

0.189 
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Table  5.2:  Range  of  Measured  Crack  Sizes  for  the  Three  Crack  Types 


a 

c 

a/c 

Type  I  Crack 

Constant  AK^ 

0.043-0.048 

0.045-0.056 

0.94-0.84 

Constant  AK^ 

0.034-0.055 

0.038-0.062 

0.88-0.71 

Constant  Load 

0.040-0.063 

0.043-0.074 

0.97-0.85 

Block  Load 

0.049-0.050 

0.065-0.071 

0.75-0.71 

Type  II  Crack 

Constant  AK^ 

0.050-0.116 

0.058-0.174 

0.87-0.67 

Constant  AlCj 

0.063-0.127 

0.075-0.148 

0.85-0.81 

Constant  Load 

0.071-0.098 

0.084-0.138 

0.84-0.69 

Block  Load 

0.052-0.135 

0.075-0.193 

0.70-0.64 

Type  III  Crack 

Constant  AK^ 

0.117-0.325 

0.180-0.650 

0.65-0.50 

Constant  AlCj 

0.127-0.292 

0.152-0.651 

0.83-0.45 

Constant  Load 

0.104-0.343 

0.149-0.826 

0.67-0.42 

Block  Load 

0.142-0.357 

0.202-0.651 

0.71-0.55 

A  general  trend  in  crack  aspect  ratio  (a/c)  is  noted  in  Table  5.2,  where  a/c  for  a  type  I 
crack  is  higher  than  a  type  II,  and  a/c  for  a  type  II  crack  is  higher  than  a  type  III.  Although  it 
can  be  reasoned  that  this  trend  is  caused  by  the  stress  being  larger  at  the  crack  free-surface 
than  at  the  depth  position  due  to  the  bending,  this  is  not  the  only  consideration.  First,  the 
stress  intensity  factor  into  the  crack  depth  (Kg,  location  A)  can  be  greater  than  at  the  crack 
free-surface  (K^j,  location  D)  for  small  cracks  (depending  on  the  crack  a/t  and  a/c  growth 
pattern).  For  the  constant  AK^  and  block  loading  tests,  AK^  was  larger  than  ABC^j  during 
precracking  when  a/t  was  less  than  0.040  and  0.061,  respectively.  Second,  even  though  the 
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crack  initiation  sight  a/c  was  0.17  (which  would  promote  crack  growth  along  the  free-surface 
if  stress  were  the  major  factor),  the  crack  did  not  grow  along  this  boundary,  but  initiated  a 
new  site  such  that  a/c  was  approximately  1.0  (that  is,  semicircular,  ref.  Figure  4.17),  and  then 
decreased  with  crack  growth. 

This  crack  growth  pattern  (where  a/c  is  initially  equal  to  or  greater  than  one  and  then 
decreases  with  crack  growth)  has  been  demonstrated  both  analytically  (5)  and  experimentally 
(40,42)  in  metals  for  both  bending  and  tension.  The  difficulty  of  growing  cracks  to  desired 
a/c  values  in  metals  was  also  described  by  Fleck,  et  al.  (18)  when  he  attempted  to  force  a 
crack  a/c  by  different  initiation  site  shapes.  He  stated  that  "By  employing  two  different 
shapes  of  slot  we  hoped  to  investigate  the  effect  of  crack  shape  on  closure  response. 
Unfortunately,  the  fatigue  cracks  quickly  grew  to  similar  a/c  ratios,  precluding  a  study  of 
crack  shape  effects." 


COD  Profiles  with  Applied  Load 

The  change  in  COD  patterns  as  load  is  applied  to  a  surface  crack  is  a  direct  function 
of  crack  type  (Figure  5.2).  For  a  type  I  crack  (whose  surfaces  are  fully  closed  under  zero 
load.  Figure  5.2a),  the  COD  pattern  is  classical  in  that  it  begins  to  open  at  the  crack  free- 
surface  (location  B),  and  continues  opening  symmetrically  toward  the  crack  tip  boundary. 
Although  the  complete  crack  tip  boundary  A  to  D  opens  at  nearly  the  same  load  level, 
location  A  was  observed  to  open  just  prior  to  location  D.  The  elastic  opening  pattern  with 
applied  load  forms  a  fringe  pattern  schematically  illustrated  in  Figure  5.4e. 

The  COD  patterns  during  load  application  for  crack  types  II  and  III  are  quite  different 
from  a  type  I  crack  as  shown  in  Figure  5.4.  Note  that  Figure  5.4  is  organized  into  "crack 
geometry"  and  "fringes"  groupings  (similar  to  Figure  5.2).  The  changes  between  Figures 
5.4a  and  5.4e  are  caused  by  increasing  the  applied  load.  Beginning  with  a  type  II  crack 
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(Figure  S.4a)  and  applying  a  load,  the  internal  crack  displacement  void  increases  in  area  and 
displacement  until  the  crack  tip  separates  at  location  A  (Figure  5.4b).  As  load  continues  to 
increase,  the  void  area  expands  to  the  specimen  hree-surface  at  location  B  (Figure  S.4c).  The 
crack  hree-surface  then  begins  to  open,  and  the  crack  internal  surfaces  continue  to  separate 
(Figure  5.4d)  until  the  crack  tip  boundary  AD  is  completely  open  (Figure  5.4e). 

The  type  III  CX)D  pattern  follows  the  same  general  trend  as  a  type  II,  but  begins  at 
Figure  5.4b  and  progresses  to  Figure  5.4e  upon  load  application.  It  should  be  noted  that  the 
fringe  patterns  in  Figures  5.4a-d  show  the  crack  maximum  displacement  to  be  internal  to  the 
crack,  whereas  it  is  at  the  crack  free-surface  in  Figure  5.4e.  These  patterns  are  described  in 
more  detail  using  COD  measurements  in  the  next  section. 

Photographic  examples  of  the  fringe  order  patterns  are  shown  in  Figures  5.5a  and 
5.5b  (reference  Figures  4.4a-d).  Note  the  crack  tip  at  location  A  is  open.  Figure  5.5a  is  an 
interferometric  photograph  of  a  type  III  crack  with  zero  applied  load  (a  =  0.195  inch  c  - 
0.300  inch).  Figure  5.5b  shows  a  type  II  crack  with  an  applied  load  level  sufficient  to 
completely  open  the  crack  surfaces  (a  =  0.168  inch  c  =  0.249  inch).  The  dark  interference 
fringe  patterns  can  be  seen  in  each  photo. 

Figure  5.5b  fringe  orders  shows  a  fully  open  crack  with  the  maximum  displacement 
internal  to  the  crack  surfaces  (i.e.,  COD  at  the  free>surface  is  less  than  the  internal 
displacement,  reference  Figure  5.4d).  The  mating  surfaces  at  location  C  are  noted  to  be  the 
last  portion  of  the  crack  to  open.  Location  C  ranged  from  12  to  15  degrees  from  the  crack 
free-surface  (angel  Figure  5.1).  This  condition  is  speculated  to  be  caused  by  an  apparent 
transition  from  a  state  of  plane  stress  to  one  of  plane  strain. 

A  three-dimensional  elastic-plastic  frnite  element  analysis  of  a  surface  crack  by 
Trantina,  deLorenze  and  Wilkening  (43)  clearly  shows  a  plane  stress/strain  transition  zone  at 
=  15  degrees.  This  effect  was  not  observed  in  the  elastic  analysis,  but  was  prevalent  in  the 
elastic-plastic  predictions.  The  paper  describes  this  variation  of  K  along  the  crack  boundary 
as  "...a  combination  of  increased  K  due  to  plasticity  effects  and  a  loss  of  plane  strain 
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Figure  5.5  Interference  photograph  of  the  crock  surface  displacement 
at  (a)  near  zero  load  snowinK  the  internal  void  formation  for 
a  ni  crack  (refvfigure  5.4b),  and  (b)  for  a  fully-open 
cr^  in  the  dosureloaid  region  showing  the  minimum 
di^lacement  at  location  C,  (|>=:12  to  15  degrees  (ref.  tigdre 
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constraint  at  the  free  surface."  Their  analysis  showed  a  rapid  decrease  in  the  plane  strain 
constraint  o-^lip2  +  at  <|>  =  15  degrees,  causing  "...increased  yielding  at  the  free  surface 
and  a  decrease  in  effective  K.  The  region  below  the  free  surface  carries  a  larger  portion  of 
the  load  and  thus  a  peak  in  the  K  variation  with  develops  at  about  15  (|i."  The  effects  of 
more  plasticity  induced  from  a  higher  strain  hardening  exponent  or  higher  strain  level  were 
found  to  reduce  the  plane  strain  constraint  at  the  free  surface  and  decrease  the  K  effectivity 
for  <^  <  15  degrees. 


COD  Variation  with  Crack  Size 

Fatigue  crack  size  is  a  dominant  factor  in  categorizing  the  three  crack  types,  and  their 
corresponding  crack  opening  displacement  (COD)  patterns.  The  question  of  crack  size  when 
the  void  forms  a  type  II  crack,  and  then  a  type  III  crack,  was  covered  by  data  presented  in 
Table  5.1.  The  question  of  void  changes  with  increasing  crack  size  and  applied  load  is 
presented  in  this  section.  First,  data  are  presented  to  track  void  growth  patterns  at  zero  load. 
Second,  data  to  track  the  void  changes  with  an  applied  load  are  presented. 

Organization  of  this  information  can  be  described  best  by  taking  a  cross-sectional 
slice  between  locations  A  and  B  (Figure  5.1,  <^  =  90  degrees)  to  visualize  the  crack  opening 
displacement  profile  as  a  function  of  applied  load  ~  shown  in  Figure  5.6.  The  type  II  crack 
data  are  from  the  constant  test  (a  =  0.12  and  c  =  0.17  inch),  and  represent  half  the  COD 
profile  sketch  shown  in  Figure  5.3b.  The  two  diamonds  located  along  the  abscissa  represent 
a  zero  fringe  order,  or  zero  displacement.  The  dotted  line  connecting  the  two  measured  zero 
fringe  order  points  is  a  projected  estimate  of  the  crack  surface  profile. 

Note  that  for  low  load  levels  the  void  maximum  displacement  is  located  internal  to 
the  crack  surfaces,  and  is  closer  to  the  crack  tip  (location  A)  than  the  crack  free-surface 
(location  B).  With  an  applied  load,  the  void  maximum  displacement  moves  away  from  the 
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location  A  crack  tip  and  toward  the  crack  free-surface  at  location  B.  The  applied  load  level 
where  the  void  maximum  displacement  transitions  hrom  being  internal,  to  being  at  the  crack 
free-surface,  is  deflned  as  tPgp  (K  value  of  187  psi*in^^^  in  Figure  5.6). 

To  track  the  crack  internal  void  displacement  and  movement  with  crack  size  and 
applied  load,  three  variables  have  been  established.  These  variables  are  schematically 
illustrated  in  Figures  S.7a  and  5.7b  such  that: 

(1)  5  defines  the  crack  maximum  fringe  order  (h,  ref  eqn.  4.1)  measured  for  a 
given  applied  load. 

(2)  and  A2  deflne  the  measured  distances  from  the  crack  tip  (along  a  line 
between  locations  A  and  B)  to  the  fringe  orders  which  encompass  the  void 
maximum  displacement.  That  is,  the  displacement  6,  and  A^  and  A2 
distances  from  the  crack  tip,  locate  the  COD  maximur”  int  for  a  given 
applied  load. 


These  definitions  are  used  in  Tables  5.3  and  5.4  lo  describe  void  COD  patterns  due  to 
crack  size  and  applied  load.  Table  5.3  presents  COD  patterns  under  zero  load  for 
representative  crack  sizes. 


Table  5.3:  Constant  AK^j  Test  Effects  on  the  Void  Maximum  Displacement  Location  with 
Zero  Load 


Crack  Size 

Applied 

Fringe  Order 

avoid 

Moment 

Distance  from 

Displacement 

a 

c 

nn-lbl 

Crack  Tin  fin'l 

ffrinee  orderl 

4l  Ai 

0.09 

0.15 

0 

— 

1 

0.14 

0.20 

0 

0.010  0.034 

1 

0.20 

0.30 

0 

0.033  0.056 

6 

0.28 

0.50 

0 

0.085  0.098 

14 

0.30 

0.55 

0 

0.080  0.113 

15 
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A  plot  of  the  Table  5.3  data  is  presented  in  Figure  5.8  to  show  relative  change  in  void 
displacement  with  increasing  crack  size.  Again,  measurements  are  made  at  zero  load,  and 
the  and  A2  distances  from  the  crack  tip  bound  the  void  maximum  displacement  point. 
Note  that  the  void  maximum  displacement  point  increases  in  height  and  moves  away  from 
the  crack  tip  and  toward  the  crack  free-surface  as  the  crack  increases  in  size.  The  crack  tip 
location  is  also  plotted  to  show  relative  effecte  with  crack  growth. 

Table  5.4  presents  the  same  type  of  data  as  Table  5.3,  but  for  several  values  of 
applied  load.  The  void  maximum  displacement  data  from  Table  5.4  are  plotted  for  two  crack 
sizes  in  Figure  5.9  to  show  general  trends  in  void  movement  with  applied  load.  Note  in 
Figure  5.9  that  the  void  position  is  plotted  as  a  function  of  the  crack  tip  instead  of  from  the 
crack  free-surface  as  in  Figure  5.8. 
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Table  5.4:  Void  Maximum  Displacement  LxKation  for  a  Constant  AK^j  Test  Effects  with  an 
Applied  Load 


Crack  Size 

Applied 

Fringe  Order 

6  -  Void  Max 

Moment 

Distance  from 

Displacement 

a 

c 

fin-lbl 

Crack  Tip  finl 

ffrinee  orderl 

Al 

Al 

0.09 

0.15 

72 

0.017 

0.030 

2 

0.09 

0.15 

88 

0.020 

0.027 

3 

0.09 

0.15 

102 

0.023 

0.039 

4 

0.09 

0.15 

132 

0.041 

0.048 

8 

0.09 

0.15 

163 

0.038 

0.048 

10 

0.14 

0.20 

18 

0.016 

0.023 

2 

0.14 

0.20 

88 

0.023 

0.050 

6 

0.14 

0.20 

102 

0.038 

0.049 

9 

0.14 

0.20 

132 

0.038 

0.094 

13 

0.20 

0.30 

18 

0.034 

0.057 

7 

0.20 

0.30 

51 

0.045 

0.059 

10 

0.20 

0.30 

88 

0.045 

0.082 

16 

0.20 

0.30 

102 

0.058 

0.082 

19 

0.20 

0.30 

32 

0.082 

0.105 

24 

0.30 

0.55 

18 

0.084 

0.115 

17 

0.30 

0.55 

51 

0.104 

0.116 

23 

0.30 

0.55 

72 

0.113 

0.142 

28 

0.30 

0.55 

88 

0.119 

0.153 

33 

95 


Three  primary  observations  can  be  made  about  the  void  maximum  displacement  point 
brom  these  data: 

(1)  The  void  maximum  displacement  location  is  closer  to  the  crack  tip  (location 
A)  for  smaller  cracks. 

(2)  The  void  height  is  lower  for  smaller  cracks  than  for  larger  cracks  at  an 
equivalent  applied  load. 

(3)  As  load  is  applied  to  the  specimen,  the  void  height  increases,  and  moves 
away  from  the  crack  tip  and  toward  the  crack  free-surface  —  eventually 
opening . 


COD  Variation  with  Load  History 

The  previous  section  described  the  void  movement  with  increased  crack  size  and  with 
an  applied  load  for  the  constant  AK^  test.  This  section  will  expand  this  data  set  to  include 
experimental  results  from  the  other  tests.  The  same  definitions  of  6,  Aj  and  A2  will  be  used 
in  tables  5.5  through  5.7  to  cover  the  constant  AKa,  constant  load  and  block  loading  tests. 


97 


Table  5.5:  Void  Maximum  Displacement  Location  for  a  Constant  AK^  Test  Effects  with  an 
Applied  Load 


Crack  Size 

Applied 

Fringe  Order 

6  -  Void  Max 

Moment 

Distance  from 

Displacement 

a 

c 

nn-lbl 

Crack  Tio  ('in'l 

ffrinee  orderl 

Ai 

4l 

0.10 

0.12 

24 

0.019 

0.031 

2 

0.10 

0.12 

36 

0.011 

0.047 

2 

0.10 

0.12 

57 

0.013 

0.058 

3 

0.10 

0.12 

73 

0.021 

0.065 

5 

0.10 

0.12 

90 

0.033 

0.082 

7 

0.13 

0.15 

0 

0.012 

0.045 

2 

0.13 

0.15 

24 

0.017 

0.045 

3 

0.13 

0.15 

36 

0.029 

0.036 

4 

0.13 

0.15 

57 

0.028 

0.050 

6 

0.13 

0.15 

90 

0.058 

0.088 

11 

0.14 

0.20 

23 

0.022 

0.060 

5 

0.14 

0.20 

36 

0.026 

0.059 

6 

0.14 

0.20 

57 

0.044 

0.059 

10 

0.14 

0.20 

73 

0.051 

0.077 

13 

0.19 

0.30 

0 

0.040 

0.089 

6 

0.19 

0.30 

2^ 

0.052 

0.091 

9 

0.19 

0.30 

57 

0.072 

0.098 

17 

0.19 

0.30 

73 

0.079 

0.119 

21 

0.19 

0.30 

81 

0.084 

0.150 

23 

0.20 

0.32 

0 

0.050 

0.091 

7 

0.20 

0.32 

23 

0.064 

0.093 

10 

0.20 

0.32 

57 

0.083 

0.107 

19 

0.20 

0.32 

81 

0.092 

0.116 

24 

0.20 

0.32 

98 

0.107 

0.122 

30 
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Table  5.6:  Void  Maximum  Displacement  Location  for  a  Constant  Load  Test  with  an 
Applied  Load 


Crack  Size 

Applied 

Moment 

Fringe  Order 
Distance  from 

6  -  Void  Max 

Displacement 

a 

c 

Hn-lb) 

Crack  Tin  fin'l 

Al  4l 

^fringe  orders 

0.08 

0.10 

132 

0.007 

0.035 

1 

0.08 

0.10 

138 

0.018 

0.026 

2 

0.08 

0.10 

154 

0.022 

0.047 

3 

0.08 

0.10 

162 

0.029 

0.047 

4 

0.09 

0.12 

100 

0.007 

0.031 

1 

0.09 

0.12 

114 

0.012 

0.024 

2 

0.09 

0.12 

132 

0.016 

0.027 

3 

0.09 

0.12 

138 

0.015 

0.036 

3 

0.09 

0.12 

154 

0.022 

0.033 

5 

0.09 

0.12 

163 

0.028 

0.035 

6 

0.09 

0.12 

195 

0.036 

0.067 

9 

0.10 

0.15 

100 

0.018 

0.027 

3 

0.10 

0.15 

114 

0.012 

0.042 

3 

0.10 

0.15 

132 

0.021 

0.036 

5 

0.10 

0.15 

138 

0.015 

0.050 

5 

0.10 

0.15 

154 

0.020 

0.050 

7 

0.10 

0.15 

163 

0.029 

0.049 

9 

0.10 

0.15 

195 

0.051 

0.082 

14 

0.12 

0.17 

46 

0.016 

0.042 

1 

0.12 

0.17 

100 

0.021 

0.031 

5 

0.12 

0.17 

114 

0.026 

0.035 

6 

0.12 

0.17 

132 

0.021 

0.047 

7 

0.12 

0.17 

138 

0.027 

0.046 

8 

0.12 

0.17 

163 

0.032 

0.055 

11 

0.12 

0.17 

195 

0.052 

0.088 

17 

0.14 

0.20 

46 

0.019 

0.039 

3 

0.14 

0.20 

82 

0.023 

0.042 

5 
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Table  5.7:  Void  Maximum  Displacement  Location  for  a  Block  Loading  Test  with  an 
Applied  Load 


Crack  Size 

Applied 

Fringe  Order 

6  -  Void  Max 

Moment 

Distance  from 

Displacement 

a 

c 

fin-lbl 

Crack  Tip 

lin) 

ffrince  order! 

4i 

4l 

0.25 

0.35 

0 

0.060 

0.105 

8 

0.25 

0.35 

23 

0.073 

0.110 

11 

0.25 

0.35 

32 

0.081 

0.107 

12 

0.25 

0.35 

47 

0.085 

0.100 

16 

0.25 

0.35 

57 

0.088 

0.099 

19 

0.25 

0.35 

73 

0.101 

0.128 

23 

0.25 

0.35 

81 

0.110 

0.120 

25 

0.25 

0.35 

98 

0.111 

0.159 

28 

To  compare  the  void  displacement  patterns  for  all  tests,  it  is  desirable  to  select  a 
common  crack  aspect  ratio  (a/c)  and  a/t  value.  TTiis  point  is  shown  in  Figure  5.10  where  a/c 
versus  a/t  is  plotted  for  all  four  tests.  The  common  point  selected  from  this  Figure  was  a/c  = 
0.69  and  a/t  =  0.19.  Comparison  of  the  void  displacement  patterns  are  shown  in  Figure  5.11 
with  applied  load  for  a  common  crack  size.  The  reference  line  segment  is  drawn  for  the 
constant  AK^  data,  and  is  found  to  fit  the  other  data  sets  quite  well. 

It  can  be  seen  from  this  plot  that  the  void  movement  toward  the  crack  free-sur&ce 
with  applied  load  (as  discussed  in  the  last  section,  ref.  Figure  5.9)  is  consistent  for  various 
load  histories.  However,  the  applied  moment  required  to  achieve  a  given  displacement 
(shown  in  parentheses)  varies  with  load  history.  That  is,  to  achieve  a  void  maximum 
displacement  of  13  fringe  orders  requires  an  applied  moment  of  132  in-lb.  for  the  constant 
AK(j  test,  but  only  73  in-lb.  for  the  constant  AKg  test  These  differences  are  also  seen  in 
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Figure  5.12,  where  the  void  maximum  displacement  is  plotted  as  a  function  of  applied  load. 
Here,  the  constant  load  test  is  shown  to  require  a  higher  applied  load  to  achieve  a  given  level 
of  COD  than  either  a  constant  AK^  or  AK^  test.  Also,  a  constant  AK<]  test  requires  a  higher 
applied  load  to  achieve  a  given  displacement  than  does  a  constant  AK^  test.  These 
differences  are  a  clear  indication  of  load  history  effects  on  closure  load,  and  will  be 
discussed  in  the  next  section  on  "Closure/Opening  Loads.” 

In  all  tests,  three  primary  observations  of  the  void  COD  pattern  were  noted: 

(1)  Under  zero  load,  the  void  maximum  displacement  is  a  direct  function  of 
crack  size 

(2)  Under  zero  load,  the  void  maximum  displacement  location  relative  to  the 
crack  tip  (location  A)  is  a  direct  function  of  crack  size 

(3)  With  applied  load,  the  void  maximum  displacement  increases,  and  moves 
away  from  the  crack  tip  and  toward  the  crack  free-surface  (location  B). 

Crack  Compliance 

Using  the  crack  opening  displacement  profile  data  from  Figure  5.6,  load  versus 
displacement  plots  (typically  used  to  define  closure  load  when  measurements  are  made  at  the 
crack  free-surface)  can  be  constructed  for  different  distances  from  the  crack  tip.  These  load- 
displacement  plots  are  presented  in  Figure  5.13  for  points  along  a  line  between  location  A 
and  B,  =  90  degrees.  The  inverse  slope  of  the  load  versus  displacement  graph  is  defined  as 
the  compliance  at  a  particular  point,  and  represents  the  crack  opening  displacement  per  unit 
load. 

For  a  linear-elastic  body,  compliance  at  a  point  is  independent  of  load  level,  and  a 
load-displacement  curve  would,  therefore,  be  linear.  Any  deviation  from  linearity  in  a 
linear-elastic  body,  such  as  those  due  to  closure,  could  be  attributed  to  changes  in  the 
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geometry  of  the  body.  This  nonlinearity  is  seen  in  Figure  5.13,  where  compliance  changes 
until  the  load  is  increased  above  a  certain  value.  This  nonlinear  behavior  in  the  load- 
displacement  plot  is  due  to  changes  in  crack  shape  during  the  loading  cycle.  The  compliance 
nonlinearity,  caused  by  changes  in  crack  geometry,  also  indicates  that  the  stress  intensity 
factor  at  any  given  point  along  the  crack  boundary  is  also  nonlinear  with  applied  load  (44). 

Besides  the  nonlinearity  of  the  load-displacement  profile,  compliance  is  also  noted  to 
be  dependent  on  measurement  location.  This  is  seen  in  Figure  5.13  where  the  compliance,  or 
the  inverse  slope  of  the  load-displacement  trace,  is  approximately  4  times  greater  at  the  crack 
firee-surface  (0.2  in  from  the  crack  tip)  than  at  a  distance  of  0.05  in  from  the  crack  tip  for 
loads  greater  than  80  in-lb.  As  the  crack  grows  in  size  due  to  fatigue  cycling,  compliance  is 
shown  to  increase  at  the  crack  free-sur£ace  (Figure  5.14).  This  pattern  follows  typical  trends 
for  thru-thickness  cracks  which  increase  in  flexibility  or  decrease  in  rigidity  with  increasing 
crack  size  (12,27). 

Closure/Opening  Loads 

The  experimental  investigation  identified  three  primary  closure  loads  which  are 
considered  important  for  the  understanding  of  fatigue  crack  growth.  How  these  closure  loads 
correlate  with  the  three  crack  types  described  above  forms  the  basis  for  relating  the 
geometric  crack  (Opening  displacement  (COD)  patterns  with  the  residual  induced  closure 
mechanisms.  This  section  will  discuss  these  relationships  by  defining  three  unique 
closure/opening  loads  measurements,  and  then  relating  them  to  the  three  COD  types,  crack 
size  and  fatigue  load  history. 


CRACK  DISPLACEMENT  MEASUREMENTS 


Figure  5.14  Load  venus  displacement  plot  at  the  crack  face  (location  B, 
4°^  degrees)  as  a  function  of  crack  size. 


Closure  Load  Definitions 


The  identification  of  three  distinct  closure  loads  raises  several  questions  concerning 
their  relationship  to  the  standard  definition  of  closure  load  (deflned  from  a  load  versus 
displacement  plots  as  described  in  Chapter  IV).  In  reality,  the  potential  for  having  more  than 
one  closure  load  for  the  same  crack  size  is  not  hypothetical.  For  example,  the  load- 
displacement  measurements  described  in  Figure  S.13  identified  the  influence  of  measurement 
location  on  closure  load.  As  can  be  seen,  several  closure  load  values  can  also  be  generated 
using  the  standard  definition  of  closure  if  different  measurement  locations  are  used.  This 
sensitivity  to  measurement  location  has  been  discussed  by  several  investigators  using  clip 
gages,  back  face  strain  and  interferometry  techniques  (11,12,24).  These  considerations  will 
be  evaltiated  more  fully  in  the  next  chapter  when  the  definitions  of  AKeff  are  considered. 

Figure  5.15  is  used  to  define  the  three  primary  closure  loads  (aPop,  bPop  and  cPop) 
of  this  investigation.  The  aPop,  bPgp  and  cPop  closure  loads  are  defined  here  as  the  loads 
needed  to  open  the  crack  surfaces  (Pop)  at  the  A,  B  and  C  locations,  respectively.  Location 
C  is  the  last  portion  of  the  crack  surface  to  completely  open  under  an  applied  load,  and  is 
typically  located  12  to  15  degrees  into  the  specimen  thickness  from  the  crack  free-surface. 
Figure  5.15  also  illustrates  typical  interference  fringe  order  patterns,  or  lines  of  constant 
displacement,  where  the  maximum  displacement  (highest  fringe  order)  is  internal  to  the 
crack  surfaces.  The  dotted  line  represents  a  section  slice  between  locations  A  and  B,  and  is 
used  to  illustrate  COD  profiles  as  shown  in  Figure  5.16. 

Correlation  of  the  bPop  and  cPgp  closure  load  measurements  with  the  geometric 
COD  profiles  is  shown  in  Figure  5.16.  This  information  was  collected  independently,  with 
the  closure  load  measurements  being  made  directly  from  the  fatigue  machine  load  cell 
readings,  the  COD  proriles  generated  using  35  mm  negatives  and  a  digitizing  board 
(described  in  Chapter  IV).  This  information  is  combined  in  Figure  5.16  by  projecting  the 
COD  profile  at  the  same  P/Pmax  ^  ^^op/^max  closure  load  ratio  of  0.22.  The 
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FREE  SURFACE 


I»ofile  projection  at  this  load  level  is  noted  to  intersect  at  the  crack  face  (free-surface),  thus 
correlating  the  definition  of  the  bPop  -*  the  load  level  which  begin  to  open  the  crack  free- 
surfrice.  Thus,  very  good  agreement  is  noted  between  the  COD  profiles  and  the  bPop  closure 
load  measurement  Projection  of  the  crack  profile  at  the  cPnp/PTpay  =  0.28  load  level 
corresponds  to  the  crack  free-surface  at  location  B  being  open  approximately  8.6E-S  inches. 
The  tPop  transition  load  is  defined  as  the  level  needed  to  move  the  void  maximum 
displacement  from  the  crack  interior  to  the  crack  free-surface.  For  a  load  ratio  of  P/Pmax  ~ 
0.30,  the  crack  maximum  displacement  is  approximately  12.8E-5  inches  at  location  B.  For 
loads  above  tPop,  the  COD  is  no  longer  in  the  closure  load  region,  and  follows  a  classical 
linear-elastic  pattern. 

Figure  5.17  schematically  illustrates  the  relationships  between  the  three  closure  loads 
(aPop,  bPgp  and  cP^p)  and  the  three  crack  opening  displacement  types  (types  I,  II  and  ni). 
For  a  type  I  crack  (Figure  S.17a),  the  crack  surfaces  are  closed  under  zero  load.  As  load  is 
applied  the  crack  remains  closed  until  a  load  level  of  bP^p  is  reached  and  the  crack  begins  to 
open  at  the  crack  free-surface  (location  B).  With  increased  load,  the  crack  opens 
symmetrically  inward  and  is  fully  open  at  locations  A  and  C  at  the  same  load  level  (aP^p  = 
cPop)*  Since  the  maximum  displacement  is  always  at  the  crack  free-surface,  tPgp  =  bPop, 
and  cPop  =  aPop  »  bPop  »  0. 

For  a  type  II  crack,  a  void  formation  has  developed  internal  to  the  crack  surfaces  as 
shown  for  zero  'oad  in  Figure  5.17b.  The  crack  boundary  and  free-surface  is  closed 
(locations  A,  B  and  C)  under  zero  load.  With  load  application,  location  A  opens  at  a  load 
level  of  aPop.  At  the  bPop  load  level,  the  crack  free-surface  begins  to  open  at  location  B  and 
is  completely  open  between  A  and  B  (^  =  90  degrees);  however,  the  crack  tip  is  closed  at 
location  C.  At  the  cPop  load  level,  the  crack  surfrices  are  completely  open  with  the 
maximum  displacement  along  the  ^  =  90  degree  line.  For  a  type  n  crack,  tPop  »  cPop  » 
bPop  ®^Op’ 
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For  a  type  III  crack,  location  A  is  always  open  while  locations  B  and  C  are  closed 
under  zero  load.  Therefore,  the  aPop  load  level  is  zero  for  a  type  III  crack.  Location  B 
begins  to  open  at  the  bPop  load  level,  and  at  cPop  the  crack  is  completely  open.  That  is, 
aPop=0  and  tpQp  »  cPop  »  bPop. 

Closure  Load  Variations  with  Crack  Size 

Closure  load  variations  with  crack  growth  are  found  to  be  directly  related  to  crack 
size.  This  is  seen  in  Figure  5.18,  where  the  three  visually  observed  closure  load  absolute 
values  are  plotted  as  a  function  of  crack  size  for  the  constant  AK^  test.  The  inverse 
relationship  of  decreasing  closure  load  with  increasing  crack  size  is  evident  in  this  Figure. 
For  a  type  I  crack,  apQp  and  cPqp  are  at  approximately  the  same  load  level,  which  is  greater 
than  bPop.  With  crack  growth,  a  type  II  crack  is  developed  due  to  the  formation  of  a  void 
close  to  the  location  A  crack  tip.  This  void  formation  causes  the  crack  tip  at  location  A  to 
open  before  the  crack  free-surface  at  location  B.  Although  aPop  is  initially  larger  than  bPop, 
the  void  formation  causes  aPop  to  decrease  at  a  faster  rate  than  bPop  with  respect  to  crack 
growth.  With  continued  crack  growth,  the  void  becomes  sufficiently  large  to  physically  open 
the  crack  tip  at  location  A,  forming  a  type  III  crack. 

Crack  sizes  (a/t)  where  transition  occurs  between  a  type  I,  n  and  III  crack  were 
determined  from  visual  observations  and  analysis  of  the  35  mm  film.  When  this  information 
is  plotted  independently  on  Figure  5.18  (Reference.  Table  5.1),  and  dashed  lines  extended 
from  the  a/t  transition  points  through  the  experimental  data  sets,  definite  changes  in  slope  are 
observed  in  the  bPop  and  cPqp  data  sets  at  the  transition  between  type  II  to  type  III  crack. 
The  aPop  closure  loads  follow  a  continually  decreasing  trend  with  increasing  crack  size,  and 
is  zero  for  a  type  III  crack. 
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Figure  5.19  presents  the  closure  load  ratio  (Pop/Pmax)  ^  ^  function  of  crack  size  for 
aPop,  bPop  and  cPop.  The  cPop  closure  ratio  (cPop/Pmax)  increase  to  an  average 

value  of  0.29  during  the  formation  of  type  I  and  II  cracks.  With  the  formation  of  a  type  III 
crack  (a/t  =  0.156),  cPop/Pj^ax  gradually  increases  and  leveJs-out  at  around  0.36.  These 
values  of  closure  ratio  correspond  to  the  general  range  of  published  values  summarized  by 
Schijve  (45)  for  mild  steel,  titanium  6-4  and  aluminum  7475-T73.  Schijve  also  reported  the 
results  of  Minakawa  and  McEvily  (46)  which  showed  K^p/Kmax  (equivalent  to  Pop/P™^^) 
increasing  with  decreasing  AK  (load  shedding).  This  trend  is  consistent  with  the  cP^p/Po^^x 
trend  in  Figure  5.19  after  the  formation  of  a  type  III  crack. 

The  bPop  closure  rat'o  (bPop/Pmax)  relatively  constant  at  a  level  of  0.21  for  a  type 
I  and  II  crack  in  Figure  5.19.  With  the  formation  of  a  type  III  crack,  bPop/Pmax  begins  to 
decrease,  reaching  a  value  of  about  0.15  at  the  test  maximum  crack  size. 

This  significant  reduction  in  closure  load  is  attributed  to  the  large  void  formation 
with  increasing  crack  size. 

Closure  Load  Variations  with  Load  History 

The  closure  load  data  presented  in  this  section  compares  the  constant  AK^,  constant 
load  and  block  loading  test  results  with  the  constant  AK^j  data  from  above.  The  closure  load 
absolute  values  (aPop,  bPop  and  cPop)  are  presented  for  the  three  tests  in  Figure  5.20 
through  5.22.  Comparisons  of  all  four  tests  are  made  in  Figure  5.23  through  5.28. 

Figures  5.20  through  5.22  appear  to  show  the  same  general  trend  in  closure  load  as 
seen  for  both  aPop  and  cPop  discussed  above  for  the  constant  AK^j  test  (Figure  5.18).  In 
Figure  5.21,  bPop  appears  to  show  a  different  trend  in  the  type  I  crack  size  region  since  it 
increases  before  decreasing  with  crack  growth.  Although  different,  the  trend  is  considered 
consistent  and  representative  of  the  crack  geometry  stabilizing  in  the  types  I  and  II  transition 
regions.  This  trend  of  closure  load  increasing  and  then  decreasing  was  also  observed  in 
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CLOSURE  LOADS  -  aPop,  bPop.  cPop  (Iba) 


CONSTANT  AKa  -  600 


Figure  5.20  Qosnre  load  (aPop,  bPop  ft  cPop)  as  function  of  crack  size 
(a/t)  from  the  constant  AlU  test 


117 


CLOSURE  LOADS  -  aPop.  bPop.  cPop  (lbs) 


CONSTANT  LOAD 


DEPTH  TO  THICKNESS  RATIO  (A/T) 


Figure  5.21  aosiue  load  (iPop,  bPop  St  cPop)  as  funcuon  of  crack  sire 
(aA)  from  the  constant  load  test 
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CLOSURE  LOADS  -  aPop,  bPop.  cPop  (lbs) 


Pigiifc  5.23  Comparison  of  the  aPop  dosure  load  for  all  four  tests  as  a 
fhnraon  of  crack  size  (aA) 
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Figure  S.26  Conqiarison  of  the  cPop/pMx  closure  ratio  as  a  function  of 
cradcsize 
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Ray's  (25)  data.  However,  when  Ray  repeated  the  same  test,  the  closvire  load  was  found  to 
continually  decrease  with  crack  size.  Also,  bP^p  measurements  made  in  this  investigation 
during  precracking  (presented  in  i^pendix  E)  ranged  from  146  to  176  lb.,  which  is 
equivalent  to  the  bP^p  levels  shown  in  Figure  S.21.  Once  a  type  II  crack  is  formed,  the  bPop 
trend  follows  the  same  decreasing  pattern  with  crack  growth  as  the  other  tests. 

The  block  loading  data  in  Figure  5.22  identifies  the  applied  load  transition  points  at 
the  top  of  the  plot  (Figure  4.12).  The  first  arrow  indicated  an  increase  in  applied  load  from  a 
delta  moment  of  644  in-lb.  to  773  in-lb.  The  next  two  arrows  represent  a  reduction  in 
applied  load  from  773  in-lb  to  644  in-lb  and  then  to  322  in-lb  for  the  remainder  of  the  test. 
These  load  transition  point  arrows  are  also  presented  on  ail  of  the  block  loading  plots  for 
reference. 

Absolute  closure  load  values  (aP^p,  bPop  and  cP^p)  for  all  four  tests  are  compared  in 
Figures  5.23  through  5.25  Each  of  the  Figures  present  characteristic  information  concerning 
the  closure  load  mechanism  as  discussed  in  items  A  through  D  below. 

A.  Overall  observations  from  Figures  5.23  through  5.25  are: 

1.  Closure  load  absolute  values  decrease  with  increasing  crack  size. 

2.  bPop  and  cPgp  load  levels  are  higher  for  the  constant  load  and  block  loading 

tests  than  for  the  constant  stress  intensity  factor  tests. 

3.  aPgp  trends  (Figure  5.25)  follow  the  same  slope  of  decreasing  closure  load 

with  increasing  crack  length  for  all  four  tests.  Therefore,  aPop  is  inversely 
proportional  to  crack  size. 


B.  Comparing  trends  from  the  constant  AK^  test  with  the  constant  AK^ 


tests: 

1.  bPop  and  cPop  are  higher  for  the  constant  AK^  tests  than  for  the  constant 

AK^  tests. 

2.  aPop,  bPop  and  cPgp  are  inversely  proportional  to  crack  size. 

C.  Comparing  the  constant  load  test  (continually  increasing  AK)  with  the 
constant  AK^  test: 

1.  The  constant  load  test  closure  loads  are  consistently  higher  than  the  constant 

AK^  levels.  This  is  anticipated  (as  discussed  in  the  approach)  since  the 
plastic  zone  size  is  continually  increasing  while  the  constant  AK^  plastic 
zone  remains  constant  after  precracking. 

2.  aPop  and  bP^p  closure  loads  follow  the  same  slope  once  a  type  II  crack  is 

formed. 


D.  Comparing  the  constant  load  test  with  the  block  loading  test: 

1.  cPgp  for  the  block  loading  test  are  only  slightly  higher  than  the  constant  load 

test  during  the  types  I  and  II  crack  formation.  Once  a  type  III  crack  is 
formed,  the  closure  loads  are  the  same. 

2.  cPgp  is  relatively  constant  for  the  constant  load  test  until  a  type  III  crack  is 

formed.  For  the  block  loading  test,  abrupt  increases  in  applied  load  cause 
cPop  to  increase,  and  decreases  in  applied  load  cause  cPop  to  decreased. 

3.  bPQp  is  only  slightly  larger  in  the  block  loading  test  than  in  the  constant  load 

test  once  a  type  II  crack  is  formed.  The  bPop  versus  a/t  slopes  are 
consistent  for  both  tests. 
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Comparisons  of  the  closure  load  ratios  (Pop^max)  shown  in 

Figures  S.26  through  5.28.  In  Figure  S.26,  three  of  the  four  tests  (constant  AKc>  load  and 
block  loading)  are  seen  to  have  average  closure  load  ratio  values  of  approximately  0.29  for 
0.0  «  a/t «  0.20  (during  transition  between  a  type  I  and  II  crack).  The  constant  AK^  test  is 
at  a  level  of  approximately  0.14. 

Closure  load  ratios  for  the  constant  AK<j  test  are  noted  to  increase  from  a  level  of 
about  0.29  before  the  formation  of  a  type  III  crack,  to  a  level  of  0.35  after  a  type  III  crack  is 
formed.  A  ratio  of  0.35  is  also  a  level  at  which  the  block  loading  test  stabilizes. 

Although  the  20%  load  change  at  a/t  =  0.07  and  0.11  during  the  block  loading  test 
caused  minor  changes  in  cPop/PiQax>  these  changes  were  noted  to  quickly  return  to  a  stable 
level  of  about  0.29  with  crack  growth.  However,  with  a  major  P^ax  reduction  of  50% 
at  a/t  =  0.25,  an  abrupt  change  in  cPop/Pj^^  occurs,  and  does  not  return  to  the  original  level 
for  the  reminder  of  the  test.  This  is  similar  to  results  observed  by  Ashbaugh  (47)  when 
testing  thru'thickness  compact  tension  specimens  of  Inconel  718  material.  He  observed  that 
with  a  60%  reduction  in  load  level  for  an  R  =  0.1,  the  crack  was  required  to  grow  twice  its 
length  in  order  to  return  to  the  original  closure  load  ratio. 

Although  the  constant  AKg  test  is  noted  to  vary  around  the  crack  type  transition 
points,  it  is  relatively  constant  for  the  complete  test.  The  constant  load  test  is  initially  at  a 
level  of  approximately  0.29,  and  reduces  to  a  level  of  0.22  with  crack  growth. 

Table  5.8  is  a  comparison  of  closure  loads  at  a  common  crack  size  (a  =  0.14,  c  = 
0.20).  Note  the  trend  of  closure  load  absolute  values  being  highest  for  the  block  loading  test 
and  lowest  for  the  constant  AK^  test.  The  closure  load  ratios  are  at  the  same  relative  level 
except  for  the  constant  AK^  test. 
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Table  5.8:  Closure  Load  Values  at  a  CommoD  Crack  Size  (a  =  0.14,  c  =  0.20) 


Test 

sEoo 

^£oo 

cPoo/ijTiax 

bPoD^ma; 

Block  Loading 

225 

152 

0.27 

0.19 

Constant  Load 

202 

139 

0.25 

0.17 

Constant 

115 

96 

0.25 

0.21 

Constant 

67 

53 

0.13 

0.09 
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VI.  EFFECTIVE  STRESS  INTENSITY  FACTOR  DEFINITIONS 


The  widely  accepted  definition  of  effective  stress  intensity  factor 
proposed  by  Elber  (1)  is  based  on  the  use  of  an  effective  load  range  (Pjjjax "  ^cl)  iitstead 
of  an  applied  load  range  (Pjnax  *  ^min)  predicting  fatigue  crack  growth  (refence 
Chapter  2,  "Plasticity-Induced  Closure").  Although  this  basic  de^nition  is  readily 
accepted,  deHning  the  closure  load  and  the  corresponding  stress  intensity  factor  (K^i)  is  a 
problem  which  has  plagued  the  fracture  mechanics  community  for  over  15  years.  For 
example,  numerous  papers  have  been  written  on  closure  load  measurement  problems 
ranging  from  instrumentation  type  to  measurement  location  (3,25,26).  Interpretation  of 
the  crack  internal  surfaces  opening  before  the  external  free-surface  for  thru-thickness 
cracks  (19,20,21)  is  another  problem  area.  These  considerations  are  all  relevant  for 
surface  flaws,  which  also  have  the  additional  complication  of  a  three-dimensional  crack 
boundary. 

An  approach  for  organizing  the  experimental  results  presented  in  Chapter  V  is 
essential  before  a  definition  of  AK^ff  can  be  proposed  for  surface  flaws.  That  rationale  is 
presented  in  this  chapter,  and  forms  the  basis  for  data  correlation  relating  to  various 
definitions  of  AK^ff  which  follow  in  the  next  two  chapters.  The  primary  variables 
associated  with  these  definitions  can  be  grouped  into  four  major  areas  of  influence: 

1.  The  closure  load  range  (Pmax "  ^cl’  ^cl  measured 

closure  loads  aPop,  bi  op  or  cPop)- 

2.  The  location  at  which  K  is  being  calculated  (i.e.,  crack  tip  location  A  or  C). 

3.  The  crack  opening  displacement  (COD)  associated  with  a  particular  crack  tip 
location  and  crack  size  (i.e.,  crack  types  1, 11  and  III). 

4.  The  stress  intensity  factor  solution  used  to  calculate  K. 


Organizing  the  above  four  areas  into  a  useful  format  requires  a  concise  set  of 
definitions.  The  first  two  areas  can  be  characterized  by  a  matrix  format  of  closure  load 
and  crack  tip  locations.  This  format  is  summarized  below,  and  will  be  used  extensively 
throughout  the  remainder  of  this  report.  The  closure  load  range  is  denoted  by 
superscripts,  the  crack  tip  locations  by  subscripts,  and  AK  and  AK^ff  are  defined  as 
follows: 

Applied  Stress  Intensity  Factor  (  KV 

AKg  =  applied  AK  at  location  A  due  to  (Pmax '  ^min) 

AK^  =  applied  AK  at  location  C  due  to  (Pmax  *  ^min) 

Effective  Stress  Intensity  Factor  (AKeff): 

AK^a  eff  “  ®^f^®ctive  AK  at  location  A  due  to  (Pjuax  ’  ^'^op) 

AK\  eff  “  ^  location  A  due  to  (Pmax  "  ^^op) 

AK^a  eff  =  effective  AK  at  location  A  due  to  (Pmax  *  cPop) 

AK^c  eff"  effective  AK  at  location  C  due  to  (Pmax  ‘  ^^^op) 

or 

AK‘j  gff  =  effective  AK  at  location  j  due  to  (Pj^ax  ‘  *^op) 

where  i  =  a,b,c  to  define  the  closure  load,  and  j  =  a  or  c  to  define  the  crack  tip  location  A 
or  C,  respectively. 

Organizing  the  third  area  from  above  (COD  influences)  requires  a  comprehensive 
evaluation  of  the  individual  crack  types,  and  their  relationship  to  crack  size,  closure  load 
and  crack  tip  location.  These  relationships  form  the  basis  for  deHning  AK^f^  and  address 
key  questions  such  as: 
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1.  How  do  the  three  crack  types  (1,  II  and  III)  influence  the  stress  intensity 

factor  around  the  crack  tip  boundary? 

2.  Which  of  the  closure  loads  (aPop,  bPop  or  cPop)  should  be  used  to  define 

Kcl? 

3.  How  are  the  crack  tip  locations,  crack  types  and  various  closure  load 

measurements  used  to  define 

The  fourth  area  for  consideration  (K-solutions)  will  be  discussed  in  each  of  the 
following  sections. 

Effective  Stress  Intensity  Factor  for  a  Type  I  Crack 

A  type  I  crack  opening  displacement  (COD)  pattern  was  described  in  Chapter  V 
("Typical  COD  Profiles")  as  being  classical  in  that  the  complete  crack  boundary  opens  at 
the  same  load  level.  That  is,  with  load  application  the  crack  initially  opens  at  the  crack 
free-surfece  (bPgp,  reference  Figure  5.15),  and  continues  to  open  inwardly  from  the  free- 
surface  until  the  internal  surfaces  and  crack  tip  boundary  are  completely  open  (such  that 
bPop  «  aPop  =  cPop»  reference.  Figures  5.2  and  5.3). 

Because  of  the  type  I  crack  opening  pattern,  classical  fracture  mechanics  K- 
solutioDS,  such  as  Newman-Raju  (5),  can  be  used  to  predict  AKgff  as  proposed  by  Elber 
(1,  eqn.  2.7).  That  is: 

^max  ■  ^cl 

AKgff  =  U  AK  =  -  AK  (6.1) 

^max  ■  ^min 

“  ^^ax  ■  ^1 
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Since  the  closure  load  (Pj,j)  needed  to  open  a  type  I  crack  tip  boundary  is  approximately 
equal  at  locations  A  and  C  (that  is,  aPQp  =  cP^p),  the  effective  load  range  is  calculated 
using  either  aPop  or  cPQp.  These  two  values  of  are  shown  in  Figure  6.1,  where 
the  ordinate  represents  either  or  K^,  and  the  abscissa  an  applied  load.  Then: 


^max  ■  ®^op 

-  («-2) 

P  -  P  • 

^max  *^min 

^max  ■  ^^op 

AK‘^ceff=  -  ^  (6-3) 


P  -  P  • 

*^max  *^min 

The  solid  line  (region  3)  defines  fully-open  COD  patterns  where  K-solutions  are 
available  for  calculating  AK^ff  A  double  dashed  line  extends  from  the  solid  line  to  the 
zero  load  to  identify  available  fully-open  K-solutions.  The  heavier  dashed  line  (extended 
horizontally  from  the  solid  line)  reflects  a  hypothesized  constant  K-level  associated  with 
the  crack  tip  boundary  being  completely  closed.  This  conclusion  is  based  on  the 
observation  that  at  cPgp  and  lower  loads  the  crack  tip  is  closed,  and  represents  two 
possible  solutions  either  the  K-level  remains  constant,  or  the  K-level  is  instantaneously 
zero  as  load  is  reduced  below  cPQp.  To  choose  the  later  would  require  K  to  be  a 
discontinuous  function,  and  therefore  the  choice  of  a  constant  K-level  was  made. 

For  loads  between  bP^p  and  cPop  (region  2),  the  crack  tip  is  closed  while  the 
crack  internal  surfaces  are  open.  For  loads  between  zero  and  bPop  (region  1),  both 
surfaces  are  closed.  For  both  region  (1)  and  (2),  the  K-level  is  considered  to  be  constant 
since  the  crack  tip  is  closed. 

Load  regions  (1),  (2),  and  (3)  are  identified  on  this  plot  to  illustrate  their  expected 
relationship  to  the  classical  definition  of  closure  on  a  load  versus  displacement  plot. 
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STRESS  INTENSITY  FACTOR 


Kc  or  Ka 


where  displacement  is  measured  at  the  crack  middle  firee-sur£ace  (Figure  6.2).  This  load- 
displacement  format  is  similar  to  data  presented  in  Figure  5.13. 

For  the  R  ratio  of  0.035  used  in  investigation,  bPop>  is  always  greater  than  PmiQ. 
Therefore,  for  an  applied  load  less  than  cP^p  the  crack  tip  boundary  is  always  closed.  If 
a  higher  R  ratio  were  used,  such  that  Pmin  were  ^eater  than  cPop,  P^.]  would  be  defined 
as  PgjjQ  instead  of  cP^p. 


Effective  Stress  Intensity  Factor  for  a  Type  II  Cracks 

A  type  II  crack  is  distinguished  by  the  formation  of  a  "void"  internal  to  the  crack 
surfaces  while  the  crack  tip  boundary  is  closed  under  zero  load  (reference  Chapter  V, 
"Typical  COD  Profiles,"  Figures  5.2  and  5.3).  This  "void"  formation  complicates  the 
definition  of  K  around  the  crack  tip  boundary  because  the  load  required  to  open  the  crack 
tip  at  location  A  (aPop)  is  different  from  the  load  required  to  open  location  C  (cPop). 
Therefore,  for  applied  loads  between  aPop  and  cPop,  the  crack  tip  is  open  at  location  A 
while  it  is  closed  at  location  C.  This  forms  a  crack  tip  boundary  which  is  "partially- 
open"  and  "partially-closed"  (as  opposed  to  "folly-open"),  for  which  there  are  no  K- 
solutions  known  to  be  available. 

The  crack  opening  displacement  (COD)  process  for  a  type  II  crack  is  shown 
schematically  in  Figure  6.3,  where  is  plotted  as  a  function  of  applied  load.  The 
opening  process  is  separated  into  four  regions  such  that  region  (1)  is  between  a  load  of 
zero  and  aPop,  region  (2)  between  aPop  and  bPop,  region  (3)  between  bP^p  and  cPop, 
and  region  (4)  between  cPop  and  Pm^x*  The  hypothesized  partially-open  K-level 
experienced  by  the  crack  tip  in  the  closure  load  region  (1-3)  is  identified  on  the  plot  by 
the  heavy  dashed  line.  The  solid  line  (region  4)  represents  available  "folly-open"  K- 


135 


LOAD 


TYPE  i  CRACK  SCHEMATIC 


Figure  6.2  Schematic  of  load  veisus  displacement  at  the  crack 

free-surface  to  relate  closure  load  regions  (1^,3)  to  figure 
6.1  for  a  type  I  crack 


PmM  aPop  bPop  cPop 


f 


APPUED  LOAD 


Fignie6J  Schematic rqaejenwiop of K is » fonction of i 
for  a  type  n  crack  « locaiioa  A  (4ii^  degrees) 
potenul  AKt  eff  definitions  are  represented 


solutions  when  the  crack  is  completely  open  for  loads  greater  than  cPop.  Each  of  these 
COD  regions  is  described  in  the  following  paragraphs. 

In  region  1  of  Figure  6.3,  the  crack  tip  boundary  is  closed  at  zero  load.  Upon  load 
application,  the  internal  void  increases  in  size  and  initially  opens  the  crack  tip  at  location 
A  at  a  load  level  of  aPop.  Since  the  crack  tip  is  closed  below  aPop,  there  is  little  change 
in  crack  tip  geometry  and  remains  relatively  constant.  The  projected  Kg  levels  are 
depicted  by  a  heavy  dashed  line  since  K-solutions  are  not  directly  available  for  this  crack 
tip  geometry. 

As  load  is  increased  above  aPQp,  the  void  increases  in  size  until  the  crack  firee- 
surface  (location  B)  at  the  bP^p  load  level.  Below  bP^p  the  crack  free>sur£ace  is  closed 
while  the  crack  tip  boundary  is  open  between  locations  A  and  C,  and  results  in  the 
specimen  geometry  being  relatively  stiff.  Therefore,  changes  in  Kg  with  increased  load 
are  expected  to  be  very  small,  as  well  as  being  nonlinear  because  of  the  continual  change 
in  crack  geometry  between  apQp  and  bpQp  (region  2).  As  load  is  increased  above  bPop, 
changes  in  Kg  are  expected  to  be  much  greater  with  aj^lied  load  since  the  geometry  is 
more  compliant  due  to  the  crack  firee-surface  being  open.  Between  bPop  and  cP^p 
(region  3),  the  crack  surfoces  and  tip  boundary  continues  to  open  until  the  cPop  load  is 
reached  and  the  crack  is  fiilly-open.  Between  cPop  and  Pq^^x  (region  4),  Kg  (solid  line) 
is  linear  with  applied  load  since  the  crack  is  hilly-open  and  the  geometry  remains 
constant  That  is,  the  slope  of  the  linear  portion  of  the  Kg  versus  applied  load  curve  is 
proportional  to  the  geometry  factor  corresponding  to  a  fiilly-open  surface  flaw,  and  Kg 
can  be  dehned  using  available  K-solutions  such  as  K  =  P*F(a,a/c,^);  where  F(a,a/c,^)  is 
the  boundary  correction  factors  conesponding  to  a  fully-open  surface  flaw  of  depth  "a", 
aspect  ratio  "a/c"  and  location  angle 

In  summary,  changes  in  actual  Kg  are  nonlinear  with  respect  to  applied  load 
between  zero  and  cP^p  because  of  the  continually  changing  crack  geometry  with 
increasing  load.  Kg  is  hypothesized  to  increase  from  a  constant  value  between  a  load  of 
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zero  and  aP^p  (region  1),  and  to  continue  along  a  nonlinear  path  (regions  2  and  3)  until  it 
reaches  the  cPop  load.  It  then  becomes  linear  (region  4)  up  to  and  including  the 
maximum  applied  load. 

The  relationship  of  actual  experienced  by  the  crack  tip  at  location  A 
(discussed  above)  and  calculated  AK^  eff  ^  Figure  6.3.  Here,  three  separate 

eff  presented  by  using  available  fully-open  K-soiutions  and  measured 

closure  loads  (aPop,  bPop  and  cPop),  such  that: 

Pmax 

eff  - - ^  C*-*) 

P  -P  • 

^max  ^min 

.  ^max  ■ 

AK**,  eff - - AK,  (6  J) 

P  -  P  ■ 

max  mm 

Pmax-‘=P0P 

eff . AKa  (6.6) 

P  -  P  • 

‘^max  *^min 

where  AKj^  is  deflned  from  the  fully-open  K-solution  (solid  line),  and  the  double-dashed 
line  extended  to  the  ordinate  and  abscissa  (for  Kjjj =  0).  Here,  AK® ^  and  AK^^  gff  is 

projected  to  bound  the  actual  AK3  value  (shown  by  heavy  dashed  line),  and  it  would 
appear  that  the  best  P^j  value  to  define  the  actual  AK^  is  between  aPop  and  bPop  when 
using  available  K-solutions.  On  the  other  hand,  not  having  accurate  K-soiutions  to 
represent  the  actual  crack  geometry  leaves  this  selection  open  for  further  evaluation 
(Chapter  VII).  Note  that  all  three  of  the  AK'^  gff  definitions  provide  values  which  are 
less  than  the  applied  AK^  valued  defined  by  -  Kjjjjjj. 

Extending  this  evaluation  to  the  classical  definition  of  closure  load  on  a  load- 
displacement  plot  can  be  seen  schematically  in  Figure  6.4.  In  this  format,  selection  of  the 
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Figure  6.4 


Schematic  of  load  versu-s  displacement  at  the  crack 
free-surface  to  relate  closure  load  regions  (1,2^, 4)  to 
Figure  6.3  for  a  type  II  crack 


aPop  to  bPop  closure  load  range  to  calculate  actual  also  corresponds  to  a  primary 
definition  of  closure  load  for  thru-thickness  cracks.  It  should  be  noted  that  although  the 
slope  changes  in  region  (2)  and  (3)  may  be  minor  and  difficult  to  detect  by  conventional 
instrumentation  (clip  gauges,  back  face  strain,  etc),  they  represent  definite  changes  in 
crack  geometry  which  were  described  in  Chapter  V  (Figure  5.13). 

In  comparison,  AK^,  at  location  C  (Figure  6.5)  is  much  simpler  to  evaluate 
than  Here,  the  crack  is  closed  at  location  C  until  the  cPop  load  level  is  reached 

(defined  as  the  last  portion  of  the  crack  surfaces  to  open).  Thereafter,  the  crack  is  fully- 
open  and  is  calculated  from  available  formulae  using  an  effective  load  range  of 

Pmax  ‘  ^^op-  The  COD  representation  on  a  load  versus  displacement  plot  is  shown  in 
Figure  6.4  as  region  (4).  This  analysis  is  valid  for  all  three  crack  types. 


Effective  Stress  Intensity  Factor  for  a  Type  III  Crack 

The  difference  between  a  type  II  and  type  III  crack  is  associated  with  the  "void" 
displacement  having  reached  a  sufficient  size  to  geometrically  open  the  crack  tip 
boundary  at  location  A  at  :^ro  load.  That  is,  the  crack  tip  at  location  A  is  open  at  zero 
load  (aPop  s  0),  while  other  locations  around  the  crack  tip  boundary  are  closed 
(reference  Figures  5.2  and  5.3). 

The  process  of  defining  and  AK^  gff  for  a  type  III  crack  at  location  A  is 
shown  schematically  in  Figure  6.6  where  is  plotted  as  a  function  of  applied  load.  It 
should  be  noted  that  for  a  type  III  crack,  is  not  zero  at  zero  load  since  the  crack  tip  is 
open,  but  is  a  variable  which  is  directly  proportional  to  crack  size.  That  is,  with  increased 
crack  growth  the  void  size  increases  (area  and  displacement),  causing  a  larger 
displaceihent  at  location  A  and  more  of  the  crack  tip  boundary  between  location  A  and  C 
to  open.  The  crack  tip  boundary  length  which  is  open  at  zero  load  has  been  measured  to 
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STRESS  INTENSITY  FACTOR  AT  t 


STRESS  INTENSITY  FACTOR  AT  A' 


Figiire6.6  Sdwmttic  representation  of  Kcs  a  function  of  m^ied  load 
for  a  type  in  crack  at  location  A  degrees).  Three 

potennal  AKa  cff  definitkxis  are  iqwesented 


range  from  17%  for  a/t  =  0.208  to  65%  for  a/t  =  0.260  (reference  Chapter  V,  "Typical 
COD  Profiles"). 

For  load  levels  between  zero  and  bPop,  the  crack  free-surface  is  closed,  the 
resultant  specimen  geometry  is  relatively  stiff,  and  changes  to  with  increased  load  are 
expected  to  be  small  and  nonlinear.  As  load  is  increased  above  bPop,  the  crack  geometry 
is  more  compliant  since  the  crack  free-surface  is  open,  and  the  rate  of  change  with 
applied  load  is  expected  to  be  much  greater.  This  process  is  identified  as  a  nonlinear 
heavy  dashed  line  between  aP^p  and  cP^p  in  Figure  6.6  since  K-solutions  for  a  type  III 
partially-open  crack  (with  an  internal  void  formation)  are  not  known  to  exist  --  thus 
precluding  from  be  determined  directly. 

For  loads  above  cPQp,  increases  linearly  with  load  since  the  ;.rack  is  fully- 
open  and  the  geometry  remains  constant.  These  values  can  be  defined  using  available 
K-solutions  such  as  K  =  P*F(a,a/c,<^),  and  is  represented  by  a  linear  solid  line  between 
cPop  and  the  maximum  applied  load.  The  correlation  of  closure  load  regions  on  a  load 
versus  displacement  plot  are  projected  to  be  represented  by  regions  2  through  4  on  Figure 
6.4.  However,  since  compliance  is  greater  for  large  cracks  (reference  Figure  5.14), 
regions  (2)  and  (3)  will  not  be  as  distinct  as  with  a  type  II  crack.  This  observation  is 
consistent  with  measurements  by  Bucci  (48),  who  presented  load-displacement  data  for 
thru-thickness  cracks  which  showed  a  smoothing  of  the  nonlinear  regions  2  and  3,  as  well 
as  a  reduction  in  closure  load  with  increased  crack  sizes  for  a  high  strength  aluminum 
alloy. 

As  can  be  seen  in  Figure  6.6,  each  of  the  three  closure  loads  (aP^p,  bPop  and 
cPop)  can  be  used  to  define  three  distinct  effective  stress  intensity  factors  at  location  A 
(AK^^  gffi  AK^g  gffi  AK^g  gff).  Again,  only  estimates  of  actual  ABCg  can  be  determined 
using  available  fully-  open  K-solutions  and  measured  load  ranges.  From  the  Figure  6.6 
schematic  it  would  appear  that  the  actual  AKg  is  best  defined  between  the  aPop  and  bPop 
closure  loads.  As  with  a  type  II  crack,  this  evaluation  is  difficult  to  substantiate  using 
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fully-open  K-solutions  when  the  crack  is  partially-open  (part  of  the  crack  tip  is  closed 
while  other  parts  are  closed). 

Note  that  and  values  are  less  than  the  classical  applied  AK^ 

defined  by  K^-K„i„.AK>a  gff  is  either  greater  than  applied  AK^  for  Pnjin  loads 
greater  than  zero,  or  equal  to  AK^  when  =  0. 

Again,  the  definition  of  AK’^  for  a  type  III  crack  is  much  simpler  than  AK*^ 
gff  as  depicted  in  Figure  6.5.  Here,  the  crack  is  closed  at  location  C  until  the  cPop  load 
level  is  reached  and  the  crack  becomes  fully-open.  Then,  K^.  is  calculated  using  available 
formulae,  and  AK^^^  ^ff  is  a  function  of  the  effective  load  range  (Pjjjax '  ‘^^op)- 


Closure  Load  Considerations 

It  is  apparent  from  the  previous  discussion  that  the  void  formation  (reference 
Figure  5.4)  is  not  modeled  by  available  K-solutions.  Although  it  is  not  apparent  this 
represents  a  calculation  problem  at  the  crack  free-surface  (reference  Figure  5.1,  K^),  it 
does  present  a  calculation  problem  at  location  A  (K^,  <(»  =  90  degrees). 

For  a  type  I  crack,  the  aPop  or  cPgp  closure  loads  (Pj.j)  can  be  used  to  define 
AK^ff  using  an  effective  load  range  (Pmax  *  ^cl)  fiiHy-open  K-solutions  at  either 
location  A  or  C  of  the  crack  tip  boundary.  For  a  type  II  crack,  it  is  not  known  whether 
aPop  or  bPop  best  represents  the  location  A  effective  load  range.  K-solutions  are  not 
known  to  exist  for  this  crack  geometry.  For  a  type  III  crack  the  closure  load  at  location 
A  is  zero,  and  the  crack  residual  opening  is  a  function  of  crack  sizes  (larger  cracks  have 
larger  internal  displacement,  reference  Figure  5.8).  Again,  K-  solutions  are  not  known  to 
exist  for  this  partially-open  crack  geometry.  These  considerations  for  selecting  a  closure 
load,  and  the  need  for  K-solutions  which  match  the  crack  geometry,  are  summarized  in 
Tables  6.1  and  6.2. 
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Table  6.1:  Hypothesized  Closure  Load  and  K-Solution  Needs  to  Define 
at  Crack  Tip  Location  A 


Crack  Tvoe 

Closure  Load 

fP.,!!  Variable 

K-solution 

I 

cPop 

OK 

II 

bPop 

Needed 

III 

bPop 

Needed 

Table  6.2:  Hypothesized  Closure  Load  and  K-Solution  Needs  to  Define 
at  Crack  Tip  Location  C 


Crack  Tvoe 

Closure  Load 

(P^il  Variable 

K-solution 

I 

cPop 

OK 

II 

cPop 

OK 

IH 

cPop 

OK 

It  is  apparent  from  Table  6.1  that  although  a  need  exists  for  partially-open  crack 
K-solutions,  available  fully-open  crack  models  represents  a  major  portion  of  the  crack 
growth  process  and  are  expected  to  provide  good  estimates  of  partially-open  cracks. 
First,  AKg  (location  C)  can  be  obtained  from  the  effective  load  range  cPop)» 
and  existing  K-solution  for  a  fiilly-open  crack.  Second,  for  both  type  II  and  type  III 
craclm,  if  is  hypoth**  zed  that  AK^  (location  A)  can  be  approximated  from  the  (Pna^x 
-  bPop)  load  range  and  available  K-solutions  for  partially-open  cracks.  Since  the  region 
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between  bPQp  and  cPqp  corresponds  to  a  crack  geometry  such  as  that  shown  in  Figures 
5.4c>d,  Kg,  although  not  known,  is  felt  to  be  representative  of  a  fiilly-open  crack.  Thus, 
the  slope  between  bP^p  and  cPop  in  Figures  6.3  and  6.6  should  be  similar  to  that  of  a 
fully-open  crack;  allowing  the  use  of  available  K-solutions  and  the  (P^ax '  ^^op) 
effective  load  range. 

For  a  type  II  crack,  the  K-solution  at  location  A  (reference  Figure  5.4b  crack 
geometry),  although  not  known,  is  expected  to  yield  smaller  values  than  for  a  fully-open 
crack  because  of  the  inherent  "stiffness"  of  that  geometry  due  to  closure  along  the  free- 
surface.  That  is,  the  slope  of  the  curve  in  Figure  6.3  between  aPop  and  bPop  is  expected 
to  be  smaller  than  for  a  fully-open  crack.  Selecting  the  (Pmax  *  ^^op)  range  for 
determining  AKg  would  appear  to  overestimate  the  crack  tip  Kg  value.  Thus,  the  use 
of  bPop  for  calculating  AKg  appears  to  be  a  reasonable  average  choice  instead  of 
aPop  which  becomes  zero  when  transition  occurs  between  a  type  II  and  type  III  crack 
(reference  Figure  5.18). 

Therefore,  the  definition  of  AK^^f  depends  on  the  crack  tip  boundary  location  (4»), 
crack  depth,  aspect  ratio  and  selection  of  a  closure  load  which  corresponds  to  a  given 
crack  type.  Considering  the  current  status  of  K-solutions,  selection  of  a  representative 
closure  load  (Pj,j)  to  define  an  effective  load  range  (P^ax  ‘  ^cl) 
hypothesized  at  location  A.  For  these  cases,  an  average  load  range  of  (P^ax  ’  ^^op)  ** 
considered  a  reasonable  choice,  such  that: 

eff  =  (Pmax  '  ‘>'’op)  F(m/c,4.=90)  (6.7) 

It  is  of  interest  to  note  that  the  degree  of  difficulty  associated  with  the  location  A 
crack  tip  has  also  been  reported  by  other  investigators.  A  detailed  review  of  surface 
crack  K-solutions  and  aspect  ratios  for  plates  in  bending  by  Mahmoud  and  Hosseini  (41), 
identified  major  prediction  errors  at  the  location  A  crack  tip.  Their  comparison  of  the 
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Newman-Raju  K-solutions  to  photoelastic  experimental  results  identified  a  56%  error  for 
one  case  at  location  A,  versus  a  ±14%  error  at  the  crack  free-surface.  The  authors 
attributed  this  enor  to  experimental  extrapolation  from  a  two-dimensional  specimen  to  a 
three-dimensional  crack.  Discrepancies  for  surface  flaws  were  also  identified  by  Carter, 
et  al.  (42)  when  testing  aluminum  in  tension.  This  study  reported  that  "It  has  not  been 
determined  whether  the  source  of  error  is  due  to  fundamental  different  crack  growth  rates 
at  surface  and  depth  or  due  to  inaccurate  stress  intensity  solutions.” 


VII.  EFFECTIVE  STRESS  INTENSITY  FACTOR 


CORRELATION  WITH  FATIGUE  CRACK  GROWTH  RATES 

Identification  of  the  three  crack  types,  and  their  influence  on  crack  opening 
displacement  (COD),  closure  loads  and  effective  stress  intensity  factor,  have  raised  a  number 
of  questions  concerning  the  prediction  of  fatigue  crack  growth  (FCG).  One  of  the  more 
significant  questions  is  based  on  selecting  the  load  range  which  "best"  defines  effective  stress 
intensity  factor  using  fiiliy-open  K-solutions.  The  approach  for  selecting  this  "best"  effective 
load  range  is  formulated  around  correlating  experimental  FCG  rate  data  —  as  described  in 
this  chapter.  The  evaluation  is  conducted  using  data  from  Chapter  V,  along  with  the 
effective  stress  intensity  factors  definitions  discussed  in  Chapter  VI. 

It  merits  re-emphasizing  that  this  chapter  will  utilize  available  K-solutions  to  bound 
the  crack  growth  rate  prediction  by  evaluating  various  effective  load  ranges.  In  effect,  not 
haveing  partially-open  K-solutions  for  a  "void"  geometry  (as  experimentally  identified  in  this 
investigation)  precludes  selecting  the  optimum  effective  load  range.  Therefore,  an 
evaluation  of ‘.y-e  "best"  correlation  of  the  data  utilizing  linear  elastic  fully-open  K-solutions 
is  presented.  This  condition  only  effects  location  A,  not  location  C  where  the  crack  is  fully- 
open  above  the  cPgp  closure  load.  The  K-soiution  differences  at  location  A  are  expected  to 
be  small  as  discussed  in  the  previous  chapter. 


Closure  Load  Stress  Intensity  Factor 

Selecting  closure  load  (Pmax  ‘  ^cl)  which  produces  the  least  data  scatter  can  be 
accomplished  by  evaluating  the  ability  of  several  AKgf^  definitions  to  correlate  fatigue  crack 
growth  rate  (FCGR)  data  into  the  specimen  depth  (da/dN)  ar  g  the  crack  free-surface 
(dc/dN).  That  is,  the  evaluation  process  for  selecting  the  mos.  ^  resentative  closure  load 


range  to  match  the  experimental  data  is  based  on  the  correlation  of  predicted  da/dN  and 
dc/dN  versus  applied  AK  and  several  definitions  of  The  matrix  of  AK^ff  definitions 

are  based  on  the  bPop  and  cpQp  closure  load  effects  at  locations  A  and  C  (Table  7.1).  These 
four  combinations  of  effective  load  ranges  are  represented  by  AK*^^;  gff,  AK'^j.  gff 

and  AK^j  gfj  and  a  cross  term  of  AK^^  AK^^g  gff  From  the  evaluation  in  Chapter 

VI,  the  logical  selections  for  correlating  da/dN  data  is  AK^^  g^^  and  for  dc/dN  is  AK^^g  gff 
(assuming  partially-open  K-solutions  are  available). 


Table  7.1:  Closure  Stress  Intensity  Factor  (Kgj)  Variables  Used  to  Evaluate  Several  AKg^ 
definitions 


FCGR 

cPop  Closure 

bPop  Closure 

Direction 

Load  Effects 

Load  Effects 

dc/dN 

da/dN 

An  example  calculation  of  .c  using  Kgj  at  locations  C  with  the  cPpp  closure  load 

follows  from  equations  2.7  and  2.c  a  that: 

AKe 

where:  Ug  =  (Pmax  *  P^cl)  ^  (^max '  ^min) 

=  (1.  -  cPop/Piuax)  ^  ~  (^*1) 
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Then: 


eff  “  max) 

“  ^  max  ■  cl  (^•^) 

Likewise,  at  location  A: 

^^a  eff  ~  ^  max "  ^^a  cl  (^•^) 

The  Kj,j  values  defined  in  Table  7.1  were  calculated  at  locations  A  (K3  ^j)  and  C  (K^. 
(,|)  using  measured  closure  loads  as  a  function  of  crack  size.  The  data  were  then  curve  fitted 
with  a  second  order  polynomial  as  shown  by  the  dotted  lines  in  Figures  7.1  through  7.4. 
Figures  7.1  and  7.2  present  the  stress  intensity  factors  at  location  C  due  to  cPqp  and  bP^p 
and  gj),  respectively.  Figures  7.3  and  7.4  present  and  j.l> 
respectively.  The  K^j  format  was  selected  over  percent  closure  load  (Pop/Pmax)  * 

more  continuous  function  -  especially  for  tiie  block  loading  data  (Figure  5.26  and  5.27). 

The  coefficients  for  these  curves  (k,l,m)  are  given  in  Table  7.2  through  7.5,  where  the 
general  equation  of  K^|  is  presented  as  a  function  of  a/t  such  that: 

K‘j  cl  =  k  +  P(a/t)  +  m*(a/t)2  (7.4) 

where  i  =  b  or  c  for  the  closure  load  bPop  or  cPgp,  respectively,  and  j  =  a  or  c  for  the  crack 
tip  locations  A  and  C,  respectively. 
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Figure  7. 1  S»re«  intensity  factors  at  the  crack  tip  free-suiface 

degrees)  due  to  the  cP<^  closure  load  level  as  a  function  of 

cradr  size  (aA)  for  all  tests 


(E/TJmUTMTSd)  SOVOl  dOdP  iV  ax 
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Table  7.2:  Polynomial  Coefficients  for  cl  “  ®  Function  of  a/t  (Figure  7.1) 


Test 

Ic 

i 

m 

Constant  AK^. 

144 

223 

-124 

Constant  AK^ 

89 

140 

-337 

Constant  Load 

116 

1227 

-2008 

Block  Loading 

183 

740 

-1059 

Table  7.3:  Polynomial  Coefficients  for  cl  “  ^  Function  of  a/t  (Figure  7.2) 


Test 

k 

1 

m 

Constant  AK^ 

115 

188 

-508 

Constant  AK^ 

95 

-51 

-381 

Constant  Load 

-3 

1953 

-4684 

Block  Loading 

123 

583 

-1010 

Table  7.4;  Polynomial  Coefficients  for  Function  of  a/t  (Figure  7.3) 


Test 

k 

1 

SI 

Constant  AK^. 

145 

159 

-371 

Constant  AK^ 

89 

18 

-128 

Constant  Load 

94 

1351 

-2657 

Block  Loading 

216 

369 

-1018 
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Table  7.5:  Polynomial  Coefficients  for  cl  ^  ^  Function  of  a/t  (Figure  7.4) 


Test 

k 

1 

Constant  AK^. 

118 

83 

-526 

Constant  AK^ 

86 

• 

00 

-262 

Constant  Load 

-6 

1884 

-4676 

Block  Loading 

154 

243 

-778 

EF^ctive  Stress  Intensity  Factor  Correlation 

The  fatigue  crack  growth  rate  (FCGR)  data  are  plotted  in  figures  7.5  through  7.7  as  a 
function  of  applied  AK  for  the  four  tests.  The  FCGR  data  (da/dN  and  dc/dN)  are  calculated 
using  a  modified  incremental  polynomial  smoothing  routine  developed  by  Larsen  (49).  The 
differences  in  data  scatter  between  locations  A  and  C  are  shown  in  Figures  7.5  and  7.6  for 
da/dN  and  dc/dN,  respectively,  while  Figure  7.7  combines  them  on  one  plot.  The  solid  line 
fitted  through  the  surface  flaw  data  has  Paris  Law  constants  of  C  =  5.737  E-22,  and  n  = 
5.5945  (reference  eqn.  2.1).  The  dashed  line  is  data  generated  by  Perez  (37)  using  thru¬ 
thickness  compact  tension  (CT)  specimens  from  the  same  sheet  of  PMMA  as  was  used  in  this 
experiment.  The  Paris  Law  constants  for  this  reference  data  is  C  =  4.807  E-20,  and  n  = 
5.0806.  This  reference  line  shows  the  surface  flaw  crack  growth  i  Ue  to  be  slower  than  the 
thru-thickness  FCGR.  It  is  important  to  note  that  both  da/dN  and  dc/dN  for  the  constand 
AKg  test  are  greater  than  the  AK<j  test.  These  differences  reflect  the  influence  of  plastic 
zones  on  FCGR. 

The  AKgff  versus  da/dN  and  dc/dN  data  are  shown  in  Figures  7.8  through  7.13  for 
the  four  closure  load  ranges  --  represented  by  AK‘^^,  ^.j,  AK^j,  ^.j,  AK'^^  and  AK^g  ^.j. 
Closure  load  effects  on  both  da/dN  and  dc/dN  are  shown  separately  in  Figures  7.8  through 
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da/dN  CRACK  GROWTH  (in/cyl) 


Figiiie7.5  Cnc±Bo^nie«locatioa A(daAiN)diietothefiofiuiMl 
qipUea  load  (Kaax'Kna)  for  all  four  leit  cy^es 


158 


dc/dN  CRACK  GROWTH  (In/cyl) 


APPLIED  DELTA  Kc  (p8iMin)>«l/2) 


Pigiiie7.6  Cncknowtfa  rale  at  looidonC  (dc/dN)  <ine  to  the  nominal 
appUed  load  (KMK-KaiH>fbr  all  four  teat  cycles 
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da/dN  s  dc/dN  CRACK  GROWTH  (In/cyl) 


APPLIED  DELTA  K  (PSIiiXN»ai/2) 


Figufe  7.7  Crack  growth  rate  at  iocadons  A  (daAIN)  and  C  (dc/dl^O  dne 
to  the  nominal  load  (Kaax-Ki^  for  all  four  test  cycles 
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da/dN  USING  bPop  (in/cyl) 


EFFECTIVE  oKb  (p8HHn*ai/2) 


Figme  7.8  Crack  me  at  location  A  (da/dN)  dne  to  the  effective 

load  ad)  for  «11  four  test  cycles 
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da/dN  USING  cPop  (In/cyl) 


dc/dN  USING' cPop  (in/cyl) 


EFFECTIVE  AKc  (psiliiniiiil/2) 


Ptgaie7.ll  Crack  gi^h  nte  at  location  C  (dc/dN)  due  to  the  effective 
load  c  d)  for  all  four  test  cyclc» 
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da/dN  fi  dc/dN  USING  bPop  (in/cyl) 


EFFECTIVE  DELTA  K  (PSlMlN*al/2) 


Figuie  7.12  Crack  grow^  rate  at  locwon  A  (da/dN)  and  C  (dc/dN)  due 
to  the  effective  load  (AK^d)  for  ail  four  test  cycles 
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da/dN  fi  dc/dN  USING  cPop  (In/cyl) 


EFFECTIVE  DELTA  K  (PSlMlNi«H/2) 


Figure  7.13  Crack  grow^  rate  at  locttion  A  (da/dN)  and  C  Cdc/dN)  due 
to  the  effective  load  (AK  cl)  for  all  four  test  cycles 
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7.11,  and  then  combined  in  Figures  7.12  and  7.13.  Thru-thickness  CT  specimen  results  are 
also  plotted  on  these  figures  for  reference  (dashed  line)  —  the  reference  line  is  for  applied  Ak 
and  is,  therefore,  not  corrected  for  closure. 

First,  consider  the  effect  of  closure  load  on  da/dN.  The  bPop  closure  load  (AK^^  cl» 
Figure  7.8)  effects  are  noted  to  make  only  minor  improvements  when  compared  with  applied 
AK  data  in  Figure  7.5.  The  effects  of  cPgp  (AK^^  j.j.  Figure  7.9)  are  found  to  be  slightly 
worse  than  applied  AK  (Figure  7.5)  --  particularly  at  the  lower  K  levels  where  the  block 
loading  data  is  separated  from  the  other  data  points. 

In  contrast  to  da/dN,  closure  load  effects  along  the  crack  free-surface  (dc/dN,  Figures 
7.10  and  7.11)  are  noted  to  show  excellent  correlation.  The  bPop  closure  load  (AK^j.  j.], 
Figure  7.10)  coalesces  both  the  constant  AK  and  the  constant  load  test  results  when 
compared  with  applied  AK  data  in  Figure  7.6.  The  effects  of  cP^p  (AK*^^.  j.j,  Figure  7.11) 
are  also  found  to  be  very  good,  but  to  produce  slightly  more  scatter  than  bP^p. 

Of  the  four  AKg^f  definitions,  the  best  correlation  (least  scatter)  is  that  of  the  dc/dN 
data  seen  in  Figure  7.10,  where  the  use  of  K^^.  (bPop  closure  load,  or  the  load  required  to 
open  the  middle  free-surface)  collapses  even  the  block  loading  data.  The  correlation  of 
dc/dN  in  Figure  7.11  using  cPop  is  also  very  good.  The  best  da/dN  correlation  is  seen  in 
Figure  7.8  (K*^^  ^.j),  where  the  bPop  closure  load  is  again  found  to  best  correlate  the 
experimental  data. 

Figures  7.12  and  7.13  combines  the  da/dN  and  dc/dN  FCGR  data  using  bPop  and 
cPop  (K^j  cl  *ttd  K^j  cl),  respectively.  From  this  evaluation,  use  of  the  bPop  is  found  to  best 
correlate  both  the  da/dN  and  dc/dn  FCGR  data  (Figure  7.12).  This  correlation  is  noted  to  be 
quite  good  except  for  the  da/dN  block  loading  data. 

It  should  be  noted  that  correlation  of  the  bPop  closure  load  is  consistent  with  thru¬ 
thickness  FCGR  investigations  which  measure  free-surface  crack  displacement  using  clip 
gauges  ahd/or  interferometry  methods. 
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Fatigue  Crack  Growth  Predictions 


The  above  correlation  of  closure  loads  for  collapsing  FCGR  data  can  also  be 
evaluated  by  using  analytical  models  to  compare  experimental  measurements  with 
predictions  of  cycles  and  crack  aspect  ratio  for  a  given  crack  length.  Although  the  approach 
is  not  as  direct  as  was  presented  in  the  previous  section  because  of  computational  influences, 
sample  calculations  are  presented  in  this  section  for  reference.  These  comparisons  are  made 
using  four  different  Paris  Law  constants  (C  and  n)  —  two  cases  using  applied  AK 
("uncorrected"  constants),  and  three  cases  using  "corrected"  constants  derived  from 
definitions  discribed  in  the  previous  section.  Since  the  closure  load  corrected  constants  are 
derived  from  this  investigation  and  are  not  independent,  the  calculations  are  presented  as  an 
example  of  how  this  approach  is  applied. 

The  analytical  program  developed  for  this  effort  calculates  crack  stress  intensity 
factors  in  the  "a"  and  ”c"  directions  (K^  and  K^)  using  the  Newman-Raju  boundary 
correction  factors  at  s  90  and  0  degrees,  respectively.  As  discussed  in  Chapter  VI, 
although  K-solutions  for  AK^  in  the  closure  load  region  are  not  known  to  be  available  for 
types  II  and  III  cracks,  the  use  of  effective  load  ranges  and  existing  K-solutions  are  expected 
to  produce  reasonable  results.  The  model  uses  either  the  applied  stress  intensity  factor  range 
of  (AK  =  Kqj^  -  KjjjJq),  or  effective  stress  intensity  factor  range  of  (AK^ff  =  K^^j^  -  Kj,j). 
The  formulation  is  an  expansion  of  the  Paris  Law  relationship  of  crack  growth  rate,  such  that 
along  the  specimen  free-surface: 

dc/dN  =  Ce(AKeeff)"c  (7.5) 

and  into  the  specimen  thickness: 

da/dN  =  Ca(AKaeff)“a  (7-6) 
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where  and  n^,  are  Paris  Law  constants  generated  from  FCGR  data  along  the  crack  free* 
surface,  and  and  n^  into  the  crack  thickness.  The  constants  can  be  generated 
independently  when  da/dN  and  dc/dN  data  are  available  --  i.e.,  Figures  7.8  through  7.11. 
Although  the  analytical  formulation  was  set-up  to  recognize  these  potential  Paris  Law 
constant  differences  into  the  specimen  thickness  and  along  the  crack  firee-surface,  this 
investigation  used  Cj^  =  C^.  and  n^  =  n^  since  they  were  approximately  equal  (Figures  7.12 
and  7.13).  This  assumption  is  typically  made  for  surface  flaws  in  metals  even  though  non- 
homogenous  condition  are  known  to  exist.  For  example,  Hodulak  (50)  suggested  that  for 
metal  plates  in  tension,  the  material  resistance  is  greater  at  the  surfaces  than  at  the  plate  mid¬ 
thickness. 

The  FCG  calculation  approach  is  based  on  an  iterative  crack  growth  model  described 
in  (SI)  and  (52),  and  follows  from  the  relationship  of; 

AN  =  Aa/ (da/dN)  (7.7) 

By  assuming  a  "Aa"  value  and  knowing  da/dN  from  an  applicable  crack  growth  model  such 
as  the  Paris  Law  or  Walker  relationships,  AN  can  be  calculated  to  match  the  assumed  Aa 
increment.  Kg  and  are  then  calculated  for  the  new  crack  size,  and  the  iterative  process 
repeated  until  a  desired  crack  size  is  reached.  An  error  term  is  built  into  the  routine  to  ensure 
only  small  increments  of  AN  are  calculated  even  when  the  crack  rate  is  high.  This  error 
function  was  established  to  maintain  efficient  conputational  speed  and  accuracy. 

The  fatigue  crack  growth  (FCG)  prediction  code  was  validated  using  thru-thickness 
specimen  data  generated  in  conjunction  with  an  ASTM  "Round  Robin"  (53)  evaluation  task 
group  E24.06.01  on  comer  cracked  holes.  During  the  "Round  Robin"  evaluation,  predictions 
were  made  by  representative  industry,  academia,  and  government  organizations.  The  error 
range  of  prediction  to  experimental  results  is  as  follows: 
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0.68  «  Life  ratio  (Np/Ng)  »  1.18 

0.84  «  Crack  shape  ratio  ((a/c)p/(a/c)g)  »  1.19 

where  Np/Ng  is  the  ratio  of  predicted  ("p”)  to  experimental  (''e'')  cycles.  Likewise,  for  the 
crack  shape  ratio  (a/c)p/(a/c)g,  Table  7.6  compares  predictions  m a  Je  by  the  code  developed 
for  this  investigation  with  the  "Round  Robin"  experimental  results.  These  results  are  all  well 
within  "Round  Robin"  prediction  ranges  shown  above. 


Table  7.6:  Comparison  of  fatigue  crack  growth  predictions  using  a  "Round  Robin" 
experimental  data  set  and  the  Newman-Raju  K-solutions 


R 

S^J,y_(k^ 

0.1 

15 

0.1 

20 

0.3 

15 

ErMf. 

(a./c)^(a/c)^ 

0.91 

0.89 

0.85 

0.94 

1.00 

1.02 

The  fatigue  crack  growth  (FCG)  computational  program  was  then  used  to  predict 
crack  growth  patterns  and  cycles  to  achieve  a  given  crack  length  for  ihis  investigation. 
Calculations  were  made  using  five  different  combinations  of  Paris  Law  constants  (C  and  n) 
and  closure  load.  The  two  uncorrected  constants  came  from  thru-thickness  data  (37),  and 
surface  flaw  data  developed  in  this  investigation  (Figure  7.7).  Two  of  the  corrected 
constants  were  developed  from  Figures  7.12  and  7.13  using  bPop  and  cPop,  respectively; 
that  is,  one  calculation  used  eff  ^*^^c  effi  eff  eff 

The  final  calculation  used  bPop  for  the  da/dN  and  cPop  for  the  dc/dN  calculations,  or 
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and  --  as  discussed  in  Chapter  VI.  The  constant  AK^  test  results  were  used  for 

these  sample  calculations  after  eliminating  precracking  effects. 

The  predicted  crack  growth  patterns  (a/c  and  a/t)  and  cycles  to  achieve  a  given  crack 
length  are  compared  to  experimental  results  in  Table  7.7,  Here,  the  cycle  prediction 
improvements  of  using  AKgff  (where  Kg|  varies  with  crack  size,  ref.  Figures  7.1  through 
7.4)  versus  an  applied  AK  are  apparent.  However,  crack  growth  pattern  results  are  not  as 
apparent,  as  shown  in  Figure  7.14.  Several  of  the  calculations  are  equivalent:  (1)  the  two 
calculations  using  the  nominal  applied  AK  for  thru-thickness  and  surface  flaws,  and  (2)  the 
two  calculations  using  bPop  (AK^^  gff  and  AK^g  gff)  or  cPgp  g^^  and  AK^g  gff)  for 
both  da/dN  and  dc/dN.  Although  all  of  the  combinations  are  quite  good,  using  bPgp  for 
da/dN  (AK^j  gf^)  and  cPgp  for  dc/dN  (AK^^g  gff)  provided  the  best  overall  match.  In 
general,  all  of  the  analytical  results  using  effective  load  range  provide  a  better  match  with 
experimental  data  than  the  nominal  applied  load  range. 
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Table  7.7:  Comparison  of  Experimental  and  Predicted  Fatigue  Crack  Growth  Patterns  and 
Cycles  Using  Uncorrected  and  Corrected  Paris  Law  Constants 


Condition 

a/t 

(a/cV(a/cL 

NnZN^. 

Apnlied  Load  fAK) 
Thru-Thickness 

0.1 

1.10 

0.16 

0.2 

1.00 

0.24 

0.4 

0.87 

0.23 

Surface  Flaw 

0.1 

1.09 

0.53 

0.2 

0.99 

0.77 

0.4 

0.84 

0.77 

Effedtive  Load  Ranee  f  AIC.ffl 

for  Surface  Flaws 

bPop 

0.1 

1.09 

0.70 

(K^effAKbeeff) 

0.2 

1.01 

0.97 

0.4 

0.91 

0.89 

cPop 

0.1 

1.09 

0.50 

(K'=aeff&K‘=oeff) 

0.2 

1.01 

0.79 

0.4 

0.92 

0.95 

bPop/cPop 

0.1 

1.02 

0.56 

(Kbaeff&KCceff) 

0.2 

0.97 

0.82 

0.4 

0.98 

0.87 

Experimental  Results: 

aZt 

(a/c')^ 

Npl'cvcles'l 

0.1 

0.72 

57,000 

0.2 

0.68 

140,000 

0.4 

0.54 

480,000 
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14  Comparisons  of  predicted  and  experimental  crack  growth 
patterns  for  die  constvit  AKd  test 


The  results  presented  in  Figure  7.14  would  indicate  that  the  fully-open  K-solutions 
are  sufficiently  accurate  to  provide  good  crack  growth  pattern  predictions  for  various 
combinations  of  Paris  Law  constants  when  using  detailed  closure  load  measurements.  The 
questions  then  arise  as  to  how  good  the  predictions  are  if  only  Paris  Law  constants  are 
available,  and  one  assumes  a  constant  (not  a  function  of  crack  length)  closure  load  ratio 
(Pgj/Pjjj3jj).  Results  of  this  approach  are  shown  in  Table  7.8,  where  applied  AK  Paris  Law 
constants  are  shifted  20%  and  30%  for  thru-thickness  (37)  and  surface  flaw  (Flgur'*.  7.7)  tests 
data.  Two  combinations  of  constant  closure  load  ratio  are  evaluated;  (1)  30%  for  both 
da/dN  and  dc/dN,  and  (2)  20%  for  da/dN  and  30%  for  dc/dN  (Figure  5.19). 


Table  7.8:  Comparison  of  Experimental  and  Predicted  Fatigue  Crack  Growth  Using  an 
Assumed  Constant  Pcl/^max  Value  for  Surface  Flaws 


Condition 

a/t 

ra/cV(aycV 

Pcl'^ max  ~ 

0.1 

1.07 

0.42 

0.2 

0.98 

0.68 

0.4 

0.87 

0.89 

Pcl^^max  ~  0.2/0.3 

0.1 

1.09 

0.52 

0.2 

0.99 

0.76 

0.4 

0.84 

0.77 

It  is,  therefore,  observed  from  these  calculations  that: 

(1)  Existing  crack  growth  pattern  predictions  are  quite  accurate. 
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(2)  Cyclic  predictions  using  an  assumed  constant  closure  load  are  reasonable, 
but  not  as  accurate  as  using  measured  variable  closure  load  data  — 

Table  7.7. 

(3)  The  best  approach  for  predicting  FCG  of  surface  flaws  in  bending  when 
only  hilly-open  K-solutions  and  closure  load  data  is  not  available  is  to:  (1) 
use  Paris  Law  constants  generated  from  thru-thickness  specimens  data  and 
assume  a  constant  closure  load  ratio  which  is  based  on  the  general  material 
class  being  analyzed. 
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VIII.  CRACK  OPENING  DISPLACEMENT  PREDICTIONS 
USING  WEIGHT  FUNCTIONS 


One  of  the  few  analytical  solutions  available  for  predicting  surface  flaw  crack 
opening  displacement  (COD)  was  developed  by  Mattheck,  Morawiety  and  Munz  (54).  Their 
approach  is  based  on  a  weight  function  method  developed  by  Petroski  and  Achenbach  (55)  to 
calculate  COD  of  thru-thickness  cracks.  The  linear-elastic  model  was  developed  for  applied 
loads,  and  does  not  reflect  closure  load  effects.  Therefore,  the  predicted  elastic  displacement 
is  a  function  of  the  applied  stress  intensity  fector. 

The  weight  function  mathematical  formulation  was  originally  introduced  by 
Bueckner  (56,57)  in  1970  to  calculate  stress  intensity  fector  for  complex  loads.  The 
formulation  uses  superposition  to  calculate  the  stress  intensity  factor  (Kq^^)  for  a  new  stress 
distribution  (a^ew)  ^hen  the  mode  I  stress  intensity  factor  (K)  and  elastic  displacement  are 
known  for  a  symmetrically  loaded  body.  The  weight  function  formula  is  given  as: 

Knew  =  (H/Kr)  Onew(x)  [  aur(a,x)  /  da]  dx  (8.1) 

where  x  is  the  distance  from  the  crack  free-surface  into  the  crack  depth  (Figure  5.1), 

H  =  E/(l-  u^)  for  plane  strain  or  E  for  plane  stress  conditions.  K^  and  u^  are  known 
reference  stress  intensity  factor  and  displacement,  respectively,  for  a  simple  symmetrical 
loading.  The  weight  function  relationship  has  been  used  successfully  by  a  number  of 
investigators  such  as;  Rice  (58),  Grandt  (59-63),  Paris  (64),  and  Labbens  (65,66). 

The  uniqueness  of  the  Petroski  and  Achenbach  weight  function  approach  for 
predicting  COD  is  their  formulation  of  a  reference  crack  displacement  function  Uj..  The 
proposed  displacement  is: 
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ui(a,x)=  [00/V2H]  [4F(a)v'aV(a-x) 
+  G(a)(a-x)3/2/Va] 


(8.2) 


where  F(a)  is  a  known  boundary  correction  factor,  and  G(a)  is  unknown.  Given  a  reference 
stress  (Of)  and  stress  intensity  factor  (Kr)  for  a  reference  displacement,  equation  8.1 
becomes: 

=  H  io®  Or(x)  [  dui.(a,x)  /  da]  dx  (8.3) 

Since  Ur(a,x)=0  when  x=a,  the  partial  derivative  can  be  moved  outside  the  integral,  and  the 
integral  expression  reduced  to: 

4*  Kr^  da  =  H  Or(x)  Ur(a,x)  dx  (8.4) 

Next,  the  reference  stress  intensity  factor  takes  the  familiar  form  of: 

Kr  =  Oq  Vna  F(a)  (8.5) 

and  solving  for  the  unknown  G(a)  term  by  substituting  equation  8.2  into  equation  8.4.  Then, 
for  a  constant  uniform  stress  field  o^: 

i/  Kr^  da  =  n  Oq^  J/  F2(a)  da 

=  [oo  /  >^2]  [  4F(a)  Va  //  Oj^x)  V(a-x)  dx 

+  G(a)/  Vsi  /q®  ar(x)  (a-x)^/^  dx]  (8.6) 

Letting: 

ll(a)  =  Ji  Oq  V2iQ®  aF2(a)  da 

12(a)  =  4*  Or(^)  >^(a-x)  dx  (8.7) 

13(a)  =  4*  ^r(x)  (a-x)3/2  dx 
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and  solving  for  G(a): 


G(a)  =  [  Ii(a)  -  4F(a)  v^a  12(a)]  Va  /  13(a)  (8.8) 

Next,  a  standard  formulation  can  be  derived  for  the  integrals  in  equation  8.7  for 
various  load  distributions.  For  a  constant  reference  load  where  =  Oq  (tension): 

12(a)  =  V  ^(®‘*) 

13(a)  =  CTq  (a-x)^^  dx  =  2/5  Oq  a^^  (8.9) 

Using  a  general  stress  distribution  equation  presented  by  Yagawa  (67)  to  represent  pure 
bending  in  a  plate: 

=  Oq  oj  (1  -  x/b)  (8.10) 

where  b  is  the  plate  half  width.  Then: 

*2(a)  =  io*  [^o  +  01  (1  -  x/b)]  (a-x)  dx 
=  2/3  a3/2  [oo  +  01  (1  -  2/5  a/b)] 

13(a)  =  io*  t^o  +  01  (1  -  xfo)]  (a-x)3/2  dx 

=  2/5  a5/2  +  01  (1  -  2/7  a/b)]  (8.1 1) 

Then,  using  equation  8.8  to  calculate  G(a),  the  crack  opening  displacement  (COD)  can  be 
determined  from  equation  8.2. 

The  above  Petroski  and  Achenbach  crack  opening  displacement  formulation  for  thru- 
cracks  was  expanded  by  Mattheck,  et  al.  to  three-dimensional  surface  cracks  under 
symmetrical  loading.  The  surface  flaw  solution  is  separated  into  three  regions  of  analysis; 
Uf(x,y=:0),  Uj(x=0,y)  and  Uf(x,y).  TTie  first  region  (x,y=0)  was  solved  by  Petroski  and 
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Achenbach  as  described  above.  The  second  region  (x=0,y)  uses  the  Griffith  crack 
relationship  of: 

Uj(x=0,y)  =  [1  -  (y/c)2  (8.12) 

where  Uj.®“  “  Uy(x=0,y=0)  from  equation  8.2.  Titus,  the  first  two  crack  regions 
displacement  are  described  for  (x,y=0)  and  (x=:0,-c<y<+c). 

For  the  complete  crack  area  (x,y),  the  crack  depth  profile  can  be  defined  in  a  similar 
maimer  as  the  crack  free-surface  where: 

a*  =  a  [1  -  (y/c)2]l/2  (8.13) 

ur(x,y)  =  [  oo/  (V2  H)]  [4Fo(a*)  V(a*)  V(a*-  x) 

+  Go(a*)  (a*-  x)3/2/  V(a*)]  (8.14) 

To  match  the  surface  displacement  profile,  the  unknown  GQ(a*)  is  defined  as: 

Go(a*)  =  L2  H  u  /  (a  Oq)]  -  4  Fo(a*)  (8.15) 

The  above  formulation  is  a  very  useful  approach  for  predicting  three-dimensional 
surfece  flaw  COD  because  of  its  flexibility  to  utilize  sophisticated  stress  intensity  factor 
solutions  and  complex  stress  fields.  The  paper  by  Mattheck,  et  al.  presented  a  closed  form 
solution  for  tension  (eqn.  8.10)  using  simple  K-solution  developed  by  Dufiesne  (54). 

This  investigation  developed  a  closed  form  solution  for  a  plate  in  bending  (eqn. 

8.11),  and  utilized  the  more  complex  Newman-Raju  K-solutions.  In  addition,  a  numerical 
computation  program  was  developed  (Appendix  A)  to  calculate  surface  flaw  COD  of  more 
complex  stress  distributions  using  an  eight  order  Newton-Cotes  integral  formulation  and  the 
Newman-Raju  K-solutions.  The  program  closed  form  solutions  for  tension  and  bending,  and 
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the  numerical  integration  routines,  were  validated  by  comparing  the  calculations  to  published 
results  (54,68,69). 


COD  Predictions  Versus  Measurements 

Comparisons  between  experimental  data  and  weight  function  (WF)  COD  predictions 
are  presented  in  this  section.  Perez  (70)  made  comparisons  between  numerous  displacement 
models,  and  thus,  this  section  will  concentrate  on  the  Petroski-Achenbach  formulation.  The 
predictions  are  made  at  a  Poisson  ratio  of  0.4,  and  average  elastic  modulus  of  4.E5  psi  (71). 

Figure  8.1  is  a  plot  of  experimental  and  predicted  COD  profiles  into  the  crack  depth 
(Figure  5.15,  between  locations  A  and  B)  at  three  applied  load  levels  —  bP^p,  tP^p  and 
tPop+.  In  this  figure,  "r"  is  the  distance  from  the  crack  tip  and  "a"  the  crack  depth.  The 
bPop  load  level  (22%  cf  Pmax)  ***  closure  load  region,  and  is  the  load  required  to  open 
the  crack  free-surface  at  location  B.  The  tPop  load  level  (30%  of  Pmax)  defined  as  the 
transition  load  where  the  "void"  maximum  displacement  is  at  the  crack  free-surface  (location 
B).  Loads  greater  than  tPgp  are  considered  to  be  in  the  COD  elastic  region,  and  loads  less 
than  tPop  are  in  the  closure  load  region.  The  third  load  level  (38%  of  Pmax)  ^  within  the 
elastic  load  range,  and  is  designated  as  tPgp-f  in  Figure  8.1. 

At  the  crack  free-surface  (r/a=1.0)  in  Figure  8.1,  the  WF  model  is  found  to 
underpredict  the  displacement  at  the  tPQp+  load  level,  match  it  quite  favorable  at  tPgp,  and 
poorly  predict  displacement  at  bP^p.  At  the  bP^p  load  level,  the  crack  fi«e-surfece  is  just 
beginning  to  open  while  the  elastic  model  predicts  a  displacement  of  0.9E-4  inches.  Note 
that  the  displacement  points  at  the  crack  free-surface  are  extrapolations  of  measured  data 
(Figure  5.16). 

The  ^dieted  displacement  profiles  show  a  fevorable  agreement  near  the  crack  tip 
(r/a  <  0.05)  for  all  three  load  levels.  At  the  lower  bP^p  load  level,  the  elliptical  profile 
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proposed  by  the  Griffith  formulation  (eqn.  8.13)  follows  the  experimental  data  for  r/a  <  0.2, 
but  then  deviated  due  to  inelastic  closure  at  the  crack  free-surface.  At  tPQp,  the  predicted 
displacement  matches  the  experimental  data  at  the  crack  free-surface,  but  does  not  match  the 
profile  for  0.05  <  r/a  >  0.8.  At  tPop+,  the  predicted  displacement  profile  only  matches  the 
experimental  measurements  near  the  crack  tip  (r/a  «■  0.0). 

The  prediction  sensitivity  to  modulus  of  elasticity,  which  ranges  from  3.5E5  to  4.9E5 
psi  for  PMMA  (69),  can  be  calculated  directly  since  displacement  is  inversely  proportional  to 
modulus.  For  this  modulus  range,  maximum  displacement  ranges  from  1.46E-4  to  1.04E-4 
at  the  tPop  load  level.  This  range  is  not  sufficient  to  account  for  differences  between 
predicted  and  experimental  data.  The  maximum  displacement  is  directly  proportional  to 
applied  load,  and  is  not  effected  by  Poissons  ratio.  Variations  in  Poissons  ratio  were  found 
to  have  minor  influences  on  displacement  profile  distribution,  but  not  of  sufficient  effect  to 
match  the  experimental  data. 

The  data  in  Figure  8.1  can  also  be  presented  in  a  load  versus  displacement  format  as 
shown  in  Figure  8.2  at  the  crack  free-surfrice  (location  B).  Here  the  predicted  elastic 
displacement  emanates  from  the  origin  and  passes  through  the  tPop  data  point  In  the  closure 
load  region  (P/Pmax  <  ^*^)>  linear  elastic  fracture  mechanic  (LEFM)  predictions  are  not 
expected  to  match  experimental  data,  and  significantly  overpredict  the  displacement  In  the 
elastic  COD  region  (F  max  >  0.3),  the  extrapolated  experimental  displacement  is  found  to 
be  greater  than  predicted  at  the  crack  free-surface.  Even  though  the  crack  free-surfece  is 
extrapolated  data,  the  measured  displacement  internal  to  the  crack  at  r/a  =  0.4  is  noted  to 
have  a  greater  displacement  than  predicted  at  the  crack  firee-surface.  (Figure  8.1).  At  Pmax’ 
these  differences  represent  a  factor  of  3  error  between  experimental  and  prediction 
displacement  at  the  free-surface. 
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LOAD  -  DISPLACEMENT  TRENDS 
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Displacement  Profile  Model  Modifications 


Differences  between  experimental  data  and  elastic  predictions  can  be  attributed  to 
inelastic  effects  not  accounted  for  in  the  model  particularly  in  the  closure  load  region.  Other 
modeling  limitations  ate  generated  from  the  assumed  elliptical  profile  distribution  not 
matching  the  experimental  displacement  profile  shape.  This  can  be  seen  in  Figure  8.1.  One 
possible  approach  is  to  modify  the  shape  factors  associated  with  the  Petroski-Achenbach 
relationship  (eqn.  8.3). 

Since  the  weight  function  model  is  elastic,  improvements  to  the  displacement  profile 
should  be  evaluated  in  the  elastic  region  (P/Pmax  >  Figure  8.2.  Near  the  crack  tip,  the 
shape  factor  term  proposed  by  Paris  and  Sih  (72): 

u(a,x)  =  4KyH  [(a-x)  /  V2J  (8. 16) 

appears  to  match  the  experimental  data.  However,  the  predicted  amplitude  is  low  by  a  factor 
of  approximately  1.4  at  r/a  =  0.1.  This  could  be  improved  by  using  a  model  other  than  the 
current  elliptical  profile  for  r/a  >  O.OS. 

In  addition,  the  "G(a)"  term  shape  factor  could  be  modified  to  reflect  experimental 
results.  As  stated  in  the  Petroski-Achenbach  publication  concerning  the  higher  order  terms 
in  (a-x),  "The  choice  of  subsequent  terms  is  largely  a  matter  of  judgment,  the  limitation 
being  that  only  one  unknown  function  can  be  determined  from  eqn.  (8.4)." 
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IX.  SUMMARY,  CONCLUSIONS  AND 


RECOMMENDATIONS 


Summary  and  Conclusions 

The  objective  of  this  investigation  was  to  determine  the  effect  of  fatigue  crack 
closure  on  the  growth  of  surface  flaws  in  bending.  A  transparent  polymethylmethacrylate 
(PMMA)  specimen  and  a  Newton  interferometry  system  were  used  to  map  the  crack  opening 
displacement  (COD)  patterns  for  various  crack  sizes.  Crack  closure/opening  load 
measurements  were  made  concurrent  with  the  COD  interference  measurements. 

Interference  fringe  order  photographs  were  made  at  crack  growth  increments  (Ac)  of 
approximately  0.05  inch  for  surface  crack  lengths  "c"  between  0.02  and  0.60  inch  long. 
These  measurements  were  duplicated  for  four  different  load  cycles  ~  two  constant  stress 
intensity  factor  (AK)  tests  (generated  by  measuring  "a"  and  "c"  crack  lengths  and  load 
shedding),  a  constant  load  and  a  block  loading  test.  COD  interference  photos  were  digitized 
using  a  photo  enlarger  and  digitizing  table.  Fatigue  crack  growth  rate  correlations  were 
made  for  various  definitions  of  effective  stress  intensity  factor  (AKgff).  Fatigue  crack 
growth  predictions  are  compared  to  experimental  cycles  and  aspect  ratios  using  AK  and 
AKgff  values.  COD  measurements  were  compared  to  analytical  predictions  using  a  weight 
function  approach  for  surface  flaws. 

The  experimentally  measured  crack  opening  displacement  (COD)  profiles  and  visual 
observations  identified  three  patterns  which  directly  influence  closure  load  level  and  crack 
tip  stress  intensity  factor  in  the  nonlinear  closure  region.  These  three  distinct  patterns  are 
referred  to  as  crack  types  I,  II  and  III  in  this  investigation.  The  COD  patterns  (i.e.  crack 
types)  were  found  to  be  a  function  of  crack  length  (a/t),  but  not  load  history.  Crack  types  II 
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and  III  COD  patterns  are  distinguishable  by  a  "void"  area  which  separates  the  crack  interior 
surfaces,  while  the  crack  ffee*surfeice  is  closed  under  zero  load.  As  a  load  is  applied  to  the 
specimen,  these  two  crack  types  begin  to  initially  open  at  the  interior  crack  tip,  and  then 
open  outward  toward  the  crack  ffee-surface.  At  the  crack  free-surface,  its  middle  begins  to 
open  first,  and  then  continues  to  open  until  the  full  surface  crack  is  open.  The  last  portion  of 
the  crack  to  open  is  not  at  the  free-surface,  but  along  the  crack  internal  surfaces,  12  to  15 
degrees  from  the  crack  free-surface.  This  surface  flaw  crack  tip  opening  pattern,  which 
occurs  in  the  closure  load  region,  is  referred  to  as  a  "partially-open"  condition  since  the  crack 
tip  boundary  is  open  at  certain  locations  while  it  is  closed  at  others. 

Closure  load  measurements  were  made  at  three  distinct  locations  on  the  crack 
surface.  These  three  locations  correspond  to  closure  load  measurements  when  (1)  the  crack 
surface  is  completely  open,  (2)  opening  at  the  crack  middle  free-surface,  and  (3)  opening  at 
the  deepest  portion  of  the  crack  tip  boundary.  The  closure  load  ratios  (Pcl^max)» 
inherent  closure  load  stress  intensity  factors  (Kc|),  were  found  to  be  a  function  of  both  crack 
size  and  load  history.  Qualitatively,  the  measured  crack  closure  ratios  compared  favorably  to 
published  experimental  and  finite  element  analytical  results. 

As  the  cracks  grow  in  size  due  to  fatigue  cycling,  the  "void"  area  also  increases  in 
size,  causing  a  decrease  in  closure/opening  load  level  and  an  increase  in  crack  displacement 
at  zero  load.  These  "void"  crack  opening  displacement  patterns  had  a  direct  impact  on  crack 
tip  boundary  in  the  closure  load  region.  There  are  no  known  K-solutions  for  this  type  of 
"partially-open"  crack  tip  boundary.  Available  K-solutions  used  in  this  investigation  are 
referred  to  as  "fully-open"  as  opposed  to  "partially-open"  geometric  COD  patterns  which 
were  measured. 

It  is  hypothesized  from  the  COD  patterns  and  crack  growth  rate  data  that  available 
"fully-open"  K-solutions  over  predict  the  actual  value  (AK  values  into  the  crack  depth) 
when  the  minimum  applied  load  is  in  the  closure  load  region.  This  difference  is 
hypothesized  to  be  caused  by  available  K-solutions  not  modeling  the  "partially-open"  crack 
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tip  boundary  developed  by  the  "void"  formation  internal  to  the  crack  surfaces.  This  is 
considered  a  conservative  discrepancy  in  that  "actual"  is  less  than  available  "fully-open" 
K-solutions  predict.  This  leads  to  higher  predicted  crack  growth  rates  and  fewer  cycles  to 
failure.  The  error  associated  with  not  having  a  "partial ly-open"  K-soIution  appears  to  vary 
with  the  number  of  cycles  to  achieve  a  given  crack  length.  That  is,  for  shorter  a/t  values,  the 
error  is  estimated  to  be  a  factor  of  2.  However,  for  larger  a/t  values  the  error  is  less  since  the 
crack  COD  pattern  better  represents  a  "fully-open"  condition,  and  the  error  is  estimated  to  be 
a  factor  of  1.2. 

The  correlation  of  fatigue  crack  growth  rates  with  several  definitions  of  AK^ff  (using 
closure  load  measurements  and  available  "fully-open”  K-solutions)  found  the  crack  middle 
free-surface  closure  load  measurement  (bPop)  produced  the  least  scatter  when  ''ompared 
with  experimental  results.  Complete  opening  of  the  crack  surfaces  (cPop)  also  produced 
good  correlation.  An  analytical  evaluation  of  ftitigue  crack  growth  using  several  AK^ff 
definitions  found  that  (1)  crack  growth  pattera  (a/c  and  a/t)  results  are  quite  accurate 
whether  closure  load  corrections  are  used  or  not,  and  (2)  fatigue  cycle  predictions  are 
improved  when  using  closure  corrections.  These  results  indicate  that  reasonably  accurate 
crack  growth  predictions  can  be  made  using  "fully-open"  K-solutions.  When  measured 
closure  load  data  are  not  available,  fatigue  crack  growth  predictions  which  assume  a  constant 
closure  load,  and  are  representative  of  the  material  being  considered,  are  shown  to  provide 
reasonably  good  results.  The  best  match  of  crack  growth  pattern  and  cycles  predictions  and 
experimental  measurements  occurs  when  a  bimodal  analysis  was  used.  That  is,  when  two 
different  closure  loads  were  used  in  the  calculation  —  one  for  crack  growth  into  the  specimen 
depth  (bPop  with  da/dN),  and  a  different  one  along  the  crack  free-surface  (cPop  with  dc/dN). 

Evaluation  of  the  elastic  weight  function  approach  for  predicting  COD  for  surface 
flaws  found  the  prediction  to  underestimate  the  data.  A  primary  cause  for  this  error  is 
associated  with  the  linear  elastic  model  not  matching  the  nonlinear  COD  pattern  in  the 
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closure  load  region.  Additional  error  occurs  from  the  assumed  elliptical  COD  profile,  which 
does  not  match  experimental  data  taken  at  loads  above  the  closure  load  level. 
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Recommendations 


The  following  recommendations  are  offered  as  suggestions  for  further  research: 

1.  There  is  a  need  to  develop  "partially open”  K-solutions  for  surface  flaw  and  thru- 
thickness  cracks  in  the  closure  load  region.  These  solutions  are  required  for  the  complete 
crack  tip  boundary,  and  would  allow  the  use  of  K-solutions  which  match  geometric  CX)D 
profiles  experienced  by  the  crack  tip  —  versus  "fiilly-open"  K-solutions  currently  in  use. 

2.  Develop  a  crack  opening  displacement  model  which  matches  the  "void"  formation 
in  the  closure  load  region. 

3.  Conduct  an  experimental  program  which  measures  crack  opening  displacement 
(COD)  in  the  elastic  load  region  in  order  to  improve  the  COD  prediction  models.  These  data 
would  provide  the  elastic  displacement  profile  shape  factors  -  as  opposed  to  the  inelastic 
closure  load  region  described  in  #2  above.  The  use  of  a  higher  frequency  monochromatic 
light  source  than  was  used  in  this  investigation  is  required  to  generate  the  increased  fringe 
order  resolution  at  the  higher  load  levels.  For  example,  by  reducing  the  number  of  fringe 
orders  by  half  (for  a  typical  crack  size  which  begins  to  open  in  the  elastic  region)  requires  a 
wavelength  of  approximately  12,000  A,  versus  the  6328  A  laser  light  source  used  in  this 
investigation.  Since  a  light  source  is  not  known  to  be  available  at  this  high  frequency,  light 
filtering  techniques  would  appear  to  be  the  most  reasonable  approach  for  this  type 
investigation. 

4.  A  need  exists  to  experimentally  evaluate  plasticity  induced  closure  effects  on  the 
plane  strain  constraint  [a^  /  (02+03)]  near  the  crack  free-surfsce.  Trantina,  et  al.  (43) 
identifies  plasticity  as  a  primary  factor  for  defining  plane  stress  and  strain  plastic  zones  size 
influences  on  the  plane  strain  constraint  Using  a  3-D  elastic-plastic  finite  element  analysis 
of  surfade  flaws,  the  authors  showed  a  rapid  loss  in  constraint  for  <^<15  degrees  which  was 
not  identified  using  elastic  solutions.  These  plane  stress/strain  effects  were  also  noted  in  this 
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Appendix  A  Constant  AKa  Test  Data 

A  complete  listing  of  the  experimental  data  is  presented  in  this  appendix  for  the 
constant  test.  The  data  sheet  column  variables  are  follows:  the  applied  delta  bending 
moment  in-lbs.),  "a"  and  "c"  crack  lengths,  a/c,  a/t,  accumulated  fatigue 

cycles,  AK,.,  plastic  zone  size  at  location  C,  AK^,  plastic  zone  size  at 

location  A,  cPop,  bPop  and  apQp  applied  load  (Ib).  Negative  one  (-1.0)  is  a  designator  to 
indicate  that  closure  load  data  were  not  recorded  at  that  crack  size.  Blank  lines  between  the 
data  separate  categories  of  precracking,  type  I,  type  II  and  type  III  cracks. 
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Appendix  B  Constant  AKd  Test  Data 


A  complete  listing  of  the  experimental  data  is  presented  in  this  appendix  for  the 
constant  test  The  data  sheet  column  variables  are  as  follows:  the  applied  delta  bending 
moment  -  Mmi^,  in-lbs.),  "a"  and  "c"  crack  lengths,  a/c,  a/t  accumulated  fatigue 

cycles,  AK^  plastic  zone  size  at  location  C,  AK^  plastic  zone  size  at 

location  A,  cP^p,  bpQp  and  aP^p  applied  load  (lb.).  Negative  one  (-1.0)  is  an  indication  that 
closure  load  data  were  not  recorded  at  that  crack  size.  Blank  lines  in  the  data  separate 
categories  of  precracking,  type  I,  type  II  and  type  III  cracks. 
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Appendix  C  Constant  Load  Test  Data 


A  complete  listing  of  the  experimental  data  is  presented  in  this  appendix  for  the  constant 
load  test.  The  data  sheet  column  variables  are  as  follows:  the  applied  delta  bending  moment 
(^max  ■  ^min>  and  "c"  crack  lengths,  a/c,  a/t,  accumulated  fatigue  cycles,  AK^., 

Kj.  jjjjjjj,  plastic  zone  size  at  location  C,  AKa,  Plastic  zone  size  at  location  A,  cPop, 

bPop  and  aPop  applied  load  (lb.).  Negative  one  (-1.0)  indicates  that  closure  load  data  were 
not  recorded  at  that  crack  size.  Blank  lines  in  the  data  separate  categories  of  precracking, 
type  I,  type  II  and  type  III  cracks. 


208 


SPeCtMlN  3-15  CilNStMIt  LOAD  ItSI  lOttIA  NOHCMI  •  A««  IN-LMI  AT  0  AC6ACCS 


^  9^^  mm9^  mmmmm 

mmm  MMttMMMM  mm 

I  1  •  1  1  1  1  1 

iJiZ 

O  1  1 

O  ••  A*  K  A*  A* 

Ml  1  g|  1  1 

M«*  0«*M0MM 
MO  ••Mm«»MM 

OM  O 

MM  » 
MMM 

MMMM  O  M 

M  mm  ••••MMMM««M«a«»M«»«aMM  MMM  ^MMM  M«» 

B  MM  ••#«<«»•#««  MM  MM 

mmmmmmmmmmmmmmm  mmm mmmmmmm mm 
u 

mMO 

•  •  • 
mmm 

0«»M 

••MM 

OMMa«OM 

M  M  M  M  M  •• 

l•ll(lTlllt•ll  ltllt••t••i• 
ummutijimmmmmmiitmM  mmiitmmmmm>^m  ^nt 

MM  MMMMMW 
OO  OOOO  OO 

i  1  1  1  1  •  1  1 

MMM 

ooo 

•  1  1 

SS5 

«•••  •• 

•  •  • 

OOO 

1  •  1  1  •  1 

rs8sx» 

MMMOM  M 

M 

MMMMMMOOMOmamMMMmMMMMMMMMMmM 

S  mm  MM  MMMMMMMM  MM  MMMMMMMMMM  MM 

M 

M 

oom^mmoo 

MO  M 

MM  M 
•mM  M 
MK  M 

K  K  M  M  M  M 

MO  MKOMKO 
MMMMOMOO 

MM  O 
•  •  • 

OM  M 
MO  M 
MMM 

M  M  O  M  ••  O 

o  ^ 

M  |(il|lliitltttiiiit«|tli|l 

Mr«MMMMMMMMMMM«aM««MMMMMBMM«*MM 

^  cmmmmmomommmmmmmmm  mmm  mmmm  mm 

M  MM  M  M  «tM  MM  M  M  mm  mm  M  MM  «•  MM  M  MM  M  MM 

mm  mmmmmm 
oo  OOMOOO 

1  1  *  1  1  1  1  1 

inOO«*OMMO 
MO  MOoMOM 

MM  mmmmmm 

MM  M 
MO  M 

»  t  1 

M«»  M 

MM  M 

A  •  • 

OO  o 

I  i  i  1  1  1 

M  M  M  W  W  W 
M  M  n  M  M  M 
M  MM  O  M  O 
M  M  M  M  n  O 

ooo 

M  at  MM  M  M  MM  MM  M  M  M  M  MM  MMM  MMM  MMM  M  MO 
M  «  MM  MMMMMMMMM^m^MMMMMMMMMMMM 

M  sommmmmmmmmmmkmmmmmmmmmmmmm 

i 

MO  MMmOOO 
MM  MMO^MM 

MO  O 
OO  M 

OK  O 

••••OMOM 

OOMOOO 

S  M  MM  mMMMMMMMMMMMMMMMMMm^^MMMO 

M  ^MMMMpiIMMMMMMOMOMMMMMMmMMMMM 

s  mommmmmmmmmmmmmmmmmmmmmmmmm 

O  M 
s 

MM  ••OOMMM 

OO  mMMOmM 

MM  O 

A  •  » 

OOO 
MK  M 

KM  M 

AA  aE 

^  MOM9MO MMO MOOOO MOM OMOmOM M OO 
•i  ^OMOOOMOOMMOOMOOMMOOMoOOOeO 
»  riM  o  •  MO  OO  O  M  MO  MO  MMM  OMO  MOM  M  M  o 
a  umkmMmo^mmmmmmmmmoommmmmmomm 
V  MM^MvaMMMMMMMiilM  MMM  OMM  ^MMM  MM 

OMMOMOMO 

O  O  MOOO  MO 
OMMOMOMO 
MMMOO«*MM 

OO  o 

A  •  • 

ooo 

MOM 

SSSSES 

^  MM  M  M  MM  MM  M  M  M  O  MO  M  MM  OMO  MOM  M  MM 

M  mmMMOMOMOMOOMmMM  mm  MMM  mMM  O  MM 
M OMMMMMM MM M MMMmM MMM M mm MOM M MM 
M  OOOMOMO OMMOMOMO MOO OMOmOM M MO 

MO  mmmmmm 
MM  OOMOOO 

MO  K 

MOO 
oo  o 
OOM 
•  •  • 

OOO 

••awMM  O  O 
MM  «•••  M  M 

s 

«|  MMMOmM  OM  M  MMO  M«*  MMM  OMMmMmM  MM 

u  MM  mmmmoomooommmooommmmmmmm 

m OMMMMmM OO M MMOMM MMM OMM mMM M MK 

^  ooooooo  MOMooooo  oomommmomooo 

MMM«»OMOO 

MM  M 
OMO 

MM  O 

OOM 

MMO  M  M  •• 
MKMM  M  O 
MOOMO  •• 

M 

#•  A 

S  mo  o  o  mo  CO  M  O  M  OMO  MOM  OMM  ooo  M  MO 

'-i  OO  ••  ••  MM  MM  M  M  M  M  MO  ••MM  MMM  ^MM  M 

m«»««mmmmmmmmmmmmmmmmmmmmmmmm 

M  momooooomooomo  OOM  mmmmooo  oo 

M  M  tfl  M  O  O  OO 
MOOMKOmM 

95S 

SSS 

AAA 

OOO 

OM«»MMM 

M 

o  OMMMMMMMMMMOMMO^MMMMmMMMMO 

w  MMMMMfiMMMMMMriM  MMmMMMmMMMMM 

C  O  OO  M  O  MO  M  O  M  O  MO  MO  MOM  MOM  MOM  M  MO 

Ul 

••M  M 
KM  fs 

OOO 

A  A  A 

OOO 

M  K  M  M  O  O 

OMOMOM 

0. 

«l  w 

^  oo  ooooooo >0000000M9^^MMMMM 

MMOOO^kK^kMMMOOOO  OOM  MMMmMMM  MM 

a  «« 

•  A  A 

MMM 
MM  M 
MM  M 

r 


209 


Appendix  D  Block  Loading  Test  Data 


A  complete  listing  of  the  experimental  data  is  presented  in  this  appendix  for  the  block 
loading  test.  The  data  sheet  column  variables  are  as  follows;  the  applied  delta  bending 
moment  (Mjjjajj  -  M^int  in-lb.),  "a"  and  “c"  crack  lengths,  a/c,  a/t,  accumulated  fatigue 
cycles,  AK^j,  K^.  plastic  zcne  size  at  location  C,  AK^,  plastic  zone  size  at 

location  A,  cPQp,  bP^^p  and  aPQp  applied  load  (lb.).  Negative  one  (-1.0)  indicates  that 
closure  load  data  were  not  recorded  at  that  crack  size.  Blank  lines  in  the  data  separate 
categories  of  precracking,  type  I,  type  II  and  type  III  cracks. 
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/^pendix  E  Precracking  Load  Cycle  Determination 


An  iterative  computer  code  was  generated  to  define  the  precracking  load  cycle 
needed  to  meet  the  ASTM  recommendations  for  thru-thickness  cracks,  plus  additional 
criteria  established  for  this  investigation  of  surfoce  flaws.  The  computational  procedure 
iterates  backward  from  the  desired  stress  intensity  factor  level  where  the  experiment  test  is  to 
begin  using  plastic  zone  size  as  the  major  load-shedding  criteria.  That  is,  given  a  desired 
beginning  stress  intensity  factor,  and  other  test  variables  discussed  in  Chapter  IV,  the 
program  calculates  backwards  to  deflne  an  appropriate  beginning  load  and  load-shedding 
cycle.  The  largest  load-shedding  step-size  of  20%  (recommended  by  ASTM)  is  initially 
assumed,  and  then  reduced  as  required  to  meet  all  other  criteria.  The  iterative  procedure 
continues  until  a  solution  is  reached,  or  it  is  determined  there  is  no  solution  within  the 
bounds  of  the  established  criteria. 

The  two  fundamental  criteria  established  by  the  ASTM-E647  (39)  committee  to 
define  a  specimen  piecracking  load  cycle  are; 

1.  Load  reduction  levels  should  be  no  greater  than  20%  of  the  applied  load. 

2.  The  crack  growth  before  load-shedding  should  be  at  least  three  times  the  plastic 
zone  size  established  by  the  prior  load  level. 

In  addition,  a  third  criterion  was  established  for  this  investigation  to  ensure  the  plastic 
zone  size  was  decreasing  with  load-shedding  steps.  That  is,  the  stress  intensity  factor  after 
each  load-shedding  step  shall  be  less  than  the  prior  stress  intensity  level. 

These  three  criteria  were  programmed  into  the  computer  code  as  follows: 

-  Step  #1:  a  load-shed  &ctor  (LSF)  is  derined  to  ensure  load-shedding  is  less 
than  20%  of  the  applied  load.  That  is: 

LSF  s  0.8  *  M/M'  (Bl) 
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where  M  is  the  new  applied  moment  to  the  specimen,  and  M*  the  applied  moment  before 
load-shedding.  Once  a  surface  crack  length  of  size  "c”  is  assumed,  and  a  desired  initial  stress 
intensity  factor  defined,  the  program  predicts  the  load  level  required  to  fatigue  the  crack  to  a 
size  6  times  the  plastic  zone  diameter.  That  is,  for  a  desired  final  K  value  (where  K  is  the 
value  when  the  fatigue  crack  growth  experimental  test  begins  after  precracking),  and  an 
assumed  crack  length  "c"  (typically  defined  as  a  measurable  value  which  is  consistent  with 
available  instrumentation  accuracy),  the  program  calculates  the  load-shedding  steps  by 
reducing  the  LSF  and  checking  to  ensure  the  constraints  are  satisfied.  In  effect,  the 
procedure  is  a  backward-marching  iterative  approach  for  meeting  the  required  criteria. 

-  Step  #2:  the  program  utilizes  the  Irwin  plastic  zone  size  prediction  for  plane 
stress  to  meet  the  ASTM  crack  growth  criteria  of  at  least  3  times  the  plastic  zone  size.  That 
is: 


Ip-  (K/OjjZ/n  (B.1) 

The  crack  growth  length  (Ac)  along  the  surface  is  defined  as: 

Ac  =  c  -  c'  (B.2) 

Letting  the  crack  stress  intensity  factor  be  K'  (where  the  prime  mark  represents  the  K  value 
before  load-shedding): 

3r’p  =  3(KVays)2/n 

K'  o’bH'jraVQ'  B'(a’,c',4)  (B.3) 

where  4=  0  degrees  and  of*]}  is  the  bending  stress  such  that: 


©•b  =  6M'/bt2  =  6P’  AU4bt2 


(B.4) 


To  meet  criterion  #3,  a  restriction  was  set  on  the  iteration  procedure  to  ensure  the 
plastic  zone  size  after  load-shedding  was  smaller  than  the  plastic  zone  size  before  load¬ 
shedding.  The  equations  to  develop  this  criteria  are: 

c'  =  c  -  3*K'2/ji  oys^  (B.5) 

c'  =  c  -  3a'M'2  (6H'BVoysbt2)2/Q'  (B.6) 

c/(1+3*36(M'H'B')2  (aVc’)/(Q'bt2oys)2)  (B.7) 

or  using  the  LSF  value: 

c’=  c/(l+3(aVc')/Q’*(6H'BV(aysbt2)2(M/LSF)2)  (B.8) 

The  plane  strain  plastic  zone  size  was  also  checked  by  using  the  Irwin  relationship  of: 


rp  =  (K/Oys)2/(3*Jt) 

(B.9) 

c'  =  (c  -  K'2/ ji  Oys^) 

(B.IO) 

=  c/(l+(6H'B'M')2/Q'(bt2ays)2) 

(B.ll) 

K  =  K'*LSF 

(B.12) 

To  initiate  the  precracking  analysis,  it  was  necessary  to  assume  an  a/c  value.  This 
value  was  established  by  reviewing  surface  flaw  experimental  data  on  PMMA  and  metals 
(from  published  literature)  to  determine  crack  growth  aspect  ratio  patterns.  The  review 
indicated  that  for  crack  where  "c"  <  0.1  in,  the  crack  aspect  ratio  (a/c)  was  in  the  range  of 
0.82  <  a/c  <  1.0.  A  sensitivity  study  to  deHne  the  range  of  K  for  this  range  of  a/c  was  then 
conducted,  showing  there  was  little  difference  between  a  value  of  K  for  a/c  between  0.9  and 
0.8.  Hierefore,  an  average  value  of  0.85  was  selected  for  the  load-shedding  iterative 
calculations. 
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^or  the  research  described  in  this  report,  the  beginning  stress  intensity  factor  was  set 
at  600  p8i*in^^  (659  KPa*M^'*^)  for  all  specimen  tests.  The  crack  growth  prediction  was 
initiated  by  assuming  an  LSF  of  20%,  and  an  a/c  value  of  0.85.  The  ASTM  crack  growth 
recommendation  of  3  times  the  plastic  zone  size  for  thru-thickness  cracks  was  initially  set  at 
6  to  provide  a  safety  factor  for  surface  flaws.  However,  if  the  program  could  not  converge 
using  a  factor  of  6,  the  crack  growth  requirement  was  reduced  by  increments  of  one  plastic 
zone  diameter  until  a  factor  of  3  plastic  zone  diameter  was  reached.  If  the  program  was 
unable  to  converge  by  meeting  all  the  above  criteria,  a  statement  was  generated  indicating 
that  a  precracking  load  cycle  could  not  be  generated  for  the  established  criteria. 


